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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

C-band microwave scattering from small balsam fir

An experiment to examine the C-band backscattering characteristics of conifer trees was conducted using a truck-mounted scatterometer. Small (1 m tall) balsam fir (Abies balsamea) were arranged at various equidistant spacings on a platform to present canopies of varying density to the radar. C-band backscattering measurements of a range of canopy densities were acquired under different polarizations and incidence angles. The measured backscattering coefficient from the tree canopies increased with increasing biomass, but approached a maximum at a LAI of 2.5 and fresh biomass of 3.3 kg/sq m. A backscatter model was implemented using measured canopy attributes and showed close agreement with scatterometer measurements over the range of canopy densities. Model results indicated that branches were the prime scatterers of the radar while needles were found to only slightly attenuate the radar signal.

Ranson, K. J.↗
DOE OSTI2017

Evolution of Carbon Clusters in the Detonation Products of the Triaminotrinitrobenzene (TATB)-Based Explosive PBX 9502

Here, the detonation of carbon-rich high explosives yields solid carbon as a major constituent of the product mixture and, depending on the thermodynamic conditions behind the shock front, a variety of carbon allotropes and morphologies may form and evolve. We applied time-resolved small angle x-ray scattering (TR-SAXS) to investigate the dynamics of carbon clustering during detonation of PBX 9502, an explosive composed of triaminotrinitrobenzene (TATB) and 5 wt% fluoropolymer binder. Solid carbon formation was probed from 0.1 to 2.0 μs behind the detonation front and revealed rapid carbon cluster growth which reached a maximum after ~200 ns. The late-time carbon clusters had a radius of gyration of 3.3 nm which is consistent with 8.4 nm diameter spherical particles and matched particle sizes of recovered products. Simulations using a clustering kinetics model were found to be in good agreement with the experimental measurements of cluster growth when invoking a freeze-out temperature, and temporal shift associated with the initial precipitation of solid carbon. Product densities from reactive flow models were compared to the electron density contrast obtained from TR-SAXS and used to approximate the carbon cluster composition as a mixture of 20% highly ordered (diamond-like) and 80% disordered carbon forms, which will inform future product equation of state models for solid carbon in PBX 9502 detonation product mixtures.

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

A mixing timescale model for TPDF simulations of turbulent premixed flames

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priori assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive scalar mixing across different flame regimes are appropriately accounted for.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

PDF approach for compressible turbulent reacting flows

The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.

Hsu, A. T.↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2015

Upper critical field and (non)-superconductivity of magnetars

We construct equilibrium models of compact stars using a realistic equation of state and obtain the density range occupied by the proton superconductor in strong B-fields. We do so by combining the density profiles of our models with microscopic calculations of proton pairing gaps and the critical unpairing field H{sub c2} above which the proton type-II superconductivity is destroyed. We find that magnetars with interior homogeneous field within the range 0.1 ≤ B{sub 16} ≤ 2, where B{sub 16} = B/10{sup 16} G, are partially superconducting, whereas those with B{sub 16} > 2 are void of superconductivity. We briefly discuss the neutrino emissivity and superfluid dynamics of magnetars in the light of their (non)-superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

The effective field theory of K-mouflage

We describe K-mouflage models of modified gravity using the effective field theory of dark energy. We show how the Lagrangian density K defining the K-mouflage models appears in the effective field theory framework, at both the exact fully nonlinear level and at the quadratic order of the effective action. We find that K-mouflage scenarios only generate the operator (δg{sup 00}{sub (u)}){sup n} at each order n. We also reverse engineer K-mouflage models by reconstructing the whole effective field theory, and the full cosmological behaviour, from two functions of the Jordan-frame scale factor in a tomographic manner. This parameterisation is directly related to the implementation of the K-mouflage screening mechanism: screening occurs when K' is large in a dense environment such as the deep matter and radiation eras. In this way, K-mouflage can be easily implemented as a calculable subclass of models described by the effective field theory of dark energy which could be probed by future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

GREEN PEA GALAXIES REVEAL SECRETS OF Lyα ESCAPE

We analyze archival Lyα spectra of 12 “Green Pea” galaxies observed with the Hubble Space Telescope, model their Lyα profiles with radiative transfer models, and explore the dependence of the Lyα escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies with [O iii] λ5007 equivalent widths (EWs) of hundreds of Å. All 12 Green Pea galaxies in our sample show Lyα lines in emission, with an Lyα EW distribution similar to high-redshift Lyα emitters. Combining the optical and UV spectra of Green Pea galaxies, we estimate their Lyα escape fractions and find correlations between Lyα escape fraction and kinematic features of Lyα profiles. The escape fraction of Lyα in these galaxies ranges from 1.4% to 67%. We also find that the Lyα escape fraction depends strongly on metallicity and moderately on dust extinction. We compare their high-quality Lyα profiles with single H i shell radiative transfer models and find that the Lyα escape fraction anticorrelates with the derived H i column densities. Single-shell models fit most Lyα profiles well, but not the ones with the highest escape fractions of Lyα. Our results suggest that low H i column density and low metallicity are essential for Lyα escape and make a galaxy an Lyα emitter.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Efficient Calibration of Expensive Computational Models

Accounting for uncertainty when calibrating expensive computational models is a common challenge faced by scientists and engineers. Often Bayesian techniques are adopted to estimate a probability density function over the model parameters given noisy empirical data. The methods used to perform this type of probabilistic calibration are computationally prohibitive in that they require a large number of evaluations of the expensive model. In these cases, surrogate modeling -- that is, using a fast-to-evaluate, lower fidelity stand-in for the original computational model -- may be the only option to alleviate this computational burden. However, the upfront cost of generating training data to build a surrogate model can itself be expensive. As such, it is important to be judicious when selecting training points at which the full-fidelity model is evaluated. Here, an active learning approach is proposed that enables efficient selection of training points using approximate samples of the calibrated parameter probability density function. In this way, the training points can be concentrated in regions where the calibration algorithm requires high model accuracy.

active learning↗
DOE OSTI2018

Wide ranging equation of state with Tartarus: A hybrid Green’s function/orbital based average atom code

Average atom models are widely used to make equation of state tables and for calculating other properties of materials over a wide range of conditions, from zero temperature isolated atom to fully ionized free electron gases. The numerical challenge of making these density functional theory based models work for any temperature, density or nuclear species is formidable. Here we present in detail a hybrid Green’s function/orbital based approach that has proved to be stable and accurate for wide ranging conditions. Algorithmic strategies are discussed. In particular the decomposition of the electron density into numerically advantageous parts is presented and a robust and rapid self consistent field method based on a quasi-Newton algorithm is given. Example application to the equation of state of lutetium (Z = 71) is explored in detail, including the effect of relativity, finite temperature exchange and correlation, and a comparison to a less approximate method. Furthermore, the hybrid scheme is found to be numerically stable and accurate for lutetium over at least 6 orders of magnitude in density and 5 orders of magnitude in temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1986

Quasi-periodic oscillations as a probe of galactic-bulge X-ray sources

Quasi-periodic oscillations (QPOs) with frequencies in the range 5-50 Hz have recently been discovered in X-rays from two bright galactic-bulge sources and from Sco X-1. A shot-noise model that provides a mathematical framework for analyzing the constraints imposed on physical models by the observed power density spectra is described. A variety of physical models are then examined briefly. One of the most promising, the beat-frequency-modulated accretion model, is adopted as a working hypothesis to consider the clues it may provide concerning the galactic-bulge sources. Finally, some key observations are noted.

Lamb, F. K.↗
NASA NTRS1988

Sinker tectonics - An approach to the surface of Miranda

Two of the proposed explanations for the coronae seen on Miranda involve mantle convection driven by density anomalies. In the sinker model, the coronae result from late-accreting large silicate bodies slowly sinking through an icy mantle toward the body's center; in the riser model, they result from a compositionally produced, low-density, rising diapir. The present study determines the surface stresses induced by such density anomalies and the expected surface expressions. The results are in good agreement with the predictions of the sinker model.

Janes, D. M.↗
DOE OSTI2020

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy–Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn–Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree–Fock and Bethe ansatz local density approximation) to show their successes and limitations. Here, we also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Nuclear ground-state masses and deformations: FRDM(2012)

We tabulate the atomic mass excesses and binding energies, ground-state shell-plus-pairing corrections, ground-state microscopic corrections, and nuclear ground-state deformations of 9318 nuclei ranging from {sup 16}O to A=339. The calculations are based on the finite-range droplet macroscopic and the folded-Yukawa single-particle microscopic nuclear-structure models, which are completely specified. Relative to our FRDM(1992) mass table in Möller et al. (1995), the results are obtained in the same model, but with considerably improved treatment of deformation and fewer of the approximations that were necessary earlier, due to limitations in computer power. The more accurate execution of the model and the more extensive and more accurate experimental mass data base now available allow us to determine one additional macroscopic-model parameter, the density-symmetry coefficient L, which was not varied in the previous calculation, but set to zero. Because we now realize that the FRDM is inaccurate for some highly deformed shapes occurring in fission, because some effects are derived in terms of perturbations around a sphere, we only adjust its macroscopic parameters to ground-state masses. The values of ten constants are determined directly from an optimization to fit ground-state masses of 2149 nuclei ranging from {sup 16}O to {sub 106}{sup 265}Sg and {sub 108}{sup 264}Hs. The error of the mass model is 0.5595 MeV for the entire region of nuclei included in the adjustment, but is only 0.3549 MeV for the region N≥65. We also provide masses in the FRLDM, which in the more accurate treatments now has an error of 0.6618 MeV, with 0.5181 MeV for nuclei with N≥65, both somewhat larger than in the FRDM. But in contrast to the FRDM, it is suitable for studies of fission and has been extensively so applied elsewhere, with FRLDM(2002) constants. The FRLDM(2012) fits 31 fission-barrier heights from {sup 70}Se to {sup 252}Cf with a root-mean-square deviation of 1.052 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

Global Characteristics of Porosity and Density Stratification Within the Lunar Crust from GRAIL Gravity and Lunar Orbiter Laser Altimeter Topography Data

The Gravity Recovery and Interior Laboratory (GRAIL) mission is providing unprecedentedly high-resolution gravity data. The gravity signal in relation to topography decreases from 100 km to 30 km wavelength, equivalent to a uniform crustal density of 2450 kg/cu m that is 100 kg/cu m smaller than the density required at 100 km. To explain such frequency-dependent behavior, we introduce rock compaction models under lithostatic pressure that yield radially stratified porosity (and thus density) and examine the depth extent of porosity. Our modeling and analysis support the assertion that the crustal density must vary from surface to deep crust by up to 500 kg/cu m. We found that the surface density of mega regolith is around 2400 kg/cu m with an initial porosity of 10-20%, and this porosity is eliminated at 10-20 km depth due to lithostatic overburden pressure. Our stratified density models provide improved fits to both GRAIL primary and extended mission data.

Density↗