Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Density”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Arrested coarsening and large density fluctuations in driven particle mixtures in two dimensions

Abstract Using molecular dynamics simulations, we study a driven, nonadditive binary mixture of spherical particles confined to move in two dimensions and immersed in an explicit solvent consisting of point particles with purely repulsive interactions. We show that, without a drive, the mixture of spherical particles phase separates and coarsens with kinetics consistent with an Ising-like conserved dynamics. Conversely, when the drive is applied, the coarsening is arrested and the system develops large density fluctuations. We show that the drive creates domains of a characteristic size which decreases with an increasing force. Furthermore, we find that these domains are anisotropic and can be oriented either parallel or perpendicular to the drive direction. Finally, we connect our findings to existing theories of strongly-driven systems, pointing out the importance of introducing the explicit solvent particles to break the Galilean invariance of the system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗
DOE OSTI2020

Collective effects and pattern formation for directional locking of disks moving through obstacle arrays

We examine directional locking effects in an assembly of disks driven through a square array of obstacles as the angle of drive rotates from 0° to 90°. For increasing disk densities, the system exhibits a series of different dynamic patterns along certain locking directions, including one-dimensional or multiple-row chain phases and density-modulated phases. For nonlocking driving directions, the disks form disordered patterns or clusters. Here, when the obstacles are small or far apart, a large number of locking phases appear; however, as the number of disks increases, the number of possible locking phases drops due to the increasing frequency of collisions between the disks and obstacles. For dense arrays or large obstacles, we find an increased clogging effect in which immobile and moving disks coexist.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Interaction of a sodium ion with the water liquid-vapor interface

Molecular dynamics results are presented for the density profile of a sodium ion near the water liquid-vapor interface at 320 K. These results are compared with the predictions of a simple dielectric model for the interaction of a monovalent ion with this interface. The interfacial region described by the model profile is too narrow and the profile decreases too abruptly near the solution interface. Thus, the simple model does not provide a satisfactory description of the molecular dynamics results for ion positions within two molecular diameters from the solution interface where appreciable ion concentrations are observed. These results suggest that surfaces associated with dielectric models of ionic processes at aqueous solution interfaces should be located at least two molecular diameters inside the liquid phase. A free energy expense of about 2 kcal/mol is required to move the ion within two molecular layers of the free water liquid-vapor interface.

NASA Center ARC↗
DOE OSTI2023

Attosecond light-field control of high-density plasmas (Final Scientific/Technical Report)

Sub-optical-cycle dynamics of dense electron bunches in relativistic-intensity laser—solid interactions (relativistic plasma mirrors) lead to the emission of high-order harmonics and attosecond light pulses. The aim of this project was to advance our understanding of the fundamental principles underlying high-energy-density physics behind this phenomenon. We achieved this goal by conducting experiments and a comprehensive theoretical analysis of harmonics driven by specially tailored light waveforms providing precise control over the sub-laser-cycle trajectories of the emitting electron bunches. The multi-color laser waveforms, generated in our experiments in a cascaded plasma mirror configuration, allowed us to manipulate (enhance and suppress) harmonic generation in relativistic laser-solid interactions by adjusting the relative phase between the colors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗
NASA NTRS2022

The Earth's Exosphere and Its Response to Space Weather

Neutral-Plasma charge exchange is a fundamental physical process that occurs ubiquitously across the universe. In geospace, charge exchange occurs in the Earth’s topside ionosphere, polar wind, plasmasphere, inner magnetosphere, and magnetosheath. Past and current space missions have profiled plasma and electromagnetic characteristics in various parts of the Earth’s magnetospheric system. However, observations of the exosphere, i.e., neutrals above 500 km altitude, are still sparse and, in some regions, non-existent, which limits our understanding of the neutral contribution to the overall dynamics of the geospace environment. Cold exospheric neutrals (< 10 eV) play an important role in the Sun-Earth interaction. Variability of exospheric density provides key information of the Earth’s atmospheric loss under dynamic space environment conditions. Various neutral species and their density variations in the polar wind can alter ion outflow patterns, modifying global magnetospheric dynamics. Exospheric neutrals also provide an energy sink for the inner magnetosphere by creating Energetic Neutral Atoms (ENAs) through charge exchange with high-energy ring current ions, which subsequently leave our geospace system unimpeded by magnetic fields. Exospheric neutrals also provide a means to observe the global interaction of the solar wind – magnetosphere, through global imaging of the system via ENAs (e.g., the TWINS and IMAGE missions) and soft X-rays (e.g., the upcoming LEXI and SMILE missions), the byproducts of neutral-plasma charge exchange. In the coming decade, we advocate that our community needs to increase our exploration of the neutral populations in the outermost reaches of the Earth’s atmosphere. It is imperative that we improve both in-situ and remote-sensing technologies for measuring key neutral species in our exosphere. We also encourage dedicated exosphere missions and to stimulate model developments of our exosphere and its interaction with the co-located magnetospheric system and neighboring Ionosphere - Thermosphere - Mesosphere system.

Hyunju Connor↗
DOE OSTI2019

Influence of Electronic Polarization on the Spectral Density

Accurate spectral densities are necessary for computing realistic exciton dynamics and nonlinear optical spectra of chromophores in condensed-phase environments, including multichromophore pigment–protein systems. However, due to the significant computational cost of computing spectral densities from first principles, requiring many thousands of excited-state calculations, most simulations of realistic systems rely on treating the environment as fixed-point charges. In this work, using a number of representative systems ranging from solvated chromophores to the photoactive yellow protein (PYP), we demonstrate that the quantum mechanical (QM) electronic polarization of the environment is key to obtaining accurate spectral densities and line shapes within the cumulant framework. We show that the QM environment can enhance or depress the coupling of fast chromophore degrees of freedom to the energy gap, altering the electronic–vibrational coupling and the resulting vibronic progressions in the absorption spectrum. In analyzing the physical origin of peaks in the spectral density, we identify vibrational modes that couple the electron and the hole as being particularly sensitive to the QM screening of the environment. For PYP, we reveal the need for careful determination of the appropriate QM region to obtain reliable spectral densities. Our findings suggest that the QM polarization of the environment can be crucial not just for excitation energies but also for electronic–vibrational coupling in complex systems with implications for the correct modeling of linear and nonlinear optical spectroscopy in the condensed phase as well as energy transfer in pigment–protein complexes.

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

Accuracy of Kohn–Sham density functional theory for warm- and hot-dense matter equation of state

We study the accuracy of Kohn–Sham density functional theory (DFT) for warm- and hot-dense matter (WDM and HDM). Specifically, considering a wide range of systems, we perform accurate ab initio molecular dynamics simulations with temperature-independent local/semilocal density functionals to determine the equations of state at compression ratios of 3x–7x and temperatures near 1 MK. We find very good agreement with path integral Monte Carlo benchmarks, while having significantly smaller error bars and smoother data, demonstrating the accuracy of DFT for the study of WDM and HDM at such conditions. In addition, using a Δ-machine learned force field scheme, we confirm that the DFT results are insensitive to the choice of exchange-correlation functional, whether local, semilocal, or nonlocal.

Suryanarayana, Phanish (ORCID:0000000151720049)↗
NASA NTRS1983

Planetary atmospheres

During the period from 1979 to 1982, Pioneer and Voyager missions to Venus, Jupiter, and Saturn combined with ongoing analysis and data collection from the Viking-Mars mission have contributed greatly to an understanding of planetary atmospheres. Research concerning the atmosphere of Venus is discussed, taking into account the atmospheric composition, the isotopic ratios of noble gases, data concerning temperature and pressure, information regarding various atmospheric layers, atmospheric density data, and models regarding the dynamical properties of Venus' atmosphere. Information concerning the atmosphere of Mars is also examined. Isotopic ratios of carbon and oxygen are considered along with altitude profiles, diurnal heating and cooling, and thermal tides and their relationship to dust storms. Attention is given to the temperature structure of the Jovian atmosphere, the morphology of Jovian cloud features, features in the Saturnian atmosphere, information with respect to the atmosphere of Titan, the rotation rates of Uranus, masses and densities of Uranus and Neptune, and the identification of Pluto as a small icy body.

Beebe, R.↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
NASA NTRS1993

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗
DOE OSTI2022

X-ray diffraction study of phase transformation dynamics of Fe and Fe-Si alloys along the shock Hugoniot using an x-ray free electron laser

The x-ray free electron laser (XFEL) enables probing a highly compressed material response at the subnanosecond timescale. Here we exploit the ultrafast XFEL pulse to combine reflection x-ray diffraction and laser-driven shock compression to perform a study of the phase transformation and stability in Fe and Fe-Si alloys. Our approach enables us to observe that solid-solid phase transformations occur in Fe and Fe- Si 8.5 wt % in ≤ 130 ps at ~ 130 GPa; no transformation is observed in Fe- Si 16 wt % up to 110 GPa. Density functional theory calculations predict similar phase relations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Faraday Rotation Measurements in High-Energy-Density Plasmas Using Shaped Laser Beams

Magnetic fields play an important role in plasma dynamics, yet it is a quantity difficult to measure accurately with physical probes, whose presence disturbs the very field they measure. The Faraday rotation of a polarized beam of light provides a mechanism to measure the magnetic field without disturbing the dynamics, and has been used with great success in astrophysics and high energy density plasma science, where physical probes cannot be used. Furthermore, the rotation is typically small, which degrades the accuracy of the measurement. Since polarization cannot be measured directly, detectors rely on a polarizer to measure a small change in beam intensity instead. In this work, we show how beam shaping can improve Faraday rotation measurements using an optical derivative setup. Since the rotation measurement is now strictly proportional to the beam shape and intensity, the system allows to improve the measurement accuracy simply by increasing the laser beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1980

Simultaneous measurements of magnetotail dynamics by IMP spacecraft

Changes in tail energy density during substorms in the magnetotail are given. In addition to plasma sheet thinnings seen prior to substorm onsets, a gradual decrease in plasma beta was detected in the deep tail which precedes onset and the more prominent plasma disappearance that typically accompanies it. The frequency of thinnings and the regions over which they occurred indicate that drastic changes in plasma sheet thickness are common features of substorms which occur at all locations across the tail.

Fairfield, D. H.↗
NASA NTRS1981

Chesapeake Bay plume dynamics from LANDSAT

LANDSAT images with enhancement and density slicing show that the Chesapeake Bay plume usually frequents the Virginia coast south of the Bay mouth. Southwestern (compared to northern) winds spread the plume easterly over a large area. Ebb tide images (compared to flood tide images) show a more dispersed plume. Flooding waters produce high turbidity levels over the shallow northern portion of the Bay mouth.

Munday, J. C., Jr.↗
NASA NTRS1993

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗