Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modulus”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pressure-driven symmetry transitions in dense H 2 O ice

X-ray diffraction and Raman spectroscopy of H 2 O (ice) structures are measured under static compression in combination with grain normalizing heat treatment via direct laser heating. Here, we report the transition from cubic ice-VII to a structure of tetragonal symmetry, ice-VII t at 5.1 ± 0.5 GPa. This is succeeded by the H-bond symmetrization transition occurring at a pressure of 30.9 ± 3 GPa. Both experimental observations are supported by simulated Raman spectra from density-functional theory quantum calculations. The transition to H-bond symmetrization is evidenced by the reversible emergence of its characteristic Raman mode and a 2.5-fold increase in bulk modulus, implying a significant increase in bonding strength.

58 GEOSCIENCES↗
DOE OSTI2023

Shock-induced melting of [100] lithium fluoride: Sound speed and Hugoniot measurements to 230 GPa

Although [100] lithium flouride (LiF) is the most widely used optical window material in dynamic compression experiments, its high stress (>100 GPa) shock compression response, including melting, is not well understood. To address this need, we measured wave profiles in plate impact experiments to determine the Hugoniot states and longitudinal sound speeds in [100] LiF crystals shock compressed to 231 GPa. The measured peak states are fitted well by a linear shock velocity-particle velocity relation, providing an accurate determination of the LiF Hugoniot curve to significantly higher stresses than previous experiments. The longitudinal sound speeds show a near linear increase with density compression to 182 GPa. Between 182 GPa and 195 GPa, the sound speed and the longitudinal modulus decrease abruptly, due to shock-induced melting. The increasing sound speeds and moduli at higher stresses suggest that shock compressed LiF is fully liquid at 195 GPa and above, allowing determination of the Gruneisen parameter for liquid LiF. Here, the melt stress determined here differs from that predicted by current multiphase equations of state for LiF. Our results provide important insight into the high stress solid and liquid states of shock compressed LiF and point to the need for an improved multiphase equation of state at high pressures and high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Real-space atomic dynamics in metallic liquids investigated by inelastic neutron scattering

Understanding the dynamics of liquids at the atomic level remains a major challenge. Even though viscosity is one of the most fundamental properties of liquids, its atomistic origin is not fully elucidated. Through inelastic neutron scattering experiment on levitated metallic liquid droplets, the time-dependent pair correlation function, the Van Hove function, was determined for Z⁢r 50 ⁢C⁢u 50 and Z⁢r 80 ⁢P⁢t 20 liquids at various temperatures. The time for change in local atomic connectivity, τ LC , which is the timescale of atomic bond cutting and forming, is estimated based on the exponential decay of the nearest neighbor peak of the Van Hove function. At high temperatures above the crossover temperature T A , τ LC is equal to the Maxwell relaxation time, τ M =η/G ∞ , where η is the macroscopic shear viscosity and G ∞ is the high-frequency shear modulus. Below T A the ratio of τ M /τ LC increases with decreasing temperature, indicating increased atomic cooperativity as predicted by molecular dynamics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elastic properties of the W 0.75 R e 0.25 alloy at high pressure up to 183 GPa

The high pressure equation of state for the W 0.75 R e 0.25 alloy is experimentally determined up to 183 GPa with synchrotron angle-dispersive powder x-ray diffraction in the diamond-anvil cell and to ∼925 GPa with density-functional theory. W-Re alloys are used in many industrial high-temperature applications and as a confining gasket material in high-pressure diamond-anvil cell research. The inclusion of 25 wt. % Re achieves the highest performance in terms of strength and ductility while also maintaining the body-centered-cubic (bcc) crystal structure, yet to date there has been no investigation into its elastic behavior at high pressure. We present the experimentally and theoretically determined volumetric and elastic pressure response and systematically compare these results to other W-Re alloys, finding that the bulk modulus of W-Re alloys varies nonlinearly with Re content and W 0.75 Re 0.25 becomes less incompressible than W at 85 GPa. Published by the American Physical Society 2025

Alloys↗
DOE OSTI2021

Asymmetric reheating by primordial black holes

Here we investigate Hawking evaporation of a population of primordial black holes (PBHs) prior to big bang nucleosynthesis as a mechanism to achieve asymmetric reheating of two sectors coupled solely by gravity. While the visible sector is reheated by the inflaton or a modulus, the dark sector is reheated by PBHs. Compared to inflationary or modular reheating of both sectors, there are two advantages: (i) inflaton or moduli mediated operators that can subsequently thermalize the dark sector with the visible sector are not relevant to the asymmetric reheating process and (ii) the mass and abundance of the PBHs provide parametric control of the thermal history of the dark sector, and in particular the ratio of the temperatures of the two sectors. Asymmetric reheating with PBHs turns out to have a particularly rich dark sector phenomenology, which we explore using a single self-interacting real scalar field in the dark sector as a template. Four thermal histories, involving nonrelativistic and relativistic dark matter at chemical equilibrium, followed by the presence or absence of cannibalism, are explored. These histories are then constrained by the observed relic abundance in the current Universe and the Bullet Cluster. The case where PBHs dominate the energy density of the Universe, and reheat both the visible as well as the dark sectors, is also treated in detail.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Anisotropic neutron star crust, solar system mountains, and gravitational waves

"Mountains" or non-axisymmetric deformations of rotating neutron stars (NS) efficiently radiate gravitational waves (GW). We consider analogies between NS mountains and surface features of solar system bodies. Both NS and moons such as Europa or Enceladus have thin crusts over deep oceans while Mercury has a thin crust over a large metallic core. Thin sheets may wrinkle in universal ways. Europa has linear features, Enceladus has "Tiger" stripes, and Mercury has lobate scarps. NS may have analogous features. The innermost inner core of the Earth is anisotropic with a shear modulus that depends on direction. If NS crust material is also anisotropic this will produce an ellipticity, when the crust is stressed, that grows with spin frequency. This yields a braking index (log derivative of spin down rate assuming only GW spin down) very different from n=5 and could explain the maximum spin observed for neutron stars and a possible minimum ellipticity of millisecond pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Microscopic ergodicity breaking governs the emergence and evolution of elasticity in glass-forming nanoclay suspensions

We report a study combining x-ray photon correlation spectroscopy (XPCS) with in situ rheology to investigate the microscopic dynamics and mechanical properties of aqueous suspensions of the synthetic hectorite clay Laponite, which is composed of charged, nanometer-scale, disk-shaped particles. The suspensions, with particle concentrations ranging from 3.25 to 3.75 wt%, evolve over time from a fluid to a soft glass that displays aging behavior. The XPCS measurements characterize the localization of the particles during the formation and aging of the soft-glass state. The fraction of localized particles f 0 increases rapidly during the early formation stage and grows more slowly during subsequent aging, while the characteristic localization length r loc steadily de- creases. Despite the strongly varying rates of aging at different concentrations, both f 0 and r loc scale with the elastic shear modulus G' in a manner independent of concentration. During the later aging stage, the scaling between r loc and G' agrees quantitatively with a prediction of naive mode coupling theory. Breakdown of agreement with the theory during the early formation stage indi- cates the prevalence of dynamic heterogeneity, suggesting the soft solid forms through precursors of dynamically localized clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure dependence of intermediate-range order and elastic properties of glassy Baltic amber

Amber is a unique example of a fragile glass that has been extensively aged below its glass transition temperature, thus reaching a state that is not accessible under normal experimental conditions. In this paper, we studied the medium-range order of Baltic amber by x-ray diffraction (XRD) at high pressures. The pressure dependences of the low-angle XRD intensity between 0 and 5 Å -1 were measured from 0 to 7.3 GPa by the energy-dispersive XRD. The first diffraction peak at 1.1 Å -1 and ambient pressure has a doublet structure consisting of the first sharp diffraction peak (FSDP) at 1.05 Å -1 and the second feature at 1.40 Å -1 . The peak position and the width of the FSDP increase as the pressure increases, while the intensity of the FSDP decreases. Below P 0 = 2.4 GPa, the rapid increase of the FSDP peak position was observed, while above P 0 , the gradual increase was observed. Below P 0 , voids and holes in a relatively low-density state are suppressed, whereas above P 0 , the suppression becomes mild. Such a change suggests the crossover from the low- to high-density state at P 0 . There is a close correlation between the pressure dependence of XRD and previously reported sound velocity results. The correlation between the mean-square fluctuation of the shear modulus on the nanometer scale and fragility in amber and other glass formers is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Different bonding type along each crystallographic axis: Computational study of poly($\textit{p}$-phenylene terephthalamide)

Poly($\textit{p}$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $\textit{a}$ axis, hydrogen bonding along its $\textit{b}$ axis, and covalent bonding along its $\textit{c}$ axis. In this work, we explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $\textit{a}$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $\textit{b}$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $\textit{ab}$ rectangle, is predicted by SCAN.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tensor renormalization group approach to critical phenomena via symmetry-twisted partition functions

The locality of field theories strongly constrains the possible behaviors of symmetry-twisted partition functions, and thus they serve as order parameters to detect low-energy realizations of global symmetries, such as spontaneous symmetry breaking (SSB). We demonstrate that the tensor renormalization group (TRG) offers an efficient framework to compute the symmetry-twisted partition functions, which enables us to detect the symmetry-breaking transition and also to study associated critical phenomena. As concrete examples of SSB, we investigate the two-dimensional (2D) classical Ising model and the three-dimensional (3D) classical 𝑂⁡(2) nonlinear sigma model, and we identify their critical points solely from the twisted partition function. By employing the finite-size scaling argument, we find the critical temperature 𝑇 𝑐 = 2.2017⁢(2) with the critical exponent 𝜈 = 0.663⁢(33) for the 3D 𝑂⁡(2) model. In addition, we also study the Berezinskii–Kosterlitz–Thouless (BKT) criticality of the 2D classical 𝑂⁡(2) model by extracting the helicity modulus from the twisted partition functions, and we obtain the BKT transition temperature, 𝑇 BKT = 0.8928⁢(2).

lattice field theory↗
DOE OSTI2024

Multidimensional Rietveld refinement of high-pressure neutron diffraction data of PbNCN

High-pressure neutron powder diffraction data from PbNCN were collected on the high-pressure diffraction beamline SNAP located at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (Tennessee, USA). The diffraction data were analyzed using the novel method of multidimensional (two dimensions for now, potentially more in the future) Rietveld refinement and, for comparison, employing the conventional Rietveld method. To achieve two-dimensional analysis, a detailed description of the SNAP instrument characteristics was created, serving as an instrument parameter file, and then yielding both cell and spatial parameters as refined under pressure for the first time for solid-state cyanamides/carbodiimides. The bulk modulus B 0 = 25.1 (15) GPa and its derivative B ′ 0 = 11.1 (8) were extracted for PbNCN following the Vinet equation of state. Surprisingly, an internal transition was observed beyond 2.0 (2) GPa, resulting from switching the bond multiplicities (and bending direction) of the NCN 2− complex anion. The results were corroborated using electronic structure calculation from first principles, highlighting both local structural and chemical bonding details.

Meinerzhagen, Yannick↗
DOE OSTI2025

Machine learning inversion from scattering for mechanically driven polymers

A machine learning inversion method is developed for analyzing scattering functions of mechanically driven polymers and extracting the corresponding feature parameters, which include energy parameters and conformation variables. The polymer is modeled as a chain of fixed-length bonds constrained by bending energy, and it is subject to external forces such as stretching and shear. We generate a data set consisting of random combinations of energy parameters, including bending modulus, stretching and shear force, along with Monte Carlo-calculated scattering functions and conformation variables such as end-to-end distance, radius of gyration and off-diagonal component of the gyration tensor. The effects of the energy parameters on the polymer are captured by the scattering function, and principal component analysis ensures the feasibility of the machine learning inversion. Finally, we train a Gaussian process regressor using part of the data set as a training set and validate the trained regressor for inversion using the rest of the data. The regressor successfully extracts the feature parameters.

Gaussian process regressors↗
DOE OSTI2022

Negative linear compressibility in Se at ultra-high pressure above 120 GPa

A series of in situ synchrotron X-ray diffraction (XRD) measurements were carried out, combined with first-principles calculations, to study structural phase transitions of selenium at high pressures and room temperature. Several phase transitions were observed, among which an isostructural phase transition was found at around 120 GPa for the first time. Evolved from the rhombohedral (space group R$\bar{3}$m) structure (Se-V), the new phase (Se-V') exhibited an interesting increase of lattice parameter a at pressures from 120 to 148 GPa, known as negative linear compressibility (NLC). The discovery of NLC behavior observed in this work is mainly attributed to the accuracy and fine steps controlled by the membrane system for in situ XRD data collected with an exposure time of 0.5 s. After 140 GPa, a body-centered cubic (b.c.c.) structure Se-VI (space group Im$\bar{3}$m) was formed, which remains stable up to 210 GPa, the highest pressure achieved in this study. The bulk moduli of phases Se-V, Se-V' and Se-VI were estimated to be 83 ± 2, 321 ± 2 and 266 ± 7 GPa, respectively, according to the P–V curve fit by the third-order Birch–Murnaghan equation of state. The Se-V' phase shows a bulk modulus almost 4 times larger than that of the Se-V phase, which is mainly due to the effect of its NLC. NLC in a higher pressure range is always more significant in terms of fundamental mechanism and new materials discovery, yet it has barely been reported at pressures above 100 GPa. This will hopefully inspire future studies on potential NLC behaviors in other materials at ultra-high pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Application of Principal Component Analysis for the Monitoring of the Aging Process of Nuclear Electrical Cable Insulation

To ensure the safe operation of nuclear power plants (NPPs), it is critical to understand how NPP electrical cable insulation will degrade under different service environments. In this study, various nondestructive examination methods were selected to evaluate degradation of electrical cable insulation after aging. Elongation at break, indenter modulus, relaxation constant, mass change, total color difference, and carbonyl index were collected stepwise on cross-linked polyethylene (XLPE) cable insulations after predetermined exposure intervals. Three different insulation aging scenarios were investigated: 1) simultaneous irradiation and heating at 150°C, 2) heating at 150°C followed by corresponding times of irradiation at ambient temperature, and 3) irradiation at ambient temperature followed by heating at 150°C for the same duration. A constant dose rate of 300 Gy/hr was selected with a total gamma irradiation absorbed dose up to 320 kGy. To gain insight into the long-term performance of NPP electrical cable insulation, principal component analysis (PCA), a data-driven approach, was utilized to identify key indicators of cable insulation degradation. By reducing the dimensionality of the data while retaining degradation information, PCA was used to highlight the changes in the measured properties under gamma irradiation according to total absorbed dose and the different aging scenarios.

Li, Donghui↗
DOE OSTI2020

Stability of Tin-Lead Halide Perovskite Solar Cells

Tin-lead perovskites have low band gaps (1.2 eV - 1.4 eV) that offer complementary absorption spectra to mixed halide perovskites with band gaps in the range of 1.6 eV - 1.85 eV. In combination, tin-lead and mixed halide perovskites can be used to make efficient all-perovskite tandem solar cells. While these can achieve high efficiency, they are hampered by the unproven longterm stability of tin-containing perovskites. We make the first demonstration of tin-lead perovskite solar cells that pass benchmark 1000-hour tests of stability under stressors of heat, light, and atmospheric exposure. We identify that mixed tin-lead perovskites oxidise by a different chemical pathway, one that is much less favourable, than the oxidation pathway followed by pure tin perovskites, making them significantly more stable. Fortuitously, mixed tin-lead perovskites also have band gaps that are close to ideal for the rear cell of an all-perovskite monolithic two-junction tandem. Beyond selecting the composition of the perovskite to be one that has inherently suppressed tendency to oxidise, we also fabricate a novel heterojunction-based architecture that eliminates the commonly used acidic PEDOT:PSS layer, which leads to enhanced stability. By packaging the resulting solar cell in a glass-on-glass package with a pliable low-elastic-modulus polyolefin encapsulant and a butyl rubber edge seal, we also demonstrate that the low band gap perovskite solar cell that passes the IEC damp heat test for thin film solar cells - a milestone that demonstrates the removal of a significant impediment to the commercialization of high efficiency perovskite tandem solar cells.

14 SOLAR ENERGY↗