Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic crystal”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗
NASA NTRS2004

Subsurface Stress Fields In Single Crystal (Anisotropic) Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent HCF failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is , presented, for evaluating the subsurface stresses in the elastic half-space, using a complex potential method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis. Effects of crystal orientation on stress response and fatigue life are examined.

Arakere, Nagaraj K.↗
NASA NTRS2003

Subsurface Stress Fields in Single Crystal (Anisotropic) Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and fatigue stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent HCF failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. Techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts are presented in this report. Figure 1 shows typical damper contact locations in a turbine blade. The subsurface stress results are used for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades.

Arakere, Nagaraj K.↗
DOE OSTI2026

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗
NASA NTRS2002

Generalized Reliability Methodology Applied to Brittle Anisotropic Single Crystals

A generalized reliability model was developed for use in the design of structural components made from brittle, homogeneous anisotropic materials such as single crystals. The model is based on the Weibull distribution and incorporates a variable strength distribution and any equivalent stress failure criteria. In addition to the reliability model, an energy based failure criterion for elastically anisotropic materials was formulated. The model is different from typical Weibull-based models in that it accounts for strength anisotropy arising from fracture toughness anisotropy and thereby allows for strength and reliability predictions of brittle, anisotropic single crystals subjected to multiaxial stresses. The model is also applicable to elastically isotropic materials exhibiting strength anisotropy due to an anisotropic distribution of flaws. In order to develop and experimentally verify the model, the uniaxial and biaxial strengths of a single crystal nickel aluminide were measured. The uniaxial strengths of the <100> and <110> crystal directions were measured in three and four-point flexure. The biaxial strength was measured by subjecting <100> plates to a uniform pressure in a test apparatus that was developed and experimentally verified. The biaxial strengths of the single crystal plates were estimated by extending and verifying the displacement solution for a circular, anisotropic plate to the case of a variable radius and thickness. The best correlation between the experimental strength data and the model predictions occurred when an anisotropic stress analysis was combined with the normal stress criterion and the strength parameters associated with the <110> crystal direction.

Salem, Jonathan A.↗
NASA NTRS1984

The effect of anisotropic crystal-melt surface tension on grain boundary groove morphology

The shape of a stationary solid-liquid interface in a temperature gradient near a grain boundary in a pure material is calculated for anisotropic crystal-melt surface tension and equal thermal conductivities of crystal and melt. Results are compared with those for the well-known problem of the two-dimensional equilibrium shape of a crystal. For small anisotropy, the resulting interface shapes have continuously turning tangents but differ in detail from the grain boundary groove shapes that have been calculated for isotropic surface tension. For larger anisotropy, the interface shapes have discontinuities in slope as a result of missing orientations.

Voorhees, P. W.↗
DOE OSTI2024

Peak effect and dynamics of stripe- and pattern-forming systems on a periodic one-dimensional substrate

Here we examine the ordering, pinning, and dynamics of two-dimensional pattern-forming systems interacting with a periodic one-dimensional substrate. In the absence of the substrate, particles with competing long-range repulsion and short-range attraction form anisotropic crystal, stripe, and bubble states. When the system is tuned across the stripe transition in the presence of a substrate, we find that there is a peak effect in the critical depinning force when the stripes align and become commensurate with the substrate. Under an applied drive, the anisotropic crystal and stripe states can exhibit soliton depinning and plastic flow. When the stripes depin plastically, they dynamically reorder into a moving stripe state that is perpendicular to the substrate trough direction. We also find that when the substrate spacing is smaller than the widths of the bubbles or stripes, the system forms pinned stripe states that are perpendicular to the substrate trough direction. The system exhibits multiple reentrant pinning effects as a function of increasing attraction, with the anisotropic crystal and large bubble states experiencing weak pinning but the stripe and smaller bubble states showing stronger pinning. We map out the different dynamic phases as a function of filling, the strength of the attractive interaction term, the substrate strength, and the drive, and demonstrate that the different phases produce identifiable features in the transport curves and particle orderings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Local Thermal Conductivity Patterning in Rotating Lattice Crystals of Anisotropic Sb 2 S 3

The ability to control material heat transport properties over space and time can drive advanced functionalities in thermal management for electronics and system-on-chip, and enable thermal circuits. Despite the technological relevance, there are limited demonstrations of local thermal property control. Rotating lattice single (RLS) crystals—formed via laser-induced crystallization of an amorphous substrate—offer a novel avenue for local crystal engineering, unlocking opportunities for microscale property patterning. Here, thermal conductivity (𝜅) imaging is applied to RLS crystals of Sb2S3 to resolve microscale 𝜅 variations across patterned regions. Amorphous areas exhibit 𝜅 as low as 0.6 Wm −1 K −1 , while crystalline regions display periodic 𝜅 variations from 0.7 to over 2.5 Wm −1 K −1 . These variations correspond to changes in crystal orientation, revealing marked 𝜅 anisotropy. The crystal out-of-plane direction (c axis)—featuring van der Waals bonds—shows amorphous-like transport, whereas in-plane directions (a, b axes) exhibit 3.5x and 1.7x larger 𝜅, respectively. First-principles calculations, in excellent agreement with experiments, suggest that the in-plane anisotropy originates from expressed Sb lone pairs, which impart a corrugation along the b axis affecting bond stiffness and 𝜅. These findings demonstrate microscale control of thermal properties via laser-processed metastructures, with significant implications for next-generation thermal management.

14 SOLAR ENERGY↗
DOE OSTI2026

i -Wave Symmetry Altermagnetism and Anomalous Hall Conductivity in Monolayer FeCl 3

In this work, the FeCl 3 monolayer is proposed to host a rare and exotic i -wave symmetry altermagnetism. Using first-principles calculations, we performed a detailed investigation of the structural stability and electronic, magnetic, and spintronic properties of monolayer FeCl 3 . Our results suggest the system is semiconducting in nature and exhibits i-wave symmetry nonrelativistic spin splitting (NRSS), originating from its hexagonal crystal structure with roto-inversion symmetries and anisotropic crystal field. The momentum-resolved spin-splitting further supports the presence of spin-dependent topological features. Further, the monolayer exhibits sizable conventional spin-Hall and anomalous Hall conductivities, highlighting its promise as a versatile platform for next-generation flexo-spintronic applications.

altermagnetism↗
DOE OSTI2025

Daily modulation of low-energy nuclear recoils from sub-GeV dark matter

At sufficiently low nuclear recoil energy, the scattering of dark matter (DM) in crystals gives rise to single phonon and multiphonon excitations. In anisotropic crystals, the scattering rate into phonons modulates over each sidereal day as the crystal rotates with respect to the DM wind. This gives a potential avenue for directional detection of DM. The daily modulation for single phonons has previously been calculated. Here we calculate the daily modulation for multiphonon excitations from DM in the mass range 1 MeV–1 GeV. We generalize previous multiphonon calculations, which made an isotropic approximation, and implement results in the DarkELF package. We find daily modulation rates up to 1%–10% for an Al 2 O 3 target and DM mass below 30 MeV, depending on the recoil energies probed. We obtain similar results for SiC, while modulation in Si, GaAs, and SiO 2 is negligible. Published by the American Physical Society 2025

Stratman, Connor↗
DOE OSTI2026

Coarse-Grained Simulations of Polyrotaxane Hydrogels under Quiescent and Shear Conditions

Cyclodextrin-based polyrotaxanes (PR) form hydrogels in water through cyclodextrin (CD) aggregation and crystallization. These networks can break under shear flows, making them versatile platforms for extrusion-based 3D printing. To optimize the material properties of 3D-printed PR gels, a microscopic understanding of the structural evolution during 3D printing is necessary. Here, we employ coarse-grained (CG) simulations to reveal the PR assembly process at the molecular level. Our simulations reproduce the experimental crystal morphologies of PR at varying concentrations and chain lengths under quiescent conditions. Using nonequilibrium simulations, we show shear flow ruptures crystalline domain connectivity in PR gels and stacks the lamellae in the gradient direction, allowing the materials to flow during 3D printing. After the cessation of flow, the anisotropic crystal alignment and absence of available dangling PR diminish the intercrystal connectivity. The printed materials are therefore mechanically weaker than the pristine hydrogels, in agreement with experimental results. Nonetheless, by relating the microscopic structural evolution with viscoelastic properties of PR gels and solutions, we elucidate how flow conditions and sample composition affect the 3D printing performance of PR hydrogels.

Smith, Cameron D. [Dartmouth College, Hanover, NH ↗
DOE OSTI2025

Directional locking and hysteresis in stripe- and bubble-forming systems on one-dimensional periodic substrates with a rotating drive

We examine the dynamics of a two-dimensional stripe, bubble, and crystal forming system interacting with a periodic one-dimensional substrate under an applied drive that is rotated with respect to the substrate periodicity direction 𝑥. We find that the stripes remain strongly directionally locked to the 𝑥 direction for an extended range of drives before undergoing motion parallel to the drive. In some cases, the stripes break apart at the unlocking transition, but they can dynamically reform into stripes aligned perpendicular to the 𝑥 direction, producing hysteresis in the directional locking and unlocking transitions. In contrast, moving anisotropic crystal and bubble phases exhibit weaker directional locking and reduced or no hysteresis. The hysteresis occurs in regimes where the particle rearrangements occur and is most pronounced near the stripe phase. In conclusion, we also show that for varied substrate strength, substrate spacing, and particle density, a number of novel dynamical patterns can form that include a combination of stripe, bubble, and crystal morphologies.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Ice crystal growth in water vapor at high saturation

A simple technique is presented for estimating the energy of formation of monolayer icelike clusters at ice-vapor interfaces. Under the assumptions that the ice surfaces are smooth and sparsely covered with monomers, dimers, etc., in near equilibrium with the vapor, and that the bond energies and configurational entropy dominate the energy of formation, it is found that the basal surfaces prefer triangular embryos with an orientation which reverses from layer to layer, whereas the most stable clusters on the prism surfaces are rectangular in configuration. The preferred prism clusters are determined to have a significantly lower critical energy of formation than the basal clusters due to differences in both corner free energy and configurational entropy. This phenomenon provides a mechanism for strongly anisotropic crystal growth at high saturations.

Bartley, D. L.↗
DOE OSTI2025

Directed Self-Assembly of the Organic Semiconductor C8-BTBT-C8 in Anodic Aluminum Oxide Nanopores

Controlling the self-assembly of organic semiconductors at the nanoscale is critical for advancing high-performance electronic and photonic devices, yet it remains challenging due to their intrinsic anisotropic crystallization and sensitivity to processing conditions. Here, in this study, we demonstrate that cylindrical nanoconfinement within anodic aluminum oxide membranes provides a versatile platform to precisely tune the molecular orientation and phase behavior of the prototypical organic semiconductor 2,7-dioctyl[1]benzothieno[3,2-b][1]benzothiophene (C8-BTBT-C8). Combining temperature-dependent high-resolution synchrotron X-ray scattering with optical birefringence measurements, we uncover that confinement geometries (pore diameters 25–180 nm) and surface chemistry govern the emergence of distinct smectic A textures, featuring molecular layers either parallel or perpendicular to the pore axis. The competition between axial and radial smectic layering is modulated by pore size, surface hydrophilicity, and thermal history, enabling reversible control over domain orientations and transitions between liquid crystalline and crystalline states. Notably, nanoconfinement stabilizes the smectic phase over an expanded temperature range compared to bulk, while inducing complex multidomain configurations owing to geometric constraints and anchoring conditions. Our results elucidate fundamental mechanisms by which anisotropic nanoscale confinement directs the self-organization of highly conjugated organic molecules, with implications for optimizing directional charge transport and anisotropic optical responses in organic–inorganic hybrid nanoarchitectures. This study establishes nanoconfinement as a powerful strategy to engineer morphology and functional properties in organic semiconducting materials with nanoscale precision.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Grain boundary self-diffusion and point defect interactions in α -U via molecular dynamics

Though metallic U-Zr fuel has been used in nuclear reactors since the 1960s, many of its fundamental and thermodynamic properties are still unknown. The a-U phase, which has a highly anisotropic crystal structure and physical properties, is present in U-Zr fuel. The character and behavior of a-U grain boundaries will strongly impact fuel thermophysical performance under irradiation. Here, we study the interaction of point defects with grain boundaries, diffusion along grain boundaries, and the predicted diffusional creep behavior of a-U via molecular dynamics. We calculate the segregation energy of vacancies and interstitials to grain boundaries and quantify the biased sink strength of the grain boundaries, and observe that this sink strength is not strongly dependent on the grain boundary orientation. We also find that grain boundary diffusivity is strongly dependent on the grain boundary energy and grain boundary orientation. The presence of point defects within the grain boundary can induce diffusion in grain boundaries with low formation energies and can enhance diffusion in high-energy grain boundaries. We also find that diffusional creep of a-U at prototypical metallic fuel operation conditions is extremely high and could help explain observed metallic fuel swelling behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗