Search NASA⌕ Search

SEARCH · Search NASA

Results for “anisotropic elasticity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spin-lattice model for cubic crystals

In this work, we present a methodology based on the Néel model to build a classical spin-lattice Hamiltonian for cubic crystals capable of describing magnetic properties induced by the spin-orbit coupling like magnetocrystalline anisotropy and anisotropic magnetostriction, as well as exchange magnetostriction. Taking advantage of the analytical solutions of the Néel model, we derive theoretical expressions for the parametrization of the exchange integrals and Néel dipole and quadrupole terms that link them to the magnetic properties of the material. This approach allows us to build accurate spin-lattice models with the desired magnetoelastic properties. We also explore a possible way to model the volume dependence of magnetic moment based on the Landau energy. This feature allows us to consider the effects of hydrostatic pressure on the saturation magnetization. We apply this method to develop a spin-lattice model for BCC Fe and FCC Ni, and we show that it accurately reproduces the experimental elastic tensor, magnetocrystalline anisotropy under pressure, anisotropic magnetostrictive coefficients, volume magnetostriction, and saturation magnetization under pressure at zero temperature. This work could constitute a step towards large-scale modeling of magnetoelastic phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗
DOE OSTI2025

Ultrasonic characterization of material heterogeneities in stainless steel components produced by laser powder bed fusion

We introduce pulse-echo ultrasound as a method for characterizing the impact of powder bed fusion parameters on the properties of additively manufactured stainless-steel components, their material anisotropy, and location-dependent heterogeneity. Our results indicate that accurate characterization requires careful selection of ultrasonic propagation paths, which must consider the direction of additive layering, variations in processing parameters, and the component's geometry. We employed two distinct methods to estimate material properties from ultrasonic data: One assumes isotropy, while the other accounts for anisotropic interactions during the propagation of elastic waves. When applied to samples fabricated with laser energy densities ranging from 24 to 42 J/mm³ , these methods revealed transverse isotropy and weak anisotropy (quantified by small Thomsen parameters, ε = 0.0651 and γ = 0.0092) and less than a ∼6 % change in acoustic impedance. The assumption of isotropy, in this case, leads to small errors (less than 4 % or 1 % for Young's modulus in the build or transverse directions) when estimating orthotropic material properties using ultrasonic data measured along just two orthogonal directions, one of which must align with the build direction. By comparing ultrasonic measurements — which aggregate the spatial variability in material properties along the length of elastic wave propagation into a single value — with localized measurements obtained from surface nanoindentation, we uncovered and spatially profiled significant differences between the surface and interior properties. Specifically, the surface Young's modulus decreased from approximately 210 GPa to 180 GPa within a depth of about 3 mm. We attribute this surface-localized heterogeneity in PBF-fabricated components to distinct thermal histories experienced by the surface and interior regions. Collectively, the results of this study establish a framework for the ultrasonic characterization of material heterogeneity and anisotropy in material properties and demonstrate its application in additively manufactured metal components.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polycrystal thermo-elasticity revisited: theory and applications

The self-consistent (SC) theory is the most commonly used mean-field homogenization method to estimate the mechanical response behavior of polycrystals based on the knowledge of the properties and orientation distribution of constituent single-crystal grains. The original elastic SC method can be extended to thermo-elasticity by adding a stress-free strain to an elastic constitutive relation that expresses stress as a linear function of strain. With the addition of this independent term, the problem remains linear. Although the thermo-elastic self-consistent (TESC) model has important theoretical implications for the development of self-consistent homogenization of non-linear polycrystals, in this paper, we focus on TESC applications to actual thermo-elastic problems involving non-cubic (i.e. thermally anisotropic) materials. To achieve this aim, here we provide a thorough description of the TESC theory, which is followed by illustrative examples involving cooling of polycrystalline non-cubic metals. The TESC model allows studying the effect of crystallographic texture and single-crystal elastic and thermal anisotropy on the effective thermo-elastic response of the aggregate and on the internal stresses that develop at the local level.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anisotropic neutron response of trans-stilbene and impact on a handheld dual particle imager

Here this manuscript details measurements of the anisotropic light-output response of trans-stilbene, henceforth referred to as stilbene, along the three primary crystal planes and analyzes how this anisotropic response impacts neutron image reconstruction in a handheld dual particle imager (H2DPI). A prototype H2DPI composed of stilbene pillars (6 x 6 x 50 mm 3 ) coupled to silicon photomultipliers was built and is capable of imaging both fast neutrons and gamma rays from kilogram quantities of special nuclear material. Stilbene was chosen as the scintillating medium for the imager because of its pulse shape discrimination capability, time resolution and relatively high light output. A drawback to using stilbene, however, is the anisotropic response. A recoiling proton from a neutron elastic-scattering event in stilbene will yield different amounts of scintillation light depending on the direction of the recoiling proton with respect to the crystal lattice. This manuscript analyzes how this anisotropic response impacts neutron image reconstruction in 4π. The light output in the three primary crystal planes of stilbene were independently measured 6–7 times for proton recoil energies ranging from 0.5–5.0 MeV using a quasi-monoenergetic time-of-flight neutron source. The measured light-output data in the three crystal planes were fit with a semi-empirical function based on Birks’ formula. These fits were used to reconstruct the location of a 252 Cf source in front of the imager (0°, -0.8°), to the left of the imager (-90°, -0.8°) and directly above the imager (0°, 90°) to determine if it is necessary to apply the directionally-dependent response of stilbene to accurately reconstruct source locations. A bootstrapping technique was applied to the measured data sets to produce 1000 images composed of 2000 cone projections for each source location and each measured light-output curve. List-mode maximum likelihood expectation maximization was applied to each image; the highest pixel location and the full width at half maximum (FWHM) were then recorded. The average and standard deviations of these parameters were taken for each set of images. The maximum difference in the average azimuthal pixel location when varying between the appropriate and incorrect light-output curves for the source locations at (0°, -0.8°) and (-90°, -0.8°) were respectively found to be 1.85 ± 1.39° and 2.03 ± 1.78°. All other pixel locations and FWHMs were within a single standard deviation of uncertainty. Taken together, the anisotropic response of stilbene has negligible impact on the neutron image reconstruction capability of the H2DPI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

An interplay between a hydrogen atmosphere and dislocation characteristics in BCC Fe from time-averaged molecular dynamics

The interplay between hydrogen and dislocations (e.g., core and elastic energies, and dislocation–dislocation interactions) has implications on hydrogen embrittlement but is poorly understood. Continuum models of hydrogen enhanced local plasticity have not considered the effect of hydrogen on dislocation core energies. Energy minimization atomistic simulations can only resolve dislocation core energies in hydrogen-free systems because hydrogen motion is omitted so hydrogen atmosphere formation can’t occur. Additionally, previous studies focused more on face-centered-cubic than body-centered-cubic metals. Discrete dislocation dynamics studies of hydrogen–dislocation interactions assume isotropic elasticity, but the validity of this assumption isn’t understood. Here, we perform time-averaged molecular dynamics simulations to study the effect of hydrogen on dislocation energies in body-centered-cubic iron for several dislocation character angles. We see atmosphere formation and highly converged dislocation energies. We find that hydrogen reduces dislocation core energies but can increase or decrease elastic energies of isolated dislocations and dislocation–dislocation interaction energies depending on character angle. We also find that isotropic elasticity can be well fitted to dislocation energies obtained from simulations if the isotropic elastic constants are not constrained to their anisotropic counterparts. These results are relevant to ongoing efforts in understanding hydrogen embrittlement and provide a foundation for future work in this field.

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

Determining elastic anisotropy of textured polycrystals using resonant ultrasound spectroscopy

Abstract Polycrystalline materials can have complex anisotropic properties depending on their crystallographic texture and crystal structure. In this study, we use resonant ultrasound spectroscopy (RUS) to nondestructively quantify the elastic anisotropy in extruded aluminum alloy 1100-O, an inherently low-anisotropy material. Further, we show that RUS can be used to indirectly provide a description of the material’s texture, which in the present case is found to be transversely isotropic. By determining the entire elastic tensor, we can identify the level and orientation of the anisotropy originated during extrusion. The relative anisotropy of the compressive (c 11 /c 33 ) and shear (c 44 /c 66 ) elastic constants is 1.5% ± 0.5% and 5.7% ± 0.5%, respectively, where the elastic constants (five independent elastic constants for transversely isotropic) are those associated with the extrusion axis that defines the symmetry of the texture. These results indicate that the texture is expected to have transversely isotropic symmetry. This finding is confirmed by two additional approaches. First, we confirm elastic constants and the degree of elastic anisotropy by direct sound velocity measurements using ultrasonic pulse echo. Second, neutron diffraction (ND) data confirm the symmetry of the bulk texture consistent with extrusion-induced anisotropy, and polycrystal elasticity simulations using the elastic self-consistent model with input from ND textures and aluminum single-crystal elastic constants render similar levels of polycrystal elastic anisotropy to those measured by RUS. We demonstrate the ability of RUS to detect texture-induced anisotropy in inherently low-anisotropy materials. Therefore, as many other common materials have intrinsically higher elastic anisotropy, this technique should be applicable for similar levels of texture, providing an efficient general diagnostic and characterization tool.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Verification of stability and unraveling the electronic and physical properties of bulk and (001)-surfaces of newly synthesized Ti 2 ZnX (X = C, N) MAX phases

MAX phase family has been extended by the addition of late transition metals at the A-site with the expectation of diverse functional properties. Here, we present our systematic density functional investigation on the thermodynamic and phonon stabilities, elastic properties, including elastic constants, elastic moduli and elastic anisotropy of newly synthesized Ti 2 ZnX (X = C, N) phases in comparison with conventional Ti 2 AlX (X = C, N). Due to the smaller size of N as compared to C, the unit cell dimension is reduced when C atoms are replaced by N atoms at the X-site. Furthermore, the Ti 2 ZnC and Ti 2 ZnN are stable at the equilibrium volume of 110.84 Å 3 and 105.70 Å 3 . The thermodynamic, mechanical and dynamical stabilities are validated by estimating the formation energies, elastic constants and phonon dispersions, respectively. The elastic properties of Ti 2 ZnN are less anisotropic as compared to those of Ti 2 ZnC. To understand the thin-film characteristics in Ti 2 ZnX, the surface properties with (001)-terminated slabs are investigated. Both Ti 2 ZnX bulk and (001)-surfaces exhibit metal-like electronic structures. There is a strong covalent bonding between Ti-X and Ti-Zn atoms confirmed by the charge density map and Mulliken population analysis. Additional states are generated at the Fermi level (EF) due to the unusual d-p states hybridization between Ti and Zn atoms. The anisotropy in chemical bonding is confirmed by the cleavage energy difference between Ti-X and Ti-Zn. Here, Ti(X)-001 and Zn-001 terminations are stable surfaces; however, in terms of chemical potentials, Zn-001 termination is the most favourable in Ti 2 ZnX.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Topology control of human fibroblast cells monolayer by liquid crystal elastomer

Eukaryotic cells in living tissues form dynamic patterns with spatially varying orientational order that affects important physiological processes such as apoptosis and cell migration. The challenge is how to impart a predesigned map of orientational order onto a growing tissue. Here, we demonstrate an approach to produce cell monolayers of human dermal fibroblasts with predesigned orientational patterns and topological defects using a photoaligned liquid crystal elastomer (LCE) that swells anisotropically in an aqueous medium. The patterns inscribed into the LCE are replicated by the tissue monolayer and cause a strong spatial variation of cells phenotype, their surface density, and number density fluctuations. Unbinding dynamics of defect pairs intrinsic to active matter is suppressed by anisotropic surface anchoring allowing the estimation of the elastic characteristics of the tissues. The demonstrated patterned LCE approach has potential to control the collective behavior of cells in living tissues, cell differentiation, and tissue morphogenesis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Water-Fe-Pressure dependent single-crystal elastic properties of wadsleyite: Implications for the seismic anisotropy in the upper Mantle Transition Zone

Recent seismic studies suggest an anisotropic Mantle Transition Zone (MTZ) in areas adjacent to subducted slabs. Wadsleyite is the main anisotropy contributor in the upper MTZ, therefore the interpretation of these seismic observations requires the knowledge of single-crystal elastic moduli (C ij s) and the deformation-induced lattice preferred orientation (LPO) of wadsleyite. Wadsleyite can host up to 3 wt% water in its crystal structure as point defects in the form of hydroxyl groups, however, the combined effect of water content, Fe content, and pressure on the C ij s of wadsleyite remains unclear. In this study, we measured the high-pressure single-crystal elasticity of a synthetic hydrous Fe-bearing wadsleyite (0.14(5) wt% water, Fe#=9.4, Fe 3+ /ΣFe=0.3) up to 18.2 (2) GPa. In combination with previous experimental data, we separated the effects of pressure, water, and Fe contents on the Cijs and intrinsic elastic anisotropy of wadsleyite. Overall, our results suggest that the intrinsic elastic anisotropy of wadsleyite decreases with pressure, water, and Fe contents. At 15 GPa, increasing the water content by 0.1 wt% or Fe# by 1 decreases the V P and V S anisotropy of wadsleyite by ~1.1-1.3%, and ~0.8-1.3% in average, respectively. Combined with the LPO determined in previous deformation experiments, we modeled the seismic anisotropy in the upper MTZ generated by a sub-vertical mantle flow near cold subducted slabs and a sub-horizontal mantle flow in the ambient mantle. In both scenarios, the LPO of wadsleyite leads to V SV (vertically polarized shear wave velocity) > V SH (horizontally polarized shear wave velocity). Our results suggest that wadsleyite may account for a weak anisotropic MTZ (<1%) on the global scale. Considering the fact that water decreases the elastic anisotropy but promotes LPO of wadsleyite, seismic anisotropy may not be a good water sensor.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fabric-enriched continuum breakage mechanics (F-CBM )

A fabric-enriched continuum breakage mechanics (F-CBM) framework is developed to examine the relation between elastic and inelastic anisotropy in granular materials subjected to high pressure. For this purpose, emphasis is given to the influence of the anisotropic fabric on the plastic cap regime – that is, a macroscopic signature of the inelasticity of granular materials widely dependent on the particle crushability. Anisotropic behaviour prior to yielding is introduced through a symmetric second-order fabric tensor embedded in the expression of the elastic energy potential. It is shown that, thanks to its energy-based formulation, the F-CBM framework provides a platform to predict the rotation and distortion of the yield cap of granular materials as an outcome of anisotropic strain energy storage prior to grain crushing. Parametric analyses indicate that both the degree of anisotropy and the orientation of the fabric have an impact on the energy release resulting from grain crushing and contribute to the resulting stress–strain behaviour. It is shown that, despite the use of very few parameters, the proposed framework can accurately predict the anisotropic yielding and inelastic stress–strain response solely on the basis of the directional properties inferred from measurements of elastic anisotropy.

Engineering↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Sliding dynamics for bubble phases on periodic modulated substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short-range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Sliding Dynamics for Bubble Phases on Periodic Modulated Substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold.

36 MATERIALS SCIENCE↗