Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic materials”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Spontaneous Twirls and Structural Frustration in Moiré Materials

Structural twirls form spontaneously in the domain wall networks of some moiré materials. We show that in heterobilayers, neighboring twirl chiralities tend to antialign, forming staggered patterns that are well described by antiferromagnetic lattice 𝜙 4 theories. In moiré systems with triangular domains, this leads to frustration in the chirality configuration of the twirls and to hysteresis with respect to variation of the average twist angle and possibly other control parameters. As a result, we expect that in typical materials, the ordering temperature of twirls is about 10 3 K, and that thermal fluctuations in individual twirl chiralities freeze below room temperature.

Elastic deformation↗
DOE OSTI2022

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗
DOE OSTI2021

The Transport Properties of Quasi–One-Dimensional Ba 3 Co 2 O 6 (CO 3 ) 0.7

We have performed combined elastic neutron diffuse, electrical transport, specific heat, and thermal conductivity measurements on the quasi–one-dimensional Ba 3 Co 2 O 6 (CO 3 ) 0.7 single crystal to characterize its transport properties. A modulated superstructure of polyatomic CO 3 2- is formed, which not only interferes the electronic properties of this compound, but also reduces the thermal conductivity along the c-axis. Furthermore, a large magnetic entropy is observed to be contributed to the heat conduction. Our investigations reveal the influence of both structural and magnetic effects on its transport properties and suggest a theoretical improvement on the thermoelectric materials by building up superlattice with conducting ionic group.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Static adhesion hysteresis in elastic structures

Adhesive interactions between elastic structures such as graphene sheets, carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to irrecoverable energy loss associated with bond breakage, even in static (rate-independent) experiments. To understand this phenomenon, we start with a minimal theory for the peeling of a thin sheet from a substrate, coupling the local event of bond breaking to the nonlocal elastic relaxation of the sheet and show that this can drive static adhesion hysteresis over a bonding/debonding cycle. Using this model we quantify hysteresis in terms of the adhesion and elasticity parameters of the system. We find that this allows us to derive a scaling relation that preserves hysteresis at different levels of granularity while resolving a seeming paradox of lattice trapping in the continuum limit of a discrete fracture process. Finally, to verify our theory, we use new experiments to demonstrate and measure adhesion hysteresis in bundled microtubules.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of tilt grain boundaries on the structural integrity of WSe 2 monolayers

Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe 2 ) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe 2 , of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe 2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe 2 monolayer with tilted GBs. WSe 2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures.

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

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Responsive‐Hydrogel Aquabots

Abstract It remains a challenge to produce soft robots that can mimic the responsive adaptability of living organisms. Rather than fabricating soft robots from bulk hydrogels,hydrogels are integrated into the interfacial assembly of aqueous two‐phase systems to generate ultra‐soft and elastic all‐aqueous aquabots that exhibit responsive adaptability, that can shrink on demand and have electrically conductive functions. The adaptive functions of the aquabots provide a new platform to develop minimally invasive surgical devices, targeted drug delivery systems, and flexible electronic sensors and actuators.

42 ENGINEERING↗
DOE OSTI2022

Magnetoelastic properties of multiferroic hexagonal ErMnO 3

Here, the strength and dynamics of magnetoelastic coupling through the paramagnetic (PM) – antiferromagnetic (AFM) – ferrimagnetic (FIM) transitions in multiferroic hexagonal ErMnO 3 have been investigated by Resonant Ultrasound Spectroscopy. Elastic stiffening by up to 2% below the PM – AFM transition at 80 K arises from biquadratic coupling between strain and the magnetic order parameter with relaxation times longer than ~10 -6 s for the response of spins to changes in strain. In contrast with YMnO 3 , the PM – AFM transition in ErMnO 3 is accompanied by a peak in acoustic loss immediately below the Néel point which is interpreted in terms of strain relaxation accompanying ordering of spins of Er 3+ at 4b sites. Changes in the magnetic ordering scheme at the AFM – FIM transition near 3 K are accompanied by elastic softening of ~0.03 %. During poling of the low temperature ferrimagnetic structure round magnetic hysteresis loops, small changes in elastic stiffness which arise due to the contribution of piezomagnetic and/or piezoelectric moduli are detected. Contributions of piezoelectric moduli to acoustic resonance frequencies also permit changes in the configuration of ferroelectric domains to be detected in response both to cycling through this transition and to application of a magnetic field. A peak in acoustic loss in the vicinity of 250 K is attributed to strain-mediated pinning/freezing of some aspect of the domain microstructure with an activation energy of ~0.25–0.3 eV. A return to the original elastic properties on heating to temperatures above ~250 K is interpreted in terms of backswitching of domains to the configuration they had at the start. These observations confirm the existence of subtle variations in magnetoelastic coupling behaviour relating to both the magnetic order parameters and magnetic domain structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
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

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geometric compatibility measure m' for twin transmission: A predictor or descriptor?

In this work, the geometric compatibility factor m' is critically analyzed to assess whether it can be used to interpret/predict twin transmission (TT) across grain boundaries (GBs). This geometric measure is widely used to relate the likelihood of TT to the misalignment of both the shear and plane-normal directions within a twin set (i.e., incoming and outgoing twin). Here, using a large set of electron back scattering diffraction (EBSD) data, a detailed statistical analysis of twin-GB interactions is performed for {${1\bar{01}}2$} tensile twins in hexagonal close-packed (HCP) metals Mg, Zr, and Ti at different strain levels. In addition, a full-field crystal plasticity model is employed to quantify the role of local stresses and the applicability of m' as a criterion for the TT process. This combined study addresses the following three main questions: (i) What is the fidelity of m' in describing experimentally observed TTs? (ii) Can m' be used as a metric to predict/anticipate TT? (iii) Does m' naturally capture local stress effects? As a descriptor, m' cannot rationalize ~25% of TT events observed in Mg or more than 50% of TT events in Zr and Ti. As a predictor, the m'-measure does not predict TT events in over ~50% of twin-GB interactions analyzed. Further, the applicability of m' to describe and predict TT events decreases with an increase in elastic anisotropy, plastic anisotropy, and macroscopic strain levels. Finally, the twinning simulations reveal that m' does not capture the key effects of local stresses on variant selection upon twin transmission. The local stress induced by the twinning shear transformation plays a dominant role in driving the TT process compared to the geometric alignment of the constituting twins, i.e., m'.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boundary conditions manipulation of polar vortex domains in BiFeO 3 membranes via phase-field simulations

Polar vortex domains have recently become an emergent research field due to the abundant physical phenomena and potential applications in high-density memories. Here, we explore the mechanisms of creating polar vortex domains in the BiFeO3 (BFO) membranes subjected to different boundary conditions using phase-field simulations. A major difference is that the vortex in membrane can be stabilized even under short-circuit electrical boundary conditions compared to vortex in other systems, such as thin film or superlattice. We found that (a) the formation of polar vortex domains at the membrane interior under bending is mainly driven by the reduction of elastic energy under short-circuit boundary condition, and the vortex chirality (namely, clockwise and counterclockwise) could be identified by n-shape and u-shape bending; (b) in the unbent open-circuit BFO membrane case, exotic trapezoid-shaped vortex nanodomains form at the terminations of 109 degrees domain walls (DWs) and partially charged 71 degrees DWs, which is driven by the local depolarization field and the interplay among electrostatic, elastic, and gradient and Landau energies. We also examine Kittel's law by establishing the dependence of vortex periods on the membrane thickness. Furthermore, these results give further understanding of the effect of boundary conditions on the formation of polar vortex domains, guiding experimental designs of vortex-based high-density memories.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Conservation of angular momentum in an elastic medium with spins

Exact conservation of the angular momentum is worked out for an elastic medium with spins. The intrinsic anharmonicity of the elastic theory is shown to be crucial for conserving the total momentum. As a result, any spin-lattice dynamics inevitably involves multiphonon processes and interaction between phonons. Furthermore, this makes transitions between spin states in a solid fundamentally different from transitions between atomic states in a vacuum governed by linear electrodynamics. Consequences for using solid-state spins as qubits are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

WHA Sphere Impacting RHA Plate: Modeling Brittle Behavior in Metal Alloys

Brittle behavior of metal alloys is often critical to modeling ballistic impact and penetration. The ALEGRA multiphysics finite element software incorporates calibrated models for the equation of state, elasticity, yield stress, plasticity and fracture, but simulations do not always capture expected metal fracture. Here we report concerted efforts to do so for one important case where experiments clearly show shear fractures: a tungsten sphere impacting a steel plate at various angles. Our best simulations show fractures that are qualitatively similar to experiments, but there are significant differences in quantitative metrics. Specifically, velocities of tracers used to quantify simulated plug parameters consistently fall short of measured plug velocities. Also, simulated plugs break apart more than expected from experimental evidence. We attribute these shortfalls to the lack of an explicit shear fracture mechanism in the material models, leading to over-estimated resistance to plug formation and movement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Low-phase-noise surface-acoustic-wave oscillator using an edge mode of a phononic band gap

Low-phase-noise microwave-frequency integrated oscillators provide compact solutions for various applications in signal processing, communications, and sensing. Surface acoustic waves (SAWs), featuring orders-of-magnitude shorter wavelength than electromagnetic waves at the same frequency, enable integrated microwave-frequency systems with much smaller footprint on chip. SAW devices also allow higher-quality (Q) factors than electronic components at room temperature. Here, we demonstrate a low-phase-noise gigahertz-frequency SAW oscillator on 128°Y-cut lithium niobate, where the SAW resonator occupies a footprint of 0.05 mm 2 . Leveraging phononic crystal band-gap-edge modes to balance between Q factors and insertion losses, our 1-GHz SAW oscillator features a low phase noise of -132.5 dBc/Hz at a 10-kHz offset frequency and an overlapping Hadamard deviation of 6.5 × 10 -10 at an analysis time of 64 ms. The SAW resonator-based oscillator holds high potential in developing low-noise sensors and acousto-optic integrated circuits.

42 ENGINEERING↗