Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATRIX STRESS CALCULATION”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

LaFeSi–LaFe 13-x Si x composites: Modulating magnetic and magnetocaloric properties through inherent stress manipulation

We examine structural and magnetic properties of a series of La–Fe–Si alloys in the region of concentrations where they naturally form two-phase LaFeSi–LaFe 13-x Si x composites with variable content and connectivity of LaFe 13-x Si x grains distributed within the LaFeSi matrix. Theoretical calculations confirm that the LaFeSi constituent is magnetically and structurally inert below room temperature and at pressures between -10 and 10 GPa. The LaFe 13-x Si x constituent, on the other hand, is magnetically and structurally active: it exhibits first-order magnetostructural transformations that, in addition to x Si , can be controlled with temperature, magnetic field, and pressure. In composites where the concentration of the inactive constituent is ~70 wt. % or greater, the standard, single-step, LaFe 13-x Si x first-order phase transformation proceeds in two steps separated by over 30 K in a zero magnetic field. Increasing the magnetic field recouples the two steps and restores the single-step phase transformation pathway. We analyze the roles of stresses caused by both thermal expansion mismatch and the first-order magnetic phase transition in LaFe 13-x Si x to rationalize the observed physical behaviors that emerge as the temperature or/and magnetic field vary.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local residual stresses and microstructure within recrystallizing grains in iron

The intragranular strains and the microstructure within recrystallizing grains in single phase pure iron have been systematically studied in 3D using synchrotron micro-diffraction. In contradiction to common knowledge we observe that residual elastic strains of up to 1–2 × 10 –3 are present in the recrystallizing grains. The local strain variations within individual grains can be similar in magnitude to those between different grains. The effects of processing parameters, including rolling reduction, recrystallized volume fraction, and cooling rate, as well as microstructural parameters, including grain size, orientation and internal misorientation, on the development of local residual elastic strains and microstructure have been quantified. It is suggested that the difference in local defect density between recrystallized grains and the surrounding deformed matrix, as well as the redistribution of medium to long-range residual stresses in the deformed grains, are the main reasons for the development of such strains in the recrystallized grains. The local residual stress levels are calculated and their possible impact on the mechanical properties and recrystallization processes are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

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

Real-Time In Situ Neutron Diffraction Investigation of Phase-Specific Load Sharing in a Cold-Rolled TRIP Sheet Steel

Real-time in situ neutron diffraction was used to investigate the loading-direction dependence of phase-specific load sharing in a cold-rolled TRIP 780 sheet steel under tension at room temperature. A slightly lower yield stress and relatively faster initial strain hardening were observed for the tension along the rolling direction (RD), than those along the transverse direction (TD) and the diagonal direction (DD). In all three loading directions, the redistribution of inter-phase and inter-granular load sharing occurs upon the yields of the BCC matrix and the FCC phase. The FCC phase is found to yield earlier along the RD than those along the TD and DD. Based on the Hooke’s law and the rule of mixture, it is found that the lattice strains of BCC-(211)/FCC-(311) reflections via single peak fit are better representative of phase elastic strains than the strains calculated via the lattice parameters by Rietveld refinement. The quantitative evolution of phase-specific stress indicates a rapid strain hardening in the FCC phase after yielding during the RD tension. Possible mechanisms of the direction-dependent phase-specific load sharing are discussed in association with the different initial textures in the RD, DD and TD.

Yu, Dunji↗
DOE OSTI2018

Structural Phase Transformation in Strained Monolayer MoWSe 2 Alloy

Two-dimensional (2D) materials exhibit different mechanical properties from their bulk counterparts owing to their monolayer atomic thickness. In this paper, we have examined the mechanical behavior of 2D molybdenum tungsten diselenide (MoWSe 2 ) precipitation alloy grown using chemical vapor deposition and composed of numerous nanoscopic MoSe 2 and WSe 2 regions. Applying a bending strain blue-shifted the MoSe 2 and WSe 2 A 1g Raman modes with the stress concentrated near the precipitate interfaces predominantly affecting the WSe 2 modes. In situ local Raman measurements suggested that the crack propagated primarily thorough MoSe 2 -rich regions in the monolayer alloy. Molecular dynamics (MD) simulations were performed to study crack propagation in an MoSe 2 monolayer containing nanoscopic WSe 2 regions akin to the experiment. Raman spectra calculated from MD trajectories of crack propagation confirmed the emergence of intermediate peaks in the strained monolayer alloy, mirroring experimental results. The simulations revealed that the stress buildup around the crack tip caused an irreversible structural transformation from the 2H to 1T phase both in the MoSe 2 matrix and WSe 2 patches. This was corroborated by high-angle annular dark-field images. Finally, crack branching and subsequent healing of a crack branch were also observed in WSe 2 , indicating the increased toughness and crack propagation resistance of the alloyed 2D MoWSe 2 over the unalloyed counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗
DOE OSTI2020

Design of a V–Ti–Ni alloy with superelastic nano-precipitates

Stress-induced martensitic transformations enable metastable alloys to exhibit enhanced strain hardening capacity, leading to improved formability and toughness. As is well-known from transformation-induced plasticity (TRIP) steels, however, the resulting martensite can limit ductility and fatigue life due to its intrinsic brittleness. In this work, we explore an alloy design strategy that utilizes stress-induced martensitic transformations but does not retain the martensite phase. This strategy is based on the introduction of superelastic nano-precipitates, which exhibit reverse transformation after initial stress-induced forward transformation. To this end, utilizing ab-initio simulations and thermodynamic calculations we designed and produced a V 45 Ti 30 Ni 25 (at%) alloy. In this alloy, TiNi is present as nano-precipitates uniformly distributed within a ductile V-rich base-centered cubic (bcc) β matrix, as well as being present as a larger matrix phase. We characterized the microstructure of the produced alloy using various scanning electron microscopy (SEM) and transmission electron microscopy (TEM) methods. Furthermore, the bulk mechanical properties of the alloy are demonstrated through tensile tests, and the reversible transformation in each of the TiNi morphologies were confirmed by in-situ TEM micro-pillar compression experiments, in-situ high-energy diffraction synchrotron cyclic tensile tests, indentation experiments, and differential scanning calorimetry experiments. The observed transformation pathways and variables impacting phase stability are critically discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

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.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Modeling of Hydraulic Fracture Propagation at the kISMET Site Using a Fully Coupled 3D Network-Flow and Quasi- Static Discrete Element Model

Aimed at supporting the design of hydraulic fracturing experiments at the kISMET site, ~1500 m below ground in a deep mine, we performed pre-experimental hydraulic fracturing simulations in order to estimate the breakdown pressure, propagation pressure, fracture geometry, and the magnitude of induced seismicity using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beams between them. A conjugate 3D flow network based on the DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams. In this paper, the 3D hydraulic fracturing model was verified against the analytic solution for a penny-shaped crack model. We applied the model to simulate fracture propagation from a vertical open borehole based on initial estimates of rock mechanical properties and in-situ stress conditions. The breakdown pressure and propagation pressure are directly obtained from the simulation. In addition, the released elastic strain energies of individual fracturing events were calculated and used as a conservative estimate for the magnitudes of the potential induced seismic activities associated with fracturing. The comparisons between model predictions and experimental results are still ongoing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Development of a gall-resistant stainless-steel hardfacing alloy

Here, this work details the development of a new cobalt-free stainless steel powder metallurgy hardfacing alloy designed to replace Stellite 6, a cobalt-based hardfacing alloy used in nuclear valve applications. The fundamental strategy centers on alloying stainless steels with up to 0.5 wt% nitrogen, which is shown to increase both the volume fraction of hard phase precipitates and the strain-hardening rate of the matrix. The resultant alloy, Nitromaxx, exhibits galling performance that is comparable to Stellite 6, up to 350 °C. This performance is attributed to the suppression of strain localization events associated with galling. In particular, transmission electron microscopy and diffraction measurements from tensile tests show that the nitrogen addition decreases the calculated matrix stacking fault energy and enhances both deformation-induced martensite transformation at room temperature and deformation twinning at elevated temperature. These strain-hardening mechanisms, coupled with the increase in precipitate volume fraction, effectively suppress localization and enhance galling resistance up to 350 °C. Lastly, the enhanced galling resistance cannot be rationalized in terms of tensile stress-strain response alone.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact of Control Blade Insertion on the Deformation Behavior of SiC-SiC Channel Boxes in Boiling Water Reactors

This report describes the analysis of distortion of a Silicon carbide fiber-reinforced silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of temperature and neutron flux dependent irradiation induced swelling of SiC, the SiC-SiC composite channel box can undergo significant distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermalhydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions, to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that the channel box will undergo time dependent lateral bowing for all the control blade positions in the assembly. The bowing behavior is dominated by the swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The lateral bending will cause a temporary interference between the channel box and the control blade. The developed stresses exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the ultimate tensile strength of the material, and therefore, development of a full, through thickness crack in the channel box is not expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Selected results from full-core hydrogen redistribution asymptotic analysis in YH-moderated heat-pipe cooled microreactor

Yttrium hydride is the main candidate for moderation of high temperature nuclear microreactors . This is due to its high thermal stability and relatively high hydrogen retention at temperatures exceeding 870 C. One of the main issues associated with the use of yttrium hydride is that, when exposed to temperature, stress, or concentration gradients, the hydrogen contained in the metallic matrix tends to redistribute and leak from the moderating elements, potentially leading to reactivity losses and power swings. This paper is focused on presenting selected results from the asymptotic hydrogen redistribution analysis performed on the Simplified Microreactor Benchmark Assessment (SiMBA) problem, a microreactor core conceptual design developed at INL. The analysis is based on full core analysis, relies uniquely on Bison to perform hydrogen redistribution calculations, and discusses the physical causes of the reactivity feedback. It was found that the feedback is negative, but it is one order of magnitude lower than the one found for the empire reactor unit-cell. This is due to the lower axial temperature gradient together with the effect of the reflector.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Impact of control blade insertion on the deformation behavior of SiC-SiC channel boxes in BWRs

This article describes the analysis of distortion of a silicon carbide fiber-reinforced, silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of the temperature and neutron flux dependent irradiation-induced swelling of SiC, the SiC-SiC composite channel box can undergo distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermal-hydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that large gradients in fast neutron flux (up to 35–40% across a single axial level) will develop across the channel box. Due to these gradients, the channel box will undergo time-dependent bending for all the control blade positions in the assembly. The time-dependent bowing behavior is dominated by the transient swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The bending will cause temporary interference between the channel box and control blade, and the interference is expected to be most severe for the fully inserted control blade position. The developed stresses due to differential swelling in the channel box exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the tensile strength of the material, and therefore, development of a full, through-thickness crack in the channel box is not expected. Finally, further work is recommended to explore and evaluate the mitigation strategies.

42 ENGINEERING↗
DOE OSTI2020

Deformation accommodating periodic computational domain for a uniform velocity gradient

Many multiscale methods for granular materials use periodic computational domains to consider particle scale interactions and then calculate the stress to drive the continuum scale calculations. For problems involving large material deformations, the computation domain often needs to be reinitialized because of the distortion, causing the loss of the history information of the system. This work introduces an algorithm to accommodate a large deformation of the material while maintaining the computational domain cuboid to avoid the domain reinitialization during the computation. The algorithm uses a rotating frame of reference, in which the velocity gradient can be represented by an upper triangular matrix. The deformation caused by the upper triangular matrix is treated by the image system implied in the periodicity to maintain the computational domain cuboid. The effect from the rotation of the reference frame is considered using the inertial forces. Finally, simulations of simple and pure shear motions are carried out to illustrate the algorithm.

42 ENGINEERING↗
DOE OSTI2019

Implicit renormalization approach to the problem of Cooper instability

In the vast majority of cases, superconducting transition takes place at exponentially low temperature T c out of the Fermi liquid regime. We discuss the problem of determining T c from known system properties at temperatures T >> T c , and stress that this cannot be done reliably by following the standard protocol of solving for the largest eigenvalue of the original gap-function equation. However, within the implicit renormalization approach, the gap-function equation can be used to formulate an alternative eigenvalue problem, where solving leads to an accurate prediction for both T c and the gap function immediately below T c . Furthermore, with diagrammatic Monte Carlo techniques, this eigenvalue problem can be solved without invoking the matrix inversion or even explicitly calculating the four-point vertex function.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗