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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

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

Trajectory of a flying plasma mirror traversing a target with density gradient

It has been proposed that laser-induced relativistic plasma mirror can accelerate if the plasma has a properly tailored density profile. Such accelerating plasma mirrors can serve as analog black holes to investigate Hawking evaporation and the associated information loss paradox. Here we reexamine the underlying dynamics of mirror motion in a graded-density plasma to provide an explicit trajectory as a function of the plasma density and its gradient. Specifically, a decreasing plasma density profile (down-ramp) along the direction of laser propagation would in general accelerate the mirror. In particular, a constant-plus-exponential density profile would generate the Davies–Fulling trajectory with a well-defined analog Hawking temperature, which is sensitive to the plasma density gradient but not to the density itself. Furthermore, we show that without invoking nano-fabricated thin-films, a much lower density gas target at, for example, ~1 × 10 17 cm -3 , would be able to induce an analog Hawking temperature, k B T H ~3.1 × 10 -2 eV, in the far-infrared region. We hope that this would help to better realize the experiment proposed by Chen and Mourou.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

How accurately can we detect the splashback radius of dark matter haloes and its correlation with accretion rate?

ABSTRACT The splashback radius (Rsp) of dark matter haloes has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure Rsp and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover Rsp by comparing it with the value measured directly from particle dynamics. We show that estimates of Rsp from density and subhalo profiles correspond to different percentiles of the underlying Rsp distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the Rsp–Γdyn relation directly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

Large Space Optics: From Hubble to JWST and Beyond

If necessity truly is the mother of invention, then advances in lightweight space mirror technology have been driven by launch vehicle mass and volume constraints. In the late 1970 s, at the start of Hubble development, the state of the art in ground based telescopes was 3 to 4 meter monolithic primary mirrors with masses of 6000 to 10,000 kg - clearly too massive for the planned space shuttle 25,000 kg capability to LEO. Necessity led Hubble to a different solution. Launch vehicle mass constraints (and cost) resulted in the development of a 2.4 meter lightweight eggcrate mirror. At 810 kg (180 kg/m2), this mirror was approximately 7.4% of HST s total 11,110 kg mass. And, the total observatory structure at 4.3 m x 13.2 m fit snuggly inside the space shuttle 4.6 m x 18.3 m payload bay. In the early 1990 s, at the start of JWST development, the state of the art in ground based telescopes was 8 meter class monolithic primary mirrors (16,000 to 23,000 kg) and 10 meter segmented mirrors (14,400 kg). Unfortunately, launch vehicles were still constrained to 4.5 meter payloads and 25,000 kg to LEO or 6,600 kg to L2. Furthermore, science now demanded a space telescope with 6 to 8 meter aperture operating at L2. Mirror technology was identified as a critical capability necessary to enable the next generation of large aperture space telescopes. Specific telescope architectures were explored via three independent design concept studies conducted during the summer of 1996 (1). These studies identified two significant architectural constraints: segmentation and areal density. Because the launch vehicle fairing payload dynamic envelop diameter is approximately 4.5 meters, the only way to launch an 8 meter class mirror is to segment it, fold it and deploy it on orbit - resulting in actuation and control requirements. And, because of launch vehicle mass limits, the primary mirror allocation was only 1000 kg - resulting in a maximum areal density of 20 kg/m2. At the inception of JWST in 1996, such a capability did not exist. A highly successful technology development program was initiated resulting in matured and demonstrated mirror technology for JWST (2, 3). Today, the JWST 6.5 meter primary mirror has an areal density of 25 kg/m2 for a total mass of 625 kg or 9.6% of the total JWST observatory mass of 6,500 kg. Looking into the future, science requires increasing larger collecting apertures. Ground based telescopes are already moving towards 30+ meter mirrors. The only way to meet this challenge for space telescopes is via even lower areal density mirrors or on-orbit assembly or larger launch vehicles (4). The planned NASA Ares V with its 10 meter fairing and 55,000 kg payload to L2 eliminates this constraint (5).

Stahl, H. Philip↗
DOE OSTI2020

High susceptibility to adiabatic shear banding and high dynamic strength in tungsten heavy alloys with a high-entropy alloy matrix

The susceptibility to adiabatic shear banding (ASB) of tungsten heavy alloys (WHAs) is a key factor to determine their penetration performance for kinetic energy penetrators (KEPs). This work reports a novel WHA with strong susceptibility to ASB that uses a precipitation-hardening high-entropy alloy (HEA) as the binding phase of W particles. The results find that nanoscale L12 precipitates and high dislocation density in the matrix contribute a high dynamic yield strength of ~2300 MPa. Also, a narrow ASB with a width of ~12 μm forms in the novel WHA, which is only 6% of the width of ASB in traditional WHAs like W-NiFe, exhibiting a great susceptibility to ASB. The redissolution of precipitates into the matrix because of temperature rise during shearing could sharply soften the shearing location, which is suggested to be a very important factor of destabilizing mechanism to promote the development of ASB.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗
DOE OSTI2021

Preventing Electrolyte Decomposition on a Ca Metal Electrode Interface Using an Artificial Solid‐Electrolyte Interphase

Abstract Calcium ion batteries are gaining attention as alternatives to lithium‐ion technology because they offer comparable properties at reduced cost and improved safety. However, progress has been limited because of the inability to efficiently and reversibly plate and strip Ca metal anodes in organic electrolytes. Moreover, the inorganic components of the solid‐electrolyte interphase (SEI) that form via decomposition of the electrolyte often do not allow for the diffusion of Ca ions. In this work, an approach combining density functional theory and ab initio molecular dynamics (AIMD) simulations is utilized to show that the use of a preformed artificial SEI layer of amorphous can potentially prevent electrolyte decomposition. First, Ca is shown to be able to intercalate into an amorphous layer (up to Ca 1.5 Al 2 O 3 ) and diffuse through on a reasonable time scale. Through calculation of the density of states, the system is found to remain insulating up to the equilibrium stoichiometry. Finally, AIMD simulations with a realistic organic electrolyte environment are used to show that this calcinated layer completely prevents the decomposition of solvent molecules. This approach can provide a route to efficient rechargeable Ca ion batteries, paving the way for cheap large‐scale energy storage.

Young, Joshua↗
DOE OSTI2021

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗
NASA NTRS1995

High-Density Amorphous Ice, the Frost on Interstellar Grains

Most water ice in the universe is in a form which does not occur naturally on Earth and of which only minimal amounts have been made in the laboratory. We have encountered this 'high-density amorphous ice' in electron diffraction experiments of low-temperature (T less than 30 K) vapor-deposited water and have subsequently modeled its structure using molecular dynamics simulations. The characteristic feature of high-density amorphous ice is the presence of 'interstitial' oxygen pair distances between 3 and 4 A. However, we find that the structure is best described as a collapsed lattice of the more familiar low-density amorphous form. These distortions are frozen in at temperatures below 38 K because, we propose, it requires the breaking of one hydrogen bond, on average, per molecule to relieve the strain and to restructure the lattice to that of low-density amorphous ice. Several features of astrophysical ice analogs studied in laboratory experiments are readily explained by the structural transition from high-density amorphous ice into low-density amorphous ice. Changes in the shape of the 3.07 gm water band, trapping efficiency of CO, CO loss, changes in the CO band structure, and the recombination of radicals induced by low-temperature UV photolysis all covary with structural changes that occur in the ice during this amorphous to amorphous transition. While the 3.07 micrometers ice band in various astronomical environments can be modeled with spectra of simple mixtures of amorphous and crystalline forms, the contribution of the high-density amorphous form nearly always dominates.

Jenniskens, P.↗
DOE OSTI2021

Impact of functional groups on the electrocatalytic hydrogenation of aromatic carbonyls to alcohols

We report that electrocatalytic hydrogenation (ECH) of biomass-derived feedstocks has a critical dependence on the molecular structure of the organic and its adsorption on the electrode surface. In this study, we investigated the role of functional groups in the adsorption of the organic molecule on the charged Pd (111) surface and its subsequent effect on organic reduction in electrochemical hydrogenation of organic molecules. With three aromatic carbonyls of benzaldehyde (BZD), acetophenone (ACE), and vanillin (VAN), we rationalize molecular-scale adsorption and interfacial charge transfer processes by employing density-functional-theory based ab initio molecular dynamics simulations. We observe that the functional group and electrode charge strongly affect the proximity of organic molecule to the Pd (111) surface, where distances of aromatic ring and carbonyl group of the organic on the electrode change distinctively with functional groups and charge state of electrode, which strongly impact reduction of organics on the surface. Calculations of differential electron density show the strongest reduction with benzaldehyde via interfacial electron transfer from the charged Pd surface. We also observe that the interaction between the functional groups and solvent (VAN > BZD > ACE) significantly influence the organic interaction with the charged electrode (BZD > VAN > ACE), resulting in the net interaction energy between the organic and the electrode in the order of BZD > ACE > VAN. Experimental measurement of ECH rate also show the same trend of the net interaction energy. These results demonstrate the significance of solvent effect on the reducibility of organic molecules on electrodes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials (Final Technical Report)

The goal of this project was to conduct dynamic compression experiments and predictive simulations to reveal the fundamental high-energy-density (HED) physics of amorphous carbon. These results are essential for assessing amorphous carbon as a potential ablator material for next generation of inertial confinement fusion (ICF) capsules. We made significant progress in exploring the HED properties of amorphous carbon through experiments at Omega EP Laser and the European XFEL, in addition to billion-atom, quantum-accurate molecular dynamics (MD) simulations. Through our joint experimental and simulation program, we mapped the phase diagram of amorphous carbon, uncovering its range of metastability and identifying phase transitions to diamond and liquid carbon along the Hugoniot and at higher pressures using double shock compression pathways. Our findings indicate that amorphous carbon melts at significantly lower shock pressures than high-density carbon (diamond). However, nanocrystalline diamond nucleates across a broad range of pressures and temperatures. This emergence of the nanocrystalline microstructure during compression can negatively impact the planarity of the shock front and potentially trigger ablator/fuel mixing during Inertial Fusion Energy (IFE) applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Enabling Industrial Re-Use of Large-Format Additive Manufacturing Molding and Tooling

Large-format additive manufacturing (LFAM) is an enabling manufacturing technology capable of producing large parts with highly complex geometries for a wide variety of applications, including automotive, infrastructure/construction, and aerospace mold and tooling. In the past decade, the LFAM industry has seen widespread use of bio-based, glass, and/or carbon fiber reinforced thermoplastic composites which, when printed, serve as a lower-cost alternative to metallic parts. One of the highest-volume materials utilized by the industry is carbon fiber (CF)-filled polycarbonate (PC), which in out-of-autoclave applications can achieve comparable mechanical performance to metal at a significantly lower cost. Previous work has shown that if this material is recovered at various points throughout the manufacturing process for both the lab and pilot scale, it can be mechanically recycled with minimal impacts on the functional performance and printability of the material while significantly reducing the feedstock costs. End-of-life (EOL) CF-PC components were processed through industrial shredding, melt compounding, and LFAM equipment, followed by evaluation of the second-life material properties. Experimental assessments included quantitative analysis of fiber length attrition, polymer molecular weight degradation using gel permeation chromatography (GPC), density changes via pycnometry, thermal performance using dynamic mechanical analysis (DMA), and mechanical performance (tensile properties) in both the X- and Z-directions. Results demonstrated a 24.6% reduction in average fiber length compared to virgin prints, accompanied by a 21% decrease in X-direction tensile strength and a 39% reduction in tensile modulus. Despite these reductions, Z-direction tensile modulus improved by 4%, density increased by 6.8%, and heat deflection temperature (HDT) under high stress retained over 97% of its original value. These findings underscore the potential for integrating mechanically recycled CF-PC into industrial LFAM applications while highlighting the need for technological innovations to mitigate fiber degradation and enhance material performance for broader adoption. This critical step toward circular material practices in LFAM offers a pathway to reducing feedstock costs and environmental impact while maintaining functional performance in industrial applications.

additive manufacturing↗
DOE OSTI2021

The LISE package: solvers for static and time-dependent superfluid local density approximation equations in three dimensions

Nuclear implementation of the density functional theory (DFT) is at present the only microscopic framework applicable to the whole nuclear landscape. The extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach, the superfluid local density approximation (SLDA) and its extension to time-dependent situations, time-dependent superfluid local density ap- proximation (TDSLDA), have been extensively used to describe various static and dynamical problems in nuclear physics, neutron star crust, and cold atom systems. In this paper, we present the codes that solve the static and time-dependent SLDA equations in three-dimensional coordinate space without any symmetry restriction. These codes are fully parallelized with the message passing interface (MPI) library and take advantage of graphic processing units (GPU) for accelerating execution. The dynamic codes have checkpoint/restart capabilities and for initial conditions one can use any generalized Slater determinant type of wave function. The code can describe a large number of physical problems: nuclear fission, collisions of heavy ions, the interaction of quantized vor- tices with nuclei in the nuclear star crust, excitation of superfluid fermion systems by time dependent external fields, quantum shock waves, domain wall generation and propagation, the dynamics of the Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex rings dynamics, generation and dynamics of quantized vortices, their crossing and recombinations and the incipient phases of quantum turbulence.

Jin, Shi↗
DOE OSTI2024

Self-Consistent Theory for Structural Relaxation, Dynamic Bond Exchange Times, and the Glass Transition in Polymeric Vitrimers

Here, we formulate a statistical mechanical theory for how dynamic bond exchange influences the activated hopping-driven relaxation of Kuhn segments in dynamically cross-linked networks or vitrimers over a wide range of temperatures and cross-link densities. The key new methodological aspect is to address in a self-consistent manner the dynamic consequences of bond exchange on the Kuhn segmental alpha relaxation, and vice versa. The predicted temperature dependence of the segmental alpha time of vitrimers at high temperatures remains the same as that of permanent networks, but at lower temperatures, a significant acceleration of relaxation occurs due to bond exchanges. From a mechanistic perspective, the vitrimer local cage barrier is very weakly affected by bond exchange, while the collective elastic barrier contribution decreases significantly in the deeply supercooled regime. The vitrimer glass transition temperature is predicted to grow linearly with the square root of the cross-link density, as previously found for permanent networks. Material-specific chemical effects such as cross-linker size relative to that of the normal Kuhn segment or special attraction of a cross-linker with polymers are crudely considered based on model calculations. The theory is quantitatively applied to recent experiments on dry ethylene vitrimers. Good agreements are found including that the bond exchange time follows an Arrhenius law at high enough temperatures but upward non-Arrhenius deviations emerge in the deeply supercooled regime, and a collapsed master curve of the segment alpha time exists over a wide range of cross-link densities and temperatures. Possible extensions to treat dynamic heterogeneity effects and penetrant transport in vitrimers are briefly discussed.

chemical calculations↗
DOE OSTI2021

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples—random charge densities and periodic charge densities—with different scaling properties. Additionally, under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗