Search NASA⌕ Search

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.

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

At least 1,081 records · Page 60

Collaborative Research: The Role of Grain Boundary Structure and Chemistry in Materials Failure

Throughout the one year duration of this award significant progress was made in understanding the effect of grain boundaries in polycrystalline materials. The starting point of this work is that each grain boundary can be characterized by an interface energy and respective interface parameter, which depend on the boundary chemical composition and structure. As this is a new point of view for treating grain boundaries, an interdisciplinary approach was followed to prove the existence of such interface terms. Throughout the project, hence, theory, experiments and simulations were used. Molecular dynamics and dislocation density-function dynamics were performed to understand dislocation-grain boundary interactions in pure Al crystalline materials, as potentials for Al are well established. The formation of dislocation pile-ups was observed and the strain profiles were in agreement with interfacial gradient plasticity, which is the theoretical framework that allows the introduction of mechanically induced interface energies. After verifying numerically the existence of mechanically induced interface energies, nanoindentation experiments were performed near grain boundaries of Fe-Si crystals that contained three grain boundaries. By indenting in close proximity to the boundary these experiments allowed for measuring the critical stress at which deformation was transmitted across the boundary to the adjacent grains, either through dislocation transmission or absorption and subsequent emission. This critical stress, also termed grain boundary yield stress, was in accordance with the value predicted by interfacial gradient plasticity, further verifying that mechanically induced interface energies can help define the material response. In the Fe-Si crystals the Si content was kept constant at 3wt%, but three specimen were considered each of which was characterized by a different misorientation at the grain boundary. By relating the experiments with the analytical gradient plasticity expressions it was possible to obtain first estimates for the interface parameter as ~200N/m. This value was not proportional to the misorientation angle, but depended also on the grain boundary plane. The continuation of this project is performing more detailed experiments on additional Fe-Si and Fe-C systems (both nanoindentation and in-situ transmission electron microscopy) to further understand the physical origin of such interface parameters. Furthermore, molecular dynamics nanoindentation simulations are also performed on similar systems to better understand the dislocation grain boundary interactions.

36 MATERIALS SCIENCE↗

Tuning Ultrafast Photoinduced Strain in Ferroelectric-Based Devices

Ferroelectric materials exhibit coupled degrees of freedom and possess a switchable electric polarization coupled to strain, making them good piezoelectrics and enabling numerous devices including nonvolatile memories, actuators, and sensors. Moreover, novel photovoltaic effects are encountered through the interplay of electric polarization with the material optical properties. Consequently, light-induced deformation in ferroelectrics, or photostriction, combining photovoltaic and converse piezoelectric effects, is under investigation in the quest for multifunctional materials. Using time-resolved X-ray diffraction, the first control of ultrafast photoinduced strain is demonstrated through in situ tuning of the polarization state in ferroelectric-based devices. Both the magnitude and the sign of the photoinduced strain strongly depend on the transient photoinduced change of the internal electric field in the ferroelectric layer, and can be actively engineered to achieve two distinct remanent photostrictive responses. These results provide fundamental insight into light–matter interaction in ferroelectrics and exciting new avenues for materials functionality engineering.

36 MATERIALS SCIENCE↗

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↗

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↗

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↗

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↗

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↗

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↗

$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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Ab Initio Flexible Force Field for Metal–Organic Frameworks Using Dummy Model Coordination Bonds

We present force fields developed from periodic density functional theory (DFT) calculations that can be used in classical molecular simulations to model M–MOF-74 (M = Co, Fe, Mg, Mn, Ni, Zn) and its extended linker analogs. Our force fields are based on cationic dummy models (CDMs). These dummy models simplify the methodology required to tune the parameters and improve the accuracy of the force fields. We used our force fields to compare mechanical properties across the M–MOF-74 series and determine that increasing the size of the linker decreases the framework rigidity. In addition, we applied our force fields to an extended linker analog of Mg–MOF-74 and characterized the free energy of a previously reported deformation pattern in which the one-dimensional hexagonal channels of the framework become irregular. The free energy profiles confirm that the deformation is adsorbate induced and impossible to access solely by a pressure stimulus. On the basis of our results, we conclude that the force fields presented here and others that may be developed using our methodology are transferable across metal–organic framework series that share a metal center topology. Lastly, we believe that these force fields have the potential to be adapted for the study of complex problems in MOF chemistry, including defects and crystal growth, that have thus far been beyond the scope of classical molecular simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-objective parametrization of interatomic potentials for large deformation pathways and fracture of two-dimensional materials

Abstract This investigation presents a generally applicable framework for parameterizing interatomic potentials to accurately capture large deformation pathways. It incorporates a multi-objective genetic algorithm, training and screening property sets, and correlation and principal component analyses. The framework enables iterative definition of properties in the training and screening sets, guided by correlation relationships between properties, aiming to achieve optimal parametrizations for properties of interest. Specifically, the performance of increasingly complex potentials, Buckingham, Stillinger-Weber, Tersoff, and modified reactive empirical bond-order potentials are compared. Using MoSe 2 as a case study, we demonstrate good reproducibility of training/screening properties and superior transferability. For MoSe 2 , the best performance is achieved using the Tersoff potential, which is ascribed to its apparent higher flexibility embedded in its functional form. These results should facilitate the selection and parametrization of interatomic potentials for exploring mechanical and phononic properties of a large library of two-dimensional and bulk materials.

Zhang, Xu↗