SEARCH · Search NASA
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.
SLIP (Surrogate Launching and Integration Platform) is a software ecosystem for downloading and running ML/AI benchmarks. SLIP automatically downloads and sets up the code and data running a benchmark. The large data sets and model code are cached locally with a specific ID within a designated cache directory.
Anyons are two-dimensional particles with fractional exchange statistics that emerge as elementary excitations of fractional quantum Hall phases. Experimentally, their exchange statistics can be measured in the edge-state Fabry-Pérot interferometer, wherein the presence of 𝑁 𝑞𝑝 localized anyons contributes a phase 𝑁 𝑞𝑝 𝜃 𝑎 to the interference pattern where 𝜃 𝑎 is twice the exchange phase. Here we report the observation of large, hysteretic phase jumps in a monolayer graphene Fabry-Pérot interferometer at 𝜈 = 1/3. When the filling factor is increased from 𝜈 < 1/3 toward the center of the plateau, we observe phase slips with magnitude Δ𝜃 ≈ 2𝜋/3, consistent with the addition of individual quasiparticles to the interferometer bulk. These phase slips occur as instantaneous jumps in the interference signal, with intervals between the jumps indicating quasiparticle equilibration times exceeding 20 min. We use this long timescale to investigate the effect of changes in interferometer area 𝐴 𝐼 and 𝑁 𝑞𝑝 independently at fixed magnetic field, revealing a striking memory effect in the phase slip magnitude. In particular, as the 𝜈 =1/3 plateau is approached from higher filling, we observed phase slips with Δ𝜃 significantly larger than 2𝜋/3 over the same range of gate voltage where quantized jumps are seen for increasing 𝜈. We discuss this asymmetry in terms of bulk-edge coupling of quasiparticles localized near the edge or in the bulk, and argue that this effect can be qualitatively reconciled with theoretical expectations for strongly interacting quasiparticles in the presence of weak disorder and strongly nonequilibrium charge dynamics. Besides providing a replication of interferometric measurements sensitive to 𝜃 𝑎 , our results highlight the key role played by charge dynamics on signatures of the anyon phase, and demonstrate that fractional quasiparticles can be indefinitely localized in nonequilibrium configurations.
The mechanical behavior of two nanocrystalline Al alloys, Al-Mg-Y and Al-Fe-Y, is investigated with in-situ micropillar compression testing. Both alloys were strengthened by a hierarchical microstructure including grain boundary segregation, nanometer-thick amorphous complexions, carbide nanorod precipitates with sizes of a few nanometers, and submicron-scale intermetallic particles. The maximum yield strength of the Al-Mg-Y system is measured to be 950 MPa, exceeding that of the Al-Fe-Y system (680 MPa), primarily due to a combination of more carbide nanorods and more amorphous complexions. Both alloys exhibited yield strengths much higher than those of commercial Al alloys, and therefore have great potential for structural applications. However, some micropillar specimens were observed to plastically soften through shear banding. Post-mortem investigation revealed that intermetallic-free deformation pathways of a few micrometers in length were responsible for this failure. Further characterization showed significant grain growth within the shear band. The coarsened grains maintained the same orientation with each other, pointing to grain boundary mechanisms for plastic flow, specifically grain rotation and/or grain boundary migration. Here, the presence of intermetallic particles makes it difficult for both matrix and intermetallic grains to rotate into the same orientation due to the different lattice parameters and slip systems. Therefore, we are able to conclude that a uniform distribution of intermetallic particles with an average spacing less than the percolation length of shear localization can effectively prevent the maturation of shear bands, offering a design strategy for high-strength nanocrystalline Al alloys with both high strength and stable plastic flow.
Cross-slip is a thermally activated process by which screw dislocation changes its glide plane to another slip plane sharing the same Burgers vector. The rate at which this process happens is determined by a Boltzmann type expression that is a function of the screw segment length and the stress acting on the dislocation. In continuum dislocation dynamics (CDD), the information regarding the length of the screw dislocation segment and local stress state on dislocations are lost due to the coarse-grained representation of the density. Here, in this work, a data driven approach to characterize the lost information by analyzing the discrete dislocation configurations is proposed to enable cross-slip modeling in the CDD framework in terms of the coarse-grained dislocation density and stress fields. The analysis showed that the screw segment length follows an exponential distribution, and the stress fluctuations, defined as the difference between the stress on the dislocations and the mean field stress in CDD, follows a Lorentzian distribution. A novel approach for cross slip implementation in CDD employing the screw segment length and stress fluctuation statistics was proposed and rigorously tested by comparing the CDD cross-slip rates with discrete dislocation dynamics (DDD) rates. This approach has been applied in conjunction with three cross-slip models used in DDD simulations differing mainly in the functional form of cross slip activation energy. It was found that different cross-slip activation energy formulations yielded different cross-slip rates, yet the effect on mechanical stress-strain response and dislocation density evolution was minimal for the [001] type loading.
Lattice oxygen redox (LOR) in P2-type layered oxides is an effective strategy to break through the limit of energy density of conventional cathodes due to its high redox potential (>4 V vs. Na + /Na) as well as extra capacity. Nevertheless, LOR induced local structure distortion and irreversible phase transitions cause serious electrochemical performance degradation, hindering the practical applications. In this work, we propose that the generation of the OP4 phase can be replaced with the Z phase by introducing Sb element with higher ionic potential and strong covalent bonds within the TMO 6 octahedron. Z phase transition is realized by constraining interlayer slipping between adjacent TM layers compared to OP4, which reduces the strain in the layered structure, lowers the Na+ diffusion energy barrier and creates more efficient Na+ diffusion channels. Consequently, Sb-substituted oxides demonstrates excellent kinetics, rate capability (79 mAh g -1 at 1 A g -1 ) in half cell and a high energy density of 487 Wh kg -1 (on cathode) in full cell.
Earthquakes can grow either monotonically from a single, stressed patch or through linking multiple stressed regions. The distinction has implications for magnitude predictability with single ruptures requiring knowledge of the local stress state, while linked ruptures require knowing the global stress and energy distribution. Here, we use a laboratory fault that allows direct observation of slip to determine which fault conditions promote linked ruptures, and how to observationally evaluate their likelihood. Higher concentrations of normal stress due to increased normal force, applied stress asperities, or larger heterogeneity between the samples, all lead to significant increases in linked rupture likelihood. The mean radiated energy enhancement factor (REEF) of large events is an excellent proxy for linked event likelihood, and a combination of REEF and duration can identify the rupture style of some but not all events. The results imply that globally observed variations of REEF can be interpreted as variations in stress concentration.
We examine the stress field in Southern California with a focus on the region near the South Central Transverse Ranges (SCTR), using a refined stress inversion methodology. We independently invert declustered and aftershocks focal mechanisms from the years 1981 to 2017. Comparisons between the stress field orientations and the stress ratios from these inversions provide information on local sources of the stress field. Regionally, the SHmax tends to be toward NNE and the stress ratios vary from transtensional stress regime near the Eastern California Shear Zone (ECSZ), to strike-slip faulting near the SCTR, and toward transpression near the Western Transverse Ranges. Detailed analysis of stress parameters near the SCTR indicates deviations from the regional strike-slip faulting including the following: (1) The San Bernardino Mountain area shows S Hmax direction toward NNW and transpressional stress components likely associated with the motions of the San Andreas Fault and ECSZ; (2) the Cajon Pass and San Gorgonio Pass show transpressional stress regime near the bottom of the seismogenic zone likely associated with the elevated topography; and (3) in Crafton Hills, rotation of the principal stress plunges and SHmax direction and transtensional stress regime below ~10 km, along with lower estimated apparent friction coefficient, suggest a weak fault possibly associated with deep creep. The performed multiscale analysis resolves effects of regional and local loadings. The average stress parameters in the analyzed 37 years do not show significant temporal variations of stress field near the SCTR region.
In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.
Additively manufactured (AM) 316 L stainless steel (SS) after hot isostatic pressing (HIP) was found to exhibit superior resistance to irradiation-assisted stress corrosion cracking (IASCC) in high-temperature water, as compared to wrought 316 L SS. The well-accepted IASCC factors of radiation-induced segregation (RIS) and radiation hardening are not accurate descriptions of IASCC susceptibility in this case. A decreased strain localization along grain boundaries, caused by dislocation channel broadening, was confirmed to suppress crack initiation. A unique distribution of irradiation-induced defects in HIP AM 316 L SS eased dislocation cross-slip compared to those in the wrought counterpart, thus increasing the channel width near the grain boundaries. Here for the first time, this study highlights the importance of dislocation channel broadening as a potential mechanism to further improve the IASCC resistance of 316 L SS and provides direct experimental evidence based on commercial-grade materials.
Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.
Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.
The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.
A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.
In order to identify the mechanism responsible for the formation of charge-density waves (CDW) in cuprate superconductors, it is important to understand which aspects of the CDW’s microscopic structure are generic and which are material-dependent. Here, we show that, at the local scale probed by NMR, long-range CDW order in YBa 2 Cu 3 O y is unidirectional with a commensurate period of three unit cells ( λ = 3 b ), implying that the incommensurability found in X-ray scattering is ensured by phase slips (discommensurations). Furthermore, NMR spectra reveal a predominant oxygen character of the CDW with an out-of-phase relationship between certain lattice sites but no specific signature of a secondary CDW with λ = 6 b associated with a putative pair-density wave. These results shed light on universal aspects of the cuprate CDW. In particular, its spatial profile appears to generically result from the interplay between an incommensurate tendency at long length scales, possibly related to properties of the Fermi surface, and local commensuration effects, due to electron-electron interactions or lock-in to the lattice.
This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.
This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.
Diamond is, by virtue of its high bonding forces and Peierls-Nabarro barrier stresses, among the hardest materials on earth due to the difficulty of generating and moving dislocations. We demonstrate, using molecular dynamics and analytical calculations, that the generation of defects is dependent on loading orientation. Shock-loading single-crystal diamond along [001] and [011] to a stress of 137 GPa did not reveal dislocations, whereas loading along [111] generated profuse dislocations. The introduction of a void at the latter orientation generated dislocations at 72.6 GPa. Two slip systems were identified: <011>{100} and <112>{111}. The threshold for plastic deformation in diamond is orientation dependent and significantly lower than reported in previous studies: 0.14 G (shear modulus). The void collapse generated localized amorphization for [001]. In conclusion, these results on the anisotropy of plasticity and void effects have relevance for the improvement of symmetry in the collapse of diamond capsules in inertial fusion experiments.
It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.