Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic fields”

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 OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2024

Direction-dependent elastic wave scattering and mode coupling in elastic plates

Asymmetric scattering is a phenomenon in which the field scattered from a discontinuity is dependent on the direction of incidence. In waveguide systems such as elastic plates, the existence of multiple propagating modes provides a platform to explore asymmetric scattering through direction-dependent mode coupling. This paper describes this concept in the context of reciprocal systems and how to utilize it in a general manner.

Haberman, Michael↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2023

A planar defect spin sensor in a two-dimensional material susceptible to strain and electric fields

Abstract The boron-vacancy spin defect ( $${\,{{\mbox{V}}}}_{{{\mbox{B}}}\,}^{-}$$ V B − ) in hexagonal boron nitride (hBN) has a great potential as a quantum sensor in a two-dimensional material that can directly probe various external perturbations in atomic-scale proximity to the quantum sensing layer. Here, we apply first-principles calculations to determine the coupling of the $${\,{{\mbox{V}}}}_{{{\mbox{B}}}\,}^{-}$$ V B − electronic spin to strain and electric fields. Our work unravels the interplay between local piezoelectric and elastic effects contributing to the final response to the electric fields. The theoretical predictions are then used to analyse optically detected magnetic resonance (ODMR) spectra recorded on hBN crystals containing different densities of $${\,{{\mbox{V}}}}_{{{\mbox{B}}}\,}^{-}$$ V B − centres. We prove that the orthorhombic zero-field splitting parameter results from local electric fields produced by surrounding charge defects. This work paves the way towards applications of $${\,{{\mbox{V}}}}_{{{\mbox{B}}}\,}^{-}$$ V B − centres for quantitative electric field imaging and quantum sensing under pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proton Electric Charge Radius from Lepton Scattering

A proton is a bound state of a strong interaction, governed by Quantum Chromodynamics (QCD). The electric charge radius of a proton, denoted by $r^{p}_{E}$, characterizes the spatial distribution of its electric charge carried by the quarks. It is an important input for bound-state Quantum Electrodynamic (QED) calculations of the hydrogen atomic energy levels. However, physicists have been puzzled by the large discrepancy between $r^{p}_{E}$ measurements from muonic hydrogen spectroscopy and those from ep elastic scattering and ordinary hydrogen spectroscopy for over a decade. Tremendous efforts, both theoretical and experimental, have been dedicated to providing various insights into this puzzle, but certain issues still remain unresolved, particularly in the field of lepton scatterings. This review will focus on lepton-scattering measurements of $r^{p}_{E}$, recent theoretical and experimental developments in this field, as well as future experiments using this technique.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Enabling resonant ultrasound spectroscopy in high magnetic fields

Resonant ultrasound spectroscopy (RUS) is a powerful method to determine elastic constants with high accuracy and precision from a single measurement of the mechanical resonances of a sample. Conventionally, the quantitative extraction of elastic moduli with RUS assumes free boundary conditions which can often lead to the adoption of unstable sample positioning between ultrasonic transducers that is incompatible with extreme environments like high magnetic fields. We show that, under specific conditions, introducing a small amount of adhesive between a RUS sample and ultrasonic transducers introduces a perturbation to the free resonance condition which can be accounted for by a simple model. This means elastic constants can be determined to within the uncertainty of conventional RUS, but with significant improvements including sample stability and control of sample orientation. We demonstrate the efficacy of this approach with measurements on a range of materials including room temperature measurements on polycrystalline metals, temperature-dependent measurements of the structural phase transition in strontium titanate single crystals, and magnetic field-dependent measurements of magnetic phase transitions in gadolinium polycrystals up to 14 T.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Development of Aluminum Scandium Nitride Molecular Dynamics Force Fields with Scalable Multi-Objective Bayesian Optimization

Scandium (Sc)-doped aluminum nitride (AlN) exhibits improved piezoelectric properties, which is favorable for sensor applications. Although many experimental studies exist to fine tune the material properties for design purposes, an atomistic level understanding of the structure-property (S-P) relationships is needed, which is the aim of this work. Molecular dynamics can be used to understand the S-P relationships. However, the limited availability of suitable force fields has been a major challenge for accurate property predictions. In this article, a robust force field calibration method using a scalable multi-objective Bayesian optimization approach is presented. Optimizations with three, six, and eight objectives are applied to calibrate aluminum scandium nitride force fields based on the piezoelectric characteristics, modulus of elasticity, and lattice parameters at different Sc-doped levels. Furthermore, the performances of the different force fields are compared, and the performance of the higher dimensional objective problems is discussed. The highly scalable molecular dynamics force field development method is successfully implemented, resulting in the creation of several aluminum scandium nitride molecular dynamics force fields for piezoelectric applications at varying Sc dope levels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lattice dynamics, sound velocities, and atomic environments of szomolnokite at high pressure

Complex mixtures of sulfates, silicates, and ice have been observed in a variety of planetary environments on Earth, Mars and the icy satellites of the solar system. Characterizing the properties of the corresponding compositional endmembers is important for understanding the interiors of a range of planetary bodies in which these phases are observed. To measure the electronic and vibrational properties of the pure ferrous iron endmember of the kieserite group, szomolnokite, (FeSO 4 ∙H 2 O), we have performed synchrotron 57 Fe nuclear resonant inelastic and forward scattering experiments in the diamond-anvil cell up to 14.5 GPa. This pressure range covers depths within Earth’s interior relevant to sulfur cycling in subduction zones and the range of pressures expected within icy satellite interiors. We find evidence of crystal lattice softening, changes in elastic properties, and changes in the electric field gradients of iron atoms associated with two structural transitions occurring within the experimental pressure range. Here we apply these findings to icy satellite interiors, including discussion of elastic properties, modeling of ice-sulfate aggregates, and implications for tidal observations.

58 GEOSCIENCES↗
DOE OSTI2024

Elastic-plastic conductor damage evaluation at over 0.4% strain using a high-stress REBCO coil

Abstract Recent reports on screening current stress simulations of high-field REBCO magnets frequently present peak stresses over 1 GPa. However, this result is probably an unrealistic artifact of purely elastic calculations, considering the macroscopic yield and fracture stresses of approximately 900 MPa and less than 1.1 GPa for Hastelloy substrate-coated conductors. Here, we evaluate elastic-plastic conductor damage at over 0.4% strain using a high-stress REBCO coil exposed to a high field to explore this elastic-plastic regime. The coil was located off-center in a low-temperature superconductor magnet so as to induce a significant screening current in the enhanced radial field. Voltage taps, a Hall sensor, and two strain gauges were used for the instrumentation. We obtained strains exceeding 0.4% near the outward edge during the coil current charge from 350 A to 390 A, where the coil was exposed to external axial and radial magnetic fields of 13 T and 0.5 T. Post mortem results showed wavy plastic deformation, electrical damage, and REBCO defects. An elastic-plastic simulation reproduced the measured strains and predicted that ∼1 GPa stress is sufficient to induce ∼0.9% strain, thus validating our initial concerns with purely elastic models. This paper provides our experimental and simulation results.

Bang, Jeseok (ORCID:0000000247616334)↗
DOE OSTI2023

A guide to the design of the virtual element methods for second- and fourth-order partial differential equations

Here we discuss the design and implementation details of two conforming virtual element methods for the numerical approximation of two partial differential equations that emerge in phase-field modeling of fracture propagation in elastic material. The two partial differential equations are: (i) a linear hyperbolic equation describing the momentum balance and (ii) a fourth-order elliptic equation modeling the damage of the material. Inspired by, we develop a new conforming VEM for the discretization of the two equations, which is implementation-friendly, i.e., different terms can be implemented by exploiting a single projection operator. We use C 0 and C 1 virtual elements for the second-and fourth-order partial differential equation, respectively. For both equations, we review the formulation of the virtual element approximation and discuss the details pertaining the implementation.

42 ENGINEERING↗
DOE OSTI2022

Degradation of Phonons in Disordered Moiré Superlattices

The elastic collective modes of a moiré superlattice arise not from vibrations of a rigid crystal but from the relative displacement between the constituent layers. Despite their similarity to acoustic phonons, these modes, called phasons, are not protected by any conservation law. Furthermore, we show that disorder in the relative orientation between the layers and thermal fluctuations associated with their sliding motion degrade the propagation of sound in the moiré superlattice. Specifically, the phason modes become overdamped at low energies and acquire a finite gap, which displays a universal dependence on the twist-angle variance. Thus, twist-angle inhomogeneity is manifested not only in the non-interacting electronic structure of moiré systems, but also in their phonon-like modes. More broadly, our results have important implications for the electronic properties of twisted moiré systems that are sensitive to the electron-phonon coupling.

2-dimensional systems↗
DOE OSTI2022

Nuclear properties with semilocal momentum-space regularized chiral interactions beyond N 2 LO

We present a comprehensive investigation of few-nucleon systems as well as light and medium-mass nuclei up to A = 48 using the current Low Energy Nuclear Physics International Collaboration two-nucleon interactions in combination with the third-order (N 2 LO) three-nucleon forces. To address the systematic overbinding of nuclei starting from A ≈ 10 found in our earlier study utilizing the N 2 LO two- and three-nucleon forces, we take into account higher-order corrections to the two-nucleon potentials up through fifth order in chiral effective field theory. Here, the resulting Hamiltonian can be completely determined using the A = 3 binding energies and selected nucleon-deuteron cross sections as input. It is then shown to predict other nucleon-deuteron scattering observables and spectra of light p-shell nuclei, for which a detailed correlated truncation error analysis is performed, in agreement with experimental data. Moreover, the predicted ground state energies of nuclei in the oxygen isotopic chain from 14 O to 26 O as well as 40 Ca and 48 Ca show a remarkably good agreement with experimental values, given that the Hamiltonian is fixed completely from the A ≤ 3 data, once the fourth-order (N 3 LO) corrections to the two-nucleon interactions are taken into account. On the other hand, the charge radii are found to be underpredicted by ≈10% for the oxygen isotopes and by almost 20% for 40 Ca and 48 Ca.

20 ≤ A ≤ 38↗
DOE OSTI2023

Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce 2 Zr 2 O 7

The pyrochlore magnet Ce 2 Zr 2 O 7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce 3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce 2 Zr 2 O 7 down to T~0.1K in a magnetic field along the [$1,\bar{1},0$] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [$1,\bar{1},0$] magnetic field at ~2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho 2 Ti 2 O 7 and Dy 2 Ti 2 O 7 , but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce 2 Zr 2 O 7 in [$1,\bar{1},0$] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [$1,\bar{1},0$] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ~3T. For both [$1,\bar{1},0$] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce 2 Zr 2 O 7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. Finally, these results support previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗