SEARCH · Search NASA

Results for “deformation function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Predictive machine learning approaches for the microstructural behavior of multiphase zirconium alloys

Abstract Zirconium alloys are widely used in harsh environments characterized by high temperatures, corrosivity, and radiation exposure. These alloys, which have a hexagonal closed packed (h.c.p.) structure thermo-mechanically degrade, when exposed to severe operating environments due to hydride formation. These hydrides have a different crystalline structure, than the matrix, which results in a multiphase alloy. To accurately model these materials at the relevant physical scale, it is necessary to fully characterize them based on a microstructural fingerprint, which is defined here as a combination of features that include hydride geometry, parent and hydride texture and crystalline structure of these multiphase alloys. Hence, this investigation will develop a reduced order modeling approach, where this microstructural fingerprint is used to predict critical fracture stress levels that are physically consistent with microstructural deformation and fracture modes. Machine Learning (ML) methodologies based on Gaussian Process Regression, random forests, and multilayer perceptrons (MLP) were used to predict material fracture critical stress states. MLPs, or neural networks, had the highest accuracy on held-out test sets across three predetermined strain levels of interest. Hydride orientation, grain orientation or texture, and hydride volume fraction had the greatest effect on critical fracture stress levels and had partial dependencies that were highly significant, and in comparison hydride length and hydride spacing have less effects on fracture stresses. Furthermore, these models were also used accurately predicted material response to nominal applied strains as a function of the microstructural fingerprint.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effect of superconductivity on the shape of flat bands

Abstract For the first time, basing both on experimental facts and our theoretical consideration, we show that Fermi systems with flat bands should be tuned with the superconducting state. Experimental measurements on magic-angle twisted bilayer graphene of the Fermi velocity V F as a function of the temperature T c of superconduction phase transition have revealed , where N s (0) is the density of states at the Fermi level. We show that the high- T c compounds exhibit the same behavior. Such observation is a challenge to theories of high- T c superconductivity, since V F is negatively correlated with T c , for . We show that the theoretical idea of forming flat bands in strongly correlated Fermi systems can explain this behavior and other experimental data collected on both and twisted bilayer graphene. Our findings place stringent constraints on theories describing the nature of high- T c superconductivity and the deformation of flat band by the superconducting phase transition.

Physics↗

DOE OSTI2017

The Tunisian Jurassic aquifer in the North African Sahara aquifer system: information derived from two-dimensional seismic reflection and well logs

Southern Tunisia is an arid area where socio-economic activities are dependent on groundwater resources. The presented study aims to better characterize the Jurassic aquifer based on geological and geophysical data, with a view to develop a rational exploitation program. Well logs are used to precisely determine the position and composition of the known Jurassic aquifer layers and to identify others able to produce good quality water. The logs show that limestones, sandstones and dolomites of the Krachoua, Techout and Foum Tataouine formations are the main Jurassic aquifers. Sixty-eight seismic-reflection sections are integrated within this study. The interpolation between the interpreted sections leads to the construction of isochronous isopach maps and geoseismic sections, and their analysis finds that compressive and extensive tectonic deformations have influenced the Jurassic aquifer geometry. The Hercynian orogeny phase manifestation is remarkable in that there are several stratigraphic gaps in the Jurassic sequence. The E–W, NW–SE, and NNW–SSE accidents, reactivated in normal faults since the Permian to Lower Cretaceous epochs, have generated the structures found in the Jurassic series, such as subsided and raised blocks. Their syn-sedimentary activity has controlled the thickness and facies of these series. The Cretaceous, Tortonian and Post-Villafranchian compressions are responsible for the Jurassic-deposits folding in some localities. The highlighted tectonic and sedimentary events have an important impact on the Jurassic aquifer function by favoring the Jurassic aquifer interconnections and their connections with the Triassic and Cretaceous permeable series.

58 GEOSCIENCES↗

DOE OSTI2025

Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides

Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Measurement of crack length in width tapered beam experiments

The width tapered beam method for measuring fracture toughness has been contemplated for use in measuring the adhesion of photovoltaic materials, is promoted by it being a viable method for use directly on modules rather than to engineered test specimens. However, precise determination of crack length is often difficult as failure often involves tendril formation, cavitation voids, and other plastic deformation at the crack front in addition to obfuscation by opaque materials. Because fracture energy varies as the square of crack length, imprecise crack tip length results in large measurement uncertainties. In this work, we develop a method to estimate the crack tip length using the relationship between the applied force and the displacement of the load frame. The crack tip is determined for every data point and can be defined from a calibration curve related to the amount of energy put into the beam, or from a theoretically idealized estimation. Once a calibration curve is set up and the equations programmed, the fracture energy as a function of position on the newly formed crack interface can be accurately estimated. The test additionally provides much easier analysis and insight into the joint failure through more precise knowledge of fracture toughness at each de-bonded area.

42 ENGINEERING↗

DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

DOE OSTI2026

Linking Domain Structure Evolution at a Grain Boundary to Piezoelectric Response via Nano‐Diffraction

Electric-field-induced domain structure switching in a 1 𝜇⁢m -thick Pb 0.99 (Zr 0.45 Ti 0.55 ) 0.98 Nb 0.02 O 3 (Nb-doped PZT) bicrystal film was characterized in situ via nano-focused synchrotron diffraction. The epitaxial film was deposited on a (100) SrTiO 3 bicrystal substrate. The changes in domain structure were mapped within a 5 𝜇⁢m × 5 𝜇⁢m area at and around an in-plane tilt-type (23.6°) grain boundary with 50 nm spatial resolution. Rocking curves collected at each point of the mapped area provided the ability to reconstruct spatially-varying three-dimensional (3D) reciprocal space maps around the grain boundary. The 3D reciprocal space maps reveal how different tilted 𝑎 -type domain variants interact with the grain boundary as a function of increasing electric field. Initially, a subset of 𝑎 -type domain variants with their polarization vectors largely orthogonal to the grain boundary were found in greater abundance within 570 nm of the grain boundary, possibly due to X-ray beam-induced local increases in the electrical conductivity, but after the coercive field was exceeded, reconfiguration of the ferroelastic domains was observed. The spatially-varying reciprocal space maps also facilitated evaluation of the strain field across the mapped area, along with 𝑑 33,𝑓 . Significant spatial heterogeneity of strain and 𝑑 33,𝑓 are observed, especially at the coercive field, which was attributed to maintaining deformation compatibility and correlated ferroelastic switching.

36 MATERIALS SCIENCE↗

DOE OSTI2017

A Hybrid Multi-Scale Model of Crystal Plasticity for Handling Stress Concentrations

Microstructural effects become important at regions of stress concentrators such as notches, cracks and contact surfaces. A multiscale model is presented that efficiently captures microstructural details at such critical regions. The approach is based on a multiresolution mesh that includes an explicit microstructure representation at critical regions where stresses are localized. At regions farther away from the stress concentration, a reduced order model that statistically captures the effect of the microstructure is employed. The statistical model is based on a finite element representation of the orientation distribution function (ODF). As an illustrative example, we have applied the multiscaling method to compute the stress intensity factor K I around the crack tip in a wedge-opening load specimen. The approach is verified with an analytical solution within linear elasticity approximation and is then extended to allow modeling of microstructural effects on crack tip plasticity.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions

A discrete De Rham complex enables compatible, structure-preserving discretizations for a broad range of partial differential equations problems. Such discretizations can correctly reproduce the physics of interface problems, provided the grid conforms to the interface. However, large deformations, complex geometries, and evolving interfaces makes generation of such grids difficult. We develop and demonstrate two formally equivalent approaches that, for a given background mesh, dynamically construct an interface-conforming discrete De Rham complex. Both approaches start by dividing cut elements into interface-conforming subelements but differ in how they build the finite element basis on these subelements. The first approach discards the existing non-conforming basis of the parent element and replaces it by a dynamic set of degrees of freedom of the same kind. The second approach defines the interface-conforming degrees of freedom on the subelements as superpositions of the basis functions of the parent element. These approaches generalize the Conformal Decomposition Finite Element Method (CDFEM) and the extended finite element method with algebraic constraints (XFEM-AC), respectively, across the De Rham complex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Bayesian inference of nuclear symmetry energy from measured and imagined neutron skin thickness in Sn 116 , 118 , 120 , 122 , 124 , 130 , 132 , Pb 208 , and Ca 48

The neutron skin thickness Δr np in heavy nuclei has been known as one of the most sensitive terrestrial probes of the nuclear symmetry energy E sym (ρ) around $\frac{2}{3}$ of the saturation density ρ 0 of nuclear matter. Existing neutron skin data mostly from hadronic observables suffer from large uncertainties and their extraction from experiments are often strongly model dependent. While waiting eagerly for the promised model-independent and high-precision neutron skin data for 208 Pb and 48 Ca from the parity-violating electron scattering experiments (PREX-II and CREX at JLab as well as MREX at MESA), within the Bayesian statistical framework using the Skyrme-Hartree-Fock model we infer the posterior probability distribution functions (PDFs) of the slope parameter L of the nuclear symmetry energy at ρ 0 from imagined Δr np ( 208 Pb)=0.15, 0.20, and 0.30 fm with a 1σ error bar of 0.02, 0.04, and 0.06 fm, respectively, as well as Δr np ( 48 Ca)=0.12, 0.15, and 0.25 fm, with different 1σ error bar of 0.01 and 0.02 fm, respectively. The results are compared with the PDFs of L inferred using the same approach from the available Δr np data for 116, 118, 120, 122, 124, 130, 132 Sn from hadronic probes. They are also compared with results from a recent Bayesian analysis of the radius and tidal deformability data of canonical neutron stars from GW170817 and NICER. The neutron skin data for Sn isotopes gives L = 45.5 $^{+ 26.5}_{-21.6}$ MeV surrounding its mean value or L = 53 . 4 $^{+ 18.6}_{ -29.5}$ MeV surrounding its maximum a posteriori value, respectively, with the latter smaller than but consistent with the L = 66 $^{+ 12}_{-20}$ MeV from the neutron star data within their 68% confidence intervals. We found that Δr np = 0.17 –0.18 fm in 208 Pb with an error bar of about 0.02 fm leads to a PDF of L compatible with that from analyzing the Sn data. To provide additionally useful information on L extracted from the Δr np of Sn isotopes, the experimental error bar of Δr np in 208 Pb should be at least smaller than 0.06 fm aimed by some current experiments. In addition, the Δr np ( 48 Ca) needs to be larger than 0.15 fm but smaller than 0.25 fm to be compatible with the Sn and/or neutron star results. To further improve our current knowledge about L and distinguish its PDFs in the examples considered, even higher precisions leading to significantly less than ±20 MeV error bars for L at 68% confidence level are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

A set of 4D pediatric XCAT reference phantoms for multimodality research

Purpose: The authors previously developed an adult population of 4D extended cardiac-torso (XCAT) phantoms for multimodality imaging research. In this work, the authors develop a reference set of 4D pediatric XCAT phantoms consisting of male and female anatomies at ages of newborn, 1, 5, 10, and 15 years. These models will serve as the foundation from which the authors will create a vast population of pediatric phantoms for optimizing pediatric CT imaging protocols. Methods: Each phantom was based on a unique set of CT data from a normal patient obtained from the Duke University database. The datasets were selected to best match the reference values for height and weight for the different ages and genders according to ICRP Publication 89. The major organs and structures were segmented from the CT data and used to create an initial pediatric model defined using nonuniform rational B-spline surfaces. The CT data covered the entire torso and part of the head. To complete the body, the authors manually added on the top of the head and the arms and legs using scaled versions of the XCAT adult models or additional models created from cadaver data. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from a template XCAT phantom (male or female 50th percentile adult) to the target pediatric model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. The masses of the organs in each phantom were matched to the reference values given in ICRP Publication 89. The new reference models were checked for anatomical accuracy via visual inspection. Results: The authors created a set of ten pediatric reference phantoms that have the same level of detail and functionality as the original XCAT phantom adults. Each consists of thousands of anatomical structures and includes parameterized models for the cardiac and respiratory motions. Based on patient data, the phantoms capture the anatomic variations of childhood, such as the development of bone in the skull, pelvis, and long bones, and the growth of the vertebrae and organs. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. Conclusions: The development of patient-derived pediatric computational phantoms is useful in providing variable anatomies for simulation. Future work will expand this ten-phantom base to a host of pediatric phantoms representative of the public at large. This can provide a means to evaluate and improve pediatric imaging devices and to optimize CT protocols in terms of image quality and radiation dose.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

Mechanical Switching of Nanoscale Multiferroic Phase Boundaries

Tuning the lattice degree of freedom in nanoscale functional crystals is critical to exploit the emerging functionalities such as piezoelectricity, shape-memory effect, or piezomagnetism, which are attributed to the intrinsic lattice-polar or lattice-spin coupling. In this study, it is reported that a mechanical probe can be a dynamic tool to switch the ferroic orders at the nanoscale multiferroic phase boundaries in BiFeO 3 with a phase mixture, where the material can be reversibly transformed between the “soft” tetragonal-like and the “hard” rhombohedral-like structures. The microscopic origin of the nonvolatile mechanical switching of the multiferroic phase boundaries, coupled with a reversible 180° rotation of the in-plane ferroelectric polarization, is the nanoscale pressure-induced elastic deformation and reconstruction of the spontaneous strain gradient across the multiferroic phase boundaries. The reversible control of the room-temperature multiple ferroic orders using a pure mechanical stimulus may bring us a new pathway to achieve the potential energy conversion and sensing applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Structural Design Variations in Beetle Elytra

Abstract Beetles typically use their protective wing coverings or elytra to shield their membranous hindwings from the environment. Elytra in some terrestrial species have evolved a greater protective role capable of shielding the organism from powerful antagonistic predators. The structure–function relationships of these biological composites identify how architectural and chemical variations of the cuticle are tuned to create light‐weight, impact resistant composites. Specifically, the elytral structures of a tree dwelling beetle capable of flight, Trypoxylus dichotomus , and a terrestrial beetle incapable of flight, Phloeodes diabolicus , are compared to understand how their varied environmental needs forged the elytra to facilitate fight or resist fatal predator strikes. Mechanical and microstructural analysis reveals P . diabolicus has a harder, stiffer elytra that incorporate through‐thickness fibers to resist greater mechanical stresses imposed by bending and puncture. Conversely, the elytra of T . dichotomus have a compliant structure with large voids that facilitates localized deformation. Variations in flexural strength and puncture resistance remain attributed to P . diabolicus possessing a thicker cuticle with a greater degree of cross‐linking and an increased amount of endocuticular layers. These findings may provide useful insight into the design and manufacturing of composite materials for use in light‐weight or energy‐absorbing applications.

Rivera, Jesus↗

DOE OSTI2014

Microstructure characterization of laser welded Ti-6Al-4V fusion zones

The as-welded microstructure of laser-welded Ti-6Al-4V is characterized as a function of CO2 key-hole mode laser welding speed. Martensitic α′ is the predominant phase, with some α and retained β. Phase transformation is affected by the cooling rate through laser welding speed. A higher welding speed of 1.6 to 2.0 m/min produced more martensite α′ and less retained β in the welds. 1.4 m/min welding speed produced small amounts of α, besides the martensite α′. A trace of δ titanium hydride phase seems to have formed in the weld fusion zone. Moiré fringes are a common feature in the TEM microstructure, due to abundance of multi-phase interfaces. Tensile twins and clusters of dislocations indicate that plastic deformation has happened in the as-welded microstructure, indicating the local stress levels to be approaching the yield stress on-cooling during laser welding.

36 MATERIALS SCIENCE↗

DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

First-principles study of crystallographic slip modes in ω-Zr

We use first-principles density functional theory to study the preferred modes of slip in the high-pressure ω phase of Zr. The generalized stacking fault energy surfaces associated with shearing on nine distinct crystallographic slip modes in the hexagonal ω-Zr crystal are calculated, from which characteristics such as ideal shear stress, the dislocation Burgers vector, and possible accompanying atomic shuffles, are extracted. Comparison of energy barriers and ideal shear stresses suggests that the favorable modes are prismatic < c >, prismatic-II <101¯0> and pyramidal-II < c+a >, which are distinct from the ground state hexagonal close packed α phase of Zr. Operation of these three modes can accommodate any deformation state. The relative preferences among the identified slip modes are examined using a mean-field crystal plasticity model and comparing the calculated deformation texture with the measurement. In conclusion, knowledge of the basic crystallographic modes of slip is critical to understanding and analyzing the plastic deformation behavior of ω-Zr or mixed α-ω phase-Zr.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Above-bandgap optical properties of biaxially strained GeSn alloys grown by molecular beam epitaxy

The complex dielectric function of biaxially strained Ge{sub 1−x}Sn{sub x} (0 ≤ x ≤ 0.17) alloys grown on Ge (100) has been determined by spectroscopic ellipsometry from 1.2 to 4.7 eV. The effect of substitutional Sn incorporation and the epitaxial strain on the energy transitions E{sub 1}, E{sub 1} + Δ{sub 1}, E{sub 0}′, and E{sub 2} of GeSn alloys is investigated. Our results indicate that the strained GeSn alloys show Ge-like electronic bandstructure with all the transitions shifted downward due to the alloying of Sn. The strain dependence of E{sub 1} and E{sub 1} + Δ{sub 1} transitions is explained using the deformation potential theory, and values of −5.4 ± 0.4 eV and 3.8 ± 0.5 eV are obtained for the hydrostatic and shear deformation potentials, respectively.

36 MATERIALS SCIENCE↗