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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The electric fields and “lightning jets” of the Sun and solar wind

A model of electric energization of the solar wind and corona is developed, including electro-magnetic particle effects precluded by traditional magneto-hydro (MHD) assumptions. Using standard 1D radial solar models for particle density and temperature, the core gravito-electric field is calculated; and the range of possible photospheric photo-electric fields is estimated. The extant DC field apparently arises from about 460 C of charge displacement, mainly caused by the immense solar energy flux pushing electrons outward. Energetically, this electric field can accelerate surface protons out of the 2 keV gravity well and up to the 4 keV energies observed in the fast solar wind. The electrical energy is released in pervasive, persistent “proton lightning jets,” which are proton beams, charge-neutralized by co-propagating electrons. The jets are formed by pinched “avalanche breakdown” of the weakly ionized photosphere, probably initiated on the down-welling edges of solar surface granulations. These energetic jets will glow as discrete filamentary surface spicules and will be observed in reflected solar light as the diffuse K-corona. Significantly, the total charge displacement and observed fast wind speed are quantitatively predicted by a novel “virial limit”: the (positive) electrical potential energy at r = 0 is limited in magnitude to the 10 keV gravitational well at r = 0.

Physics↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Comparing ICME simulations with scaled laboratory experiment

In stellar physics and astrophysics, numerical simulations and laboratory experiments are often compared to observational data to support their representation of the real world. However, there is also merit in comparing numerical simulations to properly scaled experiments, especially when the experiment and the simulation are both emulating the solar phenomena. Confirming the credibility of scaled experiments and their scaling with well-validated models is important to expand our knowledge of the associated physical phenomena. This is significant because experiments and simulations can be performed frequently, whereas observations may be limited by location, field of view, and missing data. In this work, we use the Alfvén Wave Solar atmosphere Model, a well-validated magnetohydrodynamic model, to simulate an interplanetary coronal mass ejection (ICME) and compare it to an experiment which provides a scaled analog to a physical ICME. The experiment was performed on the Big Red Ball facility and scaled using dimensionless parameters such as plasma β and magnetosonic Mach number to reproduce the main structure of an ICME. We compare the model-simulated temperature, density, and magnetic field to those from the experiment, as well as the scaling parameters used in the experiment, to those calculated from the simulation. This comparison is performed to further justify the scaling arguments made by the experiment. Additionally, the comparison would lead to the development of stronger scaling arguments for future experiments.

Bryant, K. [University of Michigan, Ann Arbor, MI ↗
DOE OSTI2021

Nonlinear nanoelectrodynamics of a Weyl metal

Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe 4 . As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their “symmetry fingerprints” in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe 4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Here, nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electronic and structural properties of single-crystal Jahn–Teller active Co 1+x Mn 2–x O 4 thin films

Recent investigations on spinel CoMn 2 O 4 have shown its potential for applications in water splitting and fuel cell technologies as it exhibits strong catalytic behavior through oxygen reduction reactivity. To further understand this material, we report for the first time the synthesis of single-crystalline Co 1+x Mn 2-x O 4 thin films using molecular beam epitaxy. By varying sample composition, we establish links between cation stoichiometry and material properties using in-situ x-ray photoelectron spectroscopy, x-ray diffraction, scanning transmission electron microscopy, x-ray absorption spectroscopy, and spectroscopic ellipsometry. Our results indicate that excess Co ions occupy interstitial tetrahedral sites up to a certain stoichiometry, then are substitutional for octahedrally-coordinated Mn beyond this threshold. Here, we compare these results with density functional theory models of stoichiometric CoMn 2 O 4 to understand how the Jahn-Teller distortion and hybridization in Mn-O bonds impact the ability to hole dope the material with excess Co. The findings provide important insights into CoMn 2 O 4 and related spinel oxides for future application as inexpensive oxygen reduction reaction catalysts.

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

Methanol carbonylation to acetaldehyde on Au particles supported by single-layer MoS 2 grown on silica

Homogenous single-layer MoS 2 films coated with sub-single layer amounts of gold are found to isolate the reaction of methanol with carbon monoxide, the fundamental step toward higher alcohols, from an array of possible surface reactions. Active surfaces were prepared from homogenous single-layer MoS 2 films coated with sub-single layer amounts of gold. These gold atoms formed clusters on the MoS2 surface. A gas mixture of carbon monoxide (CO) and methanol (CH 3 OH) was partially converted to acetaldehyde (CH 3 CHO) under mild process conditions (308 kPa and 393 K). This carbonylation of methanol to a C 2 species is a critical step toward the formation of higher alcohols. Density functional theory modeling of critical steps of the catalytic process identify a viable reaction pathway. Imaging and spectroscopic methods revealed that the single layer of MoS 2 facilitated formation of nanoscale gold islands, which appear to sinter through Ostwald ripening. Here, the formation of acetaldehyde by the catalytic carbonylation of methanol over supported gold clusters is an important step toward realizing controlled production of useful molecules from low carbon-count precursors.

2D catalyst↗
DOE OSTI2016

Comparison of weak lensing by NFW and Einasto halos and systematic errors

Recent N-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens density profile. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos (∼> 10{sup 15}M{sub ⊙}/h) can be over- and underestimated by 0∼ 1 per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, as observed in multi-wavelength observations of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be set apart either with deep observations of exceptionally massive structures (∼> 2×10{sup 15}M{sub ⊙}/h) or by stacking the shear profiles of thousands of group-sized lenses (∼> 10{sup 14}M{sub ⊙}/h)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Density and temperature profiles after low-Z and high-Z shattered pellet injections on DIII-D

Abstract In this work we utilize the recently upgraded Thomson scattering diagnostic to resolve density and temperature plasma profiles after pure deuterium and mixed neon/deuterium shattered pellet injections (SPIs) on DIII-D. This allows us to study individual components of the staggered scheme proposed for disruption mitigation on ITER, consisting of a low-Z material SPI followed by a delayed high-Z SPI. Obtained spatio-temporal density profiles exhibit very different dynamics after dominantly neon and pure deuterium SPIs. The neon SPI causes a fast radiative plasma collapse in a few milliseconds and results in almost flat density profile once the impurity mixes with the plasma during and after the thermal quench (TQ). The deuterium SPI leads to a disruption delayed by ten and more milliseconds, but very limited core fueling can be observed before the disruption. Even during and after the TQ, the edge deuterium density significantly exceeds the core density. 1D transport modeling suggests that this poor core fueling can be explained by strong outward grad-B-induced drift of the injected deuterium. Preliminary simulations show that larger pellet shards and greater injected quantity can be used to improve the penetration of the low-Z material into the core. These results call for optimization and further evaluation of the staggered SPI on ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

SDSS IV MaNGA - star-formation driven biconical outflows in face-on galaxies

ABSTRACT We find 132 face-on and low inclination galaxies with central star formation driven biconical gas outflows (FSFB) in the SDSS MaNGA (Mapping Nearby Galaxies at APO) survey. The FSFB galaxies show either double peaked or broadened emission line profiles at their centres. The peak and maximum outflow velocities are 58 and 212 km s−1, respectively. The gas velocity dispersion reveals a mild dependence on the central star formation surface density compatible with models of gas dispersion powered by the Jeans instability in gas clumps or by gas turbulence dissipation. We estimate the gas outflow rate and conclude that the central gas depletion time does not depend on galactic mass. In turn, the ratio of the gas outflow rate to the gas consumption rate by the star formation is low in massive galaxies and high in low mass objects, while the star formation is a more rapid process of the gas consumption. We compare properties of the FSFB galaxies with a control sample of 375 comparison galaxies and find that the FSFB objects have high central concentration of star formation and also younger central stellar population with respect to their periphery. We analysed the environment of the galaxies and identified nearby satellites and elements of low surface brightness structure. We see that many tidal-enhanced features that can be assigned to early and intermediate stages of galactic interaction are much more frequent in the FSFB galaxies with respect to the comparison sample. We conclude that the gas should be replenished via the accretion from small satellites.

Bizyaev, D. (ORCID:000000023601133X)↗
DOE OSTI2023

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

169 Tm ( n , γ ) cross section and statistical decay properties from measurements at the DANCE facility

Background: Radiative neutron capture on thulium, which is a monoisotopic element, plays a role in different applications such as nuclear astrophysics or nuclear burning environments. Considerable discrepancies—reaching 20%—exist between evaluations in the unresolved-resonance region. Furthermore, experimental data on statistical 𝛾 decay in odd-odd rare-earth nuclei is scarce. There are still open questions about the systematics of the so-called scissors mode in the 𝑀⁢1 photon strength function, especially in odd-odd nuclei. Purpose: This work is focused on two main topics—deriving experimental 169 Tm ⁢(𝑛,𝛾) cross section and studying statistical 𝛾 decay of 170 Tm, in particular properties of the scissors mode. Methods: The capture experiments to obtain experimental cross section were performed at the Los Alamos Neutron Science Center using the time-of-flight technique and employing the Detector for Advanced Neutron Capture Experiments. Measured coincident 𝛾-ray spectra were also compared with statistical simulations using the dicebox code to test different models of level density and photon strength functions. Results: The capture cross section was determined from 1.8 eV to 0.97 MeV, the broadest neutron-energy range ever measured for this isotope. Several new resonances have been observed. The statistical 𝛾 decay of 170 Tm cannot be reproduced without a scissors mode resonance centered at ≈ 3.3MeV. Conclusions: The measured cross section in the unresolved-resonance region is generally lower than the latest evaluations. The derived 169 Tm 𝑠-process abundance is expected to increase by a factor of 1.26, while the changes of the abundances of elements heavier than 169 Tm are in the order of 0.2%. The scissors mode properties in 170 Tm are similar to those deduced in previous analyses of neighboring nuclei 168 Er and 166 Ho .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Photoionization and photofragmentation of singly charged positive and negative Sc 3 N @ C 80 endohedral fullerene ions

Photoprocesses of the endohedral fullerene ions $\mathrm{Sc_3N@C_{80}^{+}}$ and $\mathrm{Sc_3N@C_{80}-}$ in the gas phase have been investigated in the photon energy ranges 30-50 eV and 280-420 eV. Single and double ionization as well as single ionization accompanied by the release of a C 2 dimer were observed as a function of the photon energy for the positive parent ion and double detachment was measured for the negative parent ion. The emphasis of the experiments was on the specific effects of the encapsulated trimetallic nitride cluster Sc 3 N on the observed reactions. Clear evidence of photoexcitation near the Sc L edge is obtained with the dominating contributions visible in the one- and two-electron-removal channels. K-vacancy production in the encapsulated central nitrogen atom is seen in the single ionization of $\mathrm{Sc_3N@C_{80}^{+}}$ but is much less pronounced in the photoionization-with-fragmentation channel. Comparison of the cross sections near the carbon K edge with the corresponding channels measured previously in the photoionization of $\mathrm{Lu_3N@C_{80}^+}$ reveal strong similarities. Previously predicted sharp resonance features in the ionization of $\mathrm{Sc_3N@C_{80}^{+}}$ ions below the Sc M edge are not confirmed. Finally, the experiments are accompanied by quantum-chemistry calculations in the Hartree-Fock approximation and by model calculations employing density functional theory (DFT).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Phase diagram of boron carbide with variable carbon composition

Boron carbide exhibits intrinsic substitutional disorder over a broad composition range. The structure consists of 12-atom icosahedra placed at the vertices of a rhombohedral lattice, together with a 3-atom chain along the threefold axis. In the high-carbon limit, one or two carbon atoms can replace boron atoms on the icosahedra while the chains are primarily of type C-B-C. We fit an interatomic pair interaction model to density-functional-theory total energies to investigate the substitutional carbon disorder. Monte Carlo simulations with sampling improved by replica exchange and augmented by two-dimensional multiple histogram analysis predict three phases. The low-temperature, high-carbon-composition monoclinic Cm structure disorders through a pair of phase transitions, first via an Ising-like transition to a monoclinic centrosymmetric state with space group C2/m, then via a first-order three-state Potts-like transition to the experimentally observed rhombohedral R$\bar3$ m symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Radius and equation of state constraints from massive neutron stars and GW190814

Motivated by the unknown nature of the 2.50–2.67M⊙ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. Here we demonstrate that a modification of the highly uncertain suprasaturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of 2.6M⊙ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of ≈ 2.0M⊙ neutron stars (but not necessarily the radii of 1.4 M⊙ neutron stars), which further motivates future NICER radius measurements of PSR J1614 – 2230 and PSR J0740+6620.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Impact of alternative transmission coeffcient parameterizations on Hauser-Feshbach theory

Here, we investigate different formulations of the transmission coefficient T c , including the form implied by Moldauer’s “sum rule for resonance reactions” and the Moldauer-Simonius form. Within these different formulations, we compute the neutron transmission coefficients in the resolved and unresolved resonance regions, allowing a direct comparison with the transmission coefficients computed using an optical model potential. For nuclei for which there are no measured resonances, these approaches allow one to predict the average neutron resonance parameters directly from the optical model and level densities. Some of the approaches are valid in both the strong and weak coupling limits (i.e., any value of the average width and mean level spacing). Finally, both the Moldauer-Simonius and Moldauer’s Sum Rule forms approaches suggest that superradiance, that is, the quantum chaotic enhancement of certain channels, may be a common phenomena in nuclear collisions. Our results suggest why superradiance has been previously overlooked. We apply our approach to neutron reactions on the closed shell 90 Zr nucleus and the mid-shell 197 Au nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗