Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density model”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Guidelines to Achieving High Selectivity for the Hydrogenation of α,β-Unsaturated Aldehydes with Bimetallic and Dilute Alloy Catalysts: A Review

Selective hydrogenation of unsaturated aldehydes to unsaturated alcohols is a challenging class of reactions, yielding valuable intermediates for the production of pharmaceuticals, perfumes, and flavorings. On monometallic heterogeneous catalysts, the formation of the unsaturated alcohols is thermodynamically disfavored over the saturated aldehydes. Hence, new catalysts are required to achieve the desired selectivity. In this work, the literature of three major research areas in catalysis is integrated as a step toward establishing the guidelines for enhancing the selectivity: reactor studies of complex catalyst materials at operating temperature and pressure, surface science studies of crystalline surfaces in ultrahigh vacuum, and first-principles modeling using density functional theory calculations. Aggregate analysis shows that bimetallic and dilute alloy catalysts significantly enhance the selectivity to the unsaturated alcohols compared to monometallic catalysts. This comprehensive review focuses primarily on the role of different metal surfaces as well as the factors that promote the adsorption of the unsaturated aldehyde via its C=O bond, most notably by electronic modification of the surface and formation of the electrophilic sites. Furthermore, challenges, gaps, and opportunities are identified to advance the rational design of efficient catalysts for this class of reactions, including the need for systematic studies of catalytic processes, theoretical modeling of complex materials, and model studies under ambient pressure and temperature.

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

Weak lensing mass-richness relation of redMaPPer clusters in LSST DESC DC2 simulations

Cluster scaling relations are key ingredients in cluster abundance-based cosmological studies. In optical cluster cosmology, where clusters are detected through their richness, cluster-weak gravitational lensing has proven to be a powerful tool to constrain the cluster mass-richness relation. This work is conducted as part of the Dark Energy Science Collaboration (DESC), which aims to analyze the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory, starting in 2026. Cluster properties inferred from weak lensing, such as mass, suffer from several sources of bias. In this paper, we aim to test the impact of modeling choices and observational systematics in cluster lensing on the inference of the mass-richness relation. We constrained the mass-richness relation of 3600 clusters detected by the redMaPPer algorithm in the cosmoDC2 extragalactic mock catalog of the LSST DESC DC2 simulation, covering 440 deg 2 , using number count measurements and either stacked weak lensing profiles or mean cluster masses in several intervals of richness (20 ≤ λ ≤ 200) and redshift (0.2 ≤ z ≤ 1). We provide the first constraints on the redMaPPer cluster mass-richness relation detected in cosmoDC2. We find that for an LSST-like source galaxy density, our constraints are robust to changes in the concentration-mass relation, as well as the dark matter density profile modeling choices, when source redshifts and shapes are perfectly known. We find that photometric redshift uncertainties can introduce bias at the 1σ level, which could be mitigated by an overall correction factor fitted jointly with the scaling parameters. We find that including positive shear-richness covariance in the fit shifts the results by up to 0.5σ. Our constraints also offer a fair comparison to a fiducial mass-richness relation, obtained from matching cosmoDC2 halo masses to redMaPPer-detected cluster richness results.

galaxy clusters↗
DOE OSTI2021

Granular Flow Constitutive Model Add-ons For Abaqus

The software is a list of Abaqus User MATerial subroutines (VUMAT) for modeling granular flow physics. It includes (1) density dependent Mohr-Coulomb model, (2) density dependent Drager-Prager\Cap model, (3)Gudehus-Bauer hypoplastic model, and (4) critical state-based NorSand model.

Wencheng, Wencheng↗
DOE OSTI2016

Lithium doping on covalent organic framework-320 for enhancing hydrogen storage at ambient temperature

Density Functional Theory (DFT) combines with grand canonical Monte Carlo (GCMC) simulations are performed to explore the effect of Li doping on the hydrogen storage capability of COF-320. The results show that the interaction energy between the H{sub 2} and the Li-doped COF-320 is about three times higher than that of pristine COF-320. GCMC simulations are employed to study the hydrogen uptake of Li-doped COF-320 at ambient temperature, further confirm that the lithium doping can improve the hydrogen uptake at ambient temperature. Our results demonstrate that Li-doped COFs have good potential in the field of hydrogen storage. - Graphical abstract: Fig. 1. The optimized cluster model used here to represent the COF-320 and possible adsorption sites (A, B, C) for adsorption of metals in the COF-320. The dangling bonds are terminated by H atoms. C, H, and N atoms are shown as gray, white, and blue colors, respectively. Fig. 2. The adsorption isotherm of H{sub 2} in the pristine and Li-doped COF-320 at 298 K. - Highlights: • The binding sites of single and two lithium atoms in COF-320 were studied. • The interaction energy between the H{sub 2} and the Li-doped COF-320 is about three times higher than that of pristine COF-320. • H{sub 2} uptakes on the Li-doped COFs obtain significant improvement at ambient temperature. • Lithium-doping is a successful strategy for improving hydrogen uptake.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗
NASA NTRS1974

Neutral and ion exosphere models for lunar hydrogen and helium

A general neutral exosphere model, which includes density and temperature variations at the exobase, is applied to the moon to obtain surface and radial density distributions for H, H2, and He. It is assumed that the source for these constituents derives from accretion of solar wind ions. The surface distributions are determined by requiring that the sum of the neutral and solar wind ion fluxes for a given constituent vanish at all points on the surface. On this basis, maximum dayside surface densities for H, H2, and He and maximum nightside surface densities for H, H2, and He are obtained that are consistent with either measured values or upper limits. In addition, model ion density distributions for H2(+) and He(+) are constructed. This ion exosphere is produced by ionization of the neutral exosphere in the solar wind, which efficiently sweeps the ions past or onto the lunar surface. Saturated H2(+) and He(+) densities ranging from about .001 to .015 per cu cm and .00003 to .0004 per cu cm over 1.5 to 3 selenocentric radii on the dayside, respectively.

Hartle, R. E.↗
DOE OSTI2020

Establishing the Maximum Collectivity in Highly Deformed N = Z Nuclei

The lifetimes of the first excited 2 + states in the N = Z nuclei 80 Zr, 78 Y, and 76 Sr have been measured using the γ-ray line shape method following population via nucleon-knockout reactions from intermediate-energy rare-isotope beams. The extracted reduced electromagnetic transition strengths yield new information on where the collectivity is maximized and provide evidence for a significant, and as yet unexplained, odd-odd vs even-even staggering in the observed values. As a result, the experimental results are analyzed in the context of state-of-the-art nuclear density-functional model calculations.

59 ≤ A ≤ 89↗
DOE OSTI2017

Non-singular and cyclic universe from the modified GUP

In this paper, we investigate the effects of a new version of the generalized uncertainty principle (modified GUP) on the dynamics of the Universe. As the modified GUP will modify the relation between the entropy and area of the apparent horizon, it will also deform the Friedmann equations within Jacobson's approach. We explicitly find these deformed Friedmann equations governing the modified GUP-corrected dynamics of such a Universe. It is shown that the modified GUP-deformed Jacobson's approach implies an upper bound for the density of such a Universe. The Big Bang singularity can therefore also be avoided using the modified GUP-corrections to horizons' thermodynamics. In fact, we are able to analyze the pre Big Bang state of the Universe. Furthermore, the equations imply that the expansion of the Universe will come to a halt and then will immediately be followed by a contracting phase. When the equations are extrapolated beyond the maximum rate of contraction, a cyclic Universe scenario emerges.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A priori filtered chemical source term modeling for LES of high Karlovitz number premixed flames

Assumptions behind closure models for the filtered source term are studied a priori using results from DNS of turbulent n-heptane/air premixed flames at varying Karlovitz numbers. Simulations with both detailed chemistry and tabulated chemistry, as well as unity and non-unity Lewis numbers, are used to determine if finite-rate chemistry and differential diffusion effects affect the filtered chemical source terms. While the unfiltered source term shows large fluctuations, the filtered source terms from detailed chemistry and tabulated chemistry are in good agreement at sufficiently large filter widths (Δ ≳ l F ). Using the concept of optimal estimators, it is shown that a tabulation approach using the filtered progress variable and its variance can predict accurately the filtered chemical source terms. Finally, the filtered source terms from the DNS are compared to predictions from two commonly assumed sub-filter probability density function models. Both models show deviations from the filtered DNS source terms but predict accurately the mean turbulent flame speed. In conclusion, the results illustrate the potential of using simple tabulated chemistry approaches based on presumed PDFs for LES of premixed flames in the thin and distributed reaction zones regimes.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

The Fourth-Order Symmetry Energy of Finite Nuclei

The fourth-order symmetry energy E{sub sym,4}(A) of heavy nuclei is investigated based on the Skyrme energy density functional in combination with a local density approximation. Unlike some previous works, in our method, the interferences from the other energy terms are removed since it is completely isolated from the rest of energy terms. The calculated E{sub sym,4}(A) is much less than that extracted from nuclear masses. The underlying reason for the big difference is discussed. The Brueckner theory also gives a small fourth-order symmetry energy coefficient of nuclear matter, which is also different from recent conclusions with another methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

On the Relative Importance of Convection and Temperature on the Behavior of the Ionosphere in North American during January 6-12, 1997

Measurements from a network of digisondes and an incoherent scatter radar In Eastern North American For January 6-12, 1997 have been compared with the Field Line Interhemispheric Plasma (FLIP) model which now includes the effects of electric field convective. With the exception of Bermuda, the model reproduces the daytime electron density very well most of the time. As is typical behavior for winter solar minimum on magnetically undisturbed nights, the measurements at Millstone Hill show high electron temperatures before midnight followed by a rapid decay, which is accompanied by a pronounced density enhancement in the early morning hours. The FLIP model reproduces the nighttime density enhancement well, provided the model is constrained to follow the topside electron temperature and the flux tube is full. Similar density enhancements are seen at Goose Bay, Wallops Island and Bermuda. However, the peak height variation and auroral images indicate the density enhancements at Goose Bay are most likely due to particle precipitation. Contrary to previously published work we find that the nighttime density variation at Millstone Hill is driven by the temperature behavior and not the other way around. Thus, in both the data and model, the overall nighttime density is lowered and the enhancement does not occur if the temperature remains high all night. Our calculations show that convections of plasma from higher magnetic latitudes does not cause the observed density maximum but it may enhance the density maximum if over-full flux tubes are convected over the station. On the other had, convection of flux tubes with high temperatures and depleted densities may prevent the density maximum from occurring. Despite the success in modeling the nighttime density enhancements, there remain two unresolved problems. First, the measured density decays much faster than the modeled density near sunset at Millstone Hill and Goose Bay though not at lower latitude stations. Second, we cannot fully explain the large temperatures before midnight nor the sudden decay near midnight.

Richards, P. G.↗
DOE OSTI2024

Electronic structure and magnetic tendencies of trilayer La 4 Ni 3 O 10 under pressure: Structural transition, molecular orbitals, and layer differentiation

Motivated by the recent observation of superconductivity in the pressurized trilayer Ruddlesden-Popper (RP) nickelate La 4 Ni 3 O 10 , we explore its structural, electronic, and magnetic properties as a function of hydrostatic pressure from first-principles calculations. We find that an orthorhombic (monoclinic)-to-tetragonal transition under pressure takes place concomitantly with the onset of superconductivity. The electronic structure of La 4 Ni 3 O 10 can be understood using a molecular trimer basis wherein n molecular subbands arise as the d z 2 orbitals hybridize strongly along the c axis within the trilayer. The magnetic tendencies indicate that the ground state at ambient pressure is formed by nonmagnetic inner planes and stripe-ordered outer planes that are antiferromagnetically coupled along the c axis, resulting in an unusual ↑, 0, ↓ stacking that is consistent with the spin density wave model previously suggested by neutron diffraction. Such a state is destabilized at the pressure where superconductivity arises. Despite the presence of d z 2 states at the Fermi level, the d x 2 –y 2 orbitals also play a key role in the electronic structure of La4Ni3O10. Finally, this active role of the d x 2 –y 2 states in the low-energy physics of the trilayer RP nickelate, together with the distinct electronic behavior of the inner and outer planes, resembles the physics of multilayer cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

ENERGETIC NEUTRAL HYDROGEN OBSERVATIONS DEMONSTRATE THAT VOYAGER 1 IS NOT OBSERVING THE EXTRAORDINARILY STRONG INTERSTELLAR MAGNETIC FIELD

It is generally believed that Voyager 1 ( V1 ) is now in interstellar space, having crossed the heliopause at a heliocentric distance of 121.58 au in late August of 2012. Here we use recently published spectra of energetic neutral hydrogen, and the magnetic field and energetic particles directly measured by V1 to find the average pressure in the inner heliosheath (termination shock to 122 au). This pressure turns out to be surprisingly large, (3.57 ± 0.71) × 10{sup −12} dyn cm{sup −2}, and is completely dominated by pressures of pickup ions (PUIs), created in the inner heliosheath, and their suprathermal tails (43%), and PUIs and their tails that are produced upstream of the termination shock and enter the heliosheath (46%). We compute the total particle pressure in the outer heliosheath near the heliopause from distribution functions of the interstellar plasma and locally created PUIs using profiles of proton density, proton temperature, and neutral hydrogen density from model 2 in Zank et al., and find it to be at most 7.7 × 10{sup −13} dyn cm{sup −2}. Balancing pressure across the heliopause, thus requires an unusually large magnetic pressure (2.8 × 10{sup −12} dyn cm{sup −2}). The resulting strength and 1 σ uncertainty of the draped magnetic field in the outer heliosheath near the heliopause is 0.839 ± 0.106 nT. The 3 σ lower limit field strength (0.52 nT) is greater than the field of ∼0.43 ± 0.02 nT measured by V1 , implying that there is less than 1% probability that V1 is measuring the interstellar draped field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A full picture of large lepton number asymmetries of the Universe

A large lepton number asymmetry of O(0.1−1) at present Universe might not only be allowed but also necessary for consistency among cosmological data. We show that, if a sizeable lepton number asymmetry were produced before the electroweak phase transition, the requirement for not producing too much baryon number asymmetry through sphalerons processes, forces the high scale lepton number asymmetry to be larger than about 03. Therefore a mild entropy release causing O(10-100) suppression of pre-existing particle density should take place, when the background temperature of the Universe is around T = O(10{sup −2}-10{sup 2}) GeV for a large but experimentally consistent asymmetry to be present today. We also show that such a mild entropy production can be obtained by the late-time decays of the saxion, constraining the parameters of the Peccei-Quinn sector such as the mass and the vacuum expectation value of the saxion field to be m {sub φ} ∼> O(10) TeV and φ{sub 0} ∼> O(10{sup 14}) GeV, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Testing approximate predictions of displacements of cosmological dark matter halos

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at z =0,0.5,1, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods, and this is true even for mesh sizes as large as the inter-particle distance. This test sets an upper limit on the ability of these methods to reproduce the clustering of halos, for the cases when these objects are reconstructed at the object-by-object level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Improved image of intrusive bodies at Newberry Volcano, Oregon, based on 3D gravity modelling

Beneath Newberry Volcano is one of the largest geothermal heat reservoirs in the western United States and it has been extensively studied for the last 40 years. Several magmatic intrusions have been recognized at depths between 2.5 and 8 km and some of them identified as suitable targets for enhanced geothermal energy and tested during two previous EGS campaigns. These subsurface structures have been intersected by three deep wells and imaged by various geophysical methods including seismic tomography and magnetotellurics. Although three high quality gravity surveys were completed between 2006 and 2010 as part of various projects, a complete synthesis and interpretation of the gravity data has not yet been performed. Regional gravity data also exist in the vicinity of the Newberry volcano and have been added to these surveys to constitute a dataset with a total of 1418 gravity measurements. When coupled with existing geologic and geophysical data and models, this new gravity dataset provides important constraints on the depth and contours of the magmatic bodies previously identified by other methods and thus greatly contributing to facilitate any future drilling and stimulation works. Using the initial structures discovered by seismic tomography, inversion of gravity data has been performed. Shape, density values and depths of various bodies were allowed to vary and three main bodies have been identified. Densities of the middle and lower intrusive bodies (~2.6-2.7 g/cm3) are consistent with rhyolite, basalt or granites. Modeled density of the near-surface caldera body match that of a low density tephra material and the density of the shallow ring structures contained in the upper kilometer correspond to that of welded tuff or low-density rhyolites. Modeled bodies are in reality a composite of thin layers; however, average densities of the modeled gravity bodies are in good agreement with the density log obtained in one well located on the western flank (well 55-29). Final gravity data residuals show that most of the observed gravity anomalies at the surface can be explained by the modeled gravity bodies and are consistent with other site characterization information.

Cascade Range↗