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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Strength and rigidity of the elastic lunar lithosphere and implications for present-day mantle convection in the moon

Lunar gravity data are used to estimate the thickness of the elastically defined lunar lithosphere, and the load-bearing capability of lithosphere with respect to lunar mascons is analyzed by representing the elastic lithosphere as a plane elastic plate subject to spatially concentrated loads from above and buoyant support from below. Gravity perturbations are computed and compared with observations, the spectrum of lunar gravity perturbations is derived, and effects of curvature are considered. The magnitude of the elastic stresses associated with a lithospheric flexure model maintaining the irregularities of the lunar gravity field is shown to be quite similar to that of the lithospheric stresses found on earth. It is concluded that associated differences in gravity, heat flow, and lithospheric thickness all combine to make the lunar lithosphere much more rigid and stronger than that of earth. Implications for present-day lunar convection are discussed.

Kuckes, A. F.↗
NASA NTRS1987

A substructure technique for dynamics of flexible mechanical systems with contact-impact

A substructure synthesis method is proposed to account for contact-impact effects in flexible components of mechanical systems. Components that may come into contact is divided into substructures, on each of which local deformation modes are defined to described deformation fields of components. Constraint modes and fixed interface normal modes are used to account for elastic deformation within each substructure. A constraint addition-deletion technique is used to account for contact between impacting bodies. Lagrange multipliers associated with the constraints, which represent constraint forces, are used to determine separation of contacting nodes. Use of the method is illustrated on problems of longitudinal and transverse impact of bodies.

Wu, S.-C.↗
NASA NTRS2015

Test and Analysis Correlation for a Y-Joint Specimen for a Composite Cryotank

The Composite Cryotank Technology Demonstration (CCTD) project under NASA's Game Changing Development Program (GCDP) developed space technologies using advanced composite materials. Under CCTD, NASA funded the Boeing Company to design and test a number of element-level joint specimens as a precursor to a 2.4-m diameter composite cryotank. Preliminary analyses indicated that the y-joint in the cryotank had low margins of safety; hence the y-joint was considered to be a critical design region. The y-joint design includes a softening strip wedge to reduce localized shear stresses at the skirt/dome interface. In this paper, NASA-developed analytical models will be correlated with the experimental results of a series of positive-peel y-joint specimens from Boeing tests. Initial analytical models over-predicted the experimental strain gage readings in the far-field region by approximately 10%. The over-prediction was attributed to uncertainty in the elastic properties of the laminate and a mismatch between the thermal expansion of the strain gages and the laminate. The elastic properties of the analytical model were adjusted to account for the strain gage differences. The experimental strain gages also indicated a large non-linear effect in the softening strip region that was not predicted by the analytical model. This non-linear effect was attributed to delamination initiating in the softening strip region at below 20% of the failure load for the specimen. Because the specimen was contained in a thermally insulated box during cryogenic testing to failure, delamination initiation and progression was not visualized during the test. Several possible failure initiation locations were investigated, and a most likely failure scenario was determined that correlated well with the experimental data. The most likely failure scenario corresponded to damage initiating in the softening strip and delamination extending to the grips at final failure.

Mason, Brian H.↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cross Sections for Coherent Elastic and Inelastic Neutrino-Nucleus Scattering

The prospects of extracting new physics signals in coherent elastic neutrino–nucleus scattering (CEνNS) processes are limited by the precision with which the underlying nuclear structure physics, embedded in the weak nuclear form factor, is known. We present calculations of charge and weak nuclear form factors and CEνNS cross sections on 12C, 16O, 40Ar, 56Fe and 208Pb nuclei. We obtain the proton and neutron densities, and charge and weak form factors by solving Hartree–Fock (HF) equations with a Skyrme (SkE2) nuclear potential. We validate our approach by comparing 208Pb and 40Ar charge form factor predictions with available elastic electron scattering data. Since CEνNS experiments at stopped-pion sources are also well suited to measure inelastic charged–current and neutral–current neutrino–nucleus cross sections, we also present calculations for these processes, incorporating a continuum Random Phase Approximation (CRPA) description on top of the HF–SkE2 picture of the nucleus. Providing both coherent as well as inelastic cross sections in a consistent framework, we aim at obtaining a reliable and detailed comparison of the strength of these processes in the energy region below 100 MeV. Furthermore, we attempt to gauge the level of theoretical uncertainty pertaining to the description of the 40Ar form factor and CEνNS cross sections by comparing relative differences between recent microscopic nuclear theory and widely-used phenomenological form factor predictions. Future precision measurements of CEνNS will potentially help in constraining these nuclear structure details that will in turn improve prospects of extracting new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

In‐Operando Spatiotemporal Imaging of Coupled Film‐Substrate Elastodynamics During an Insulator‐to‐Metal Transition

Abstract The drive toward non‐von Neumann device architectures has led to an intense focus on insulator‐to‐metal (IMT) and the converse metal‐to‐insulator (MIT) transitions. Studies of electric field‐driven IMT in the prototypical VO 2 thin‐film channel devices are largely focused on the electrical and elastic responses of the films, but the response of the corresponding TiO 2 substrate is often overlooked, since it is nominally expected to be electrically passive and elastically rigid. Here, in‐operando spatiotemporal imaging of the coupled elastodynamics using X‐ray diffraction microscopy of a VO 2 film channel device on TiO 2 substrate reveals two new surprises. First, the film channel bulges during the IMT, the opposite of the expected shrinking in the film undergoing IMT. Second, a microns thick proximal layer in the substrate also coherently bulges accompanying the IMT in the film, which is completely unexpected. Phase‐field simulations of coupled IMT, oxygen vacancy electronic dynamics, and electronic carrier diffusion incorporating thermal and strain effects suggest that the observed elastodynamics can be explained by the known naturally occurring oxygen vacancies that rapidly ionize (and deionize) in concert with the IMT (MIT). Fast electrical‐triggering of the IMT via ionizing defects and an active “IMT‐like” substrate layer are critical aspects to consider in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗
DOE OSTI2022

Elastic Solutions to 2D Plane Strain Problems: Nonlinear Contact and Settlement Analysis for Shallow Foundations

The classical Neumann boundary value problem of an isotropic, homogeneous elastic half-plane under plane strain conditions is readdressed as the limiting case of the fully three-dimensional problem. Analytical solutions of the stress and strain tensors are obtained by taking the limit from known three-dimensional solutions. It is shown that the displacement fields for the plane strain problem are not well defined. A small number of simple expressions are developed, which provide a general solution for linearly-varying traction over arbitrary regions on the boundary. A simple, efficient, and rapidly convergent algorithm is developed which uses these solutions as analytic elements and provides a solution approach to the general boundary value problem. The method is verified against known solutions for Hertzian contact between parallel cylinders. Two numerical examples are presented for the analysis of shallow foundation systems. In the first, the boundary conditions are informed by analytical elastoplastic calculations and a strain influence analysis is performed and compared with the Schmertmann method. Subsequently, empirical laboratory contact traction distributions measured by Bauer et al., in both the normal and tangential directions are employed as boundary conditions for an analysis of the underlying stress field.

42 ENGINEERING↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2025

Magnetic properties of a staggered S = 1 chain with an alternating single-ion anisotropy direction

Materials composed of spin-1 antiferromagnetic (AFM) chains are known to adopt complex ground states that are sensitive to the single-ion-anisotropy (SIA) energy ( D ), and intrachain ( J 0 ) and interchain ( J 1 , 2 ′ ) exchange energy scales. While theoretical and experimental studies have extended this model to include various other energy scales, the effect of the lack of a common SIA axis is not well explored. Here we investigate the magnetic properties of Ni ( pyrimidine ) ( H 2 O ) 2 ( NO 3 ) 2 , a chain compound where the tilting of Ni octahedra leads to a twofold alternation of the easy-axis directions along the chain. Muon-spin relaxation measurements indicate a transition to long-range order at T N = 2.3 K and the magnetic structure is initially determined to be antiferromagnetic and collinear using elastic neutron diffraction experiments. Inelastic neutron scattering measurements were used to find J 0 = 5.107 ( 7 ) K , D = 2.79 ( 1 ) K , J 1 ′ = 0.00 ( 5 ) K , J 2 ′ = 0.18 ( 3 ) K , and a rhombic anisotropy energy E = 0.19 ( 9 ) K . Mean-field modeling reveals that the ground state structure hosts spin canting of ϕ ≈ 6 . 5 ∘ , which is not detectable above the noise floor of the elastic neutron diffraction data. Monte Carlo simulation of the powder-averaged magnetization, M ( H ) , is then used to confirm these Hamiltonian parameters, while single-crystal M ( H ) simulations provide insight into features observed in the data. Published by the American Physical Society 2025

Vaidya, S. (ORCID:000000033975497X)↗
NASA NTRS1980

Noise transmission and attenuation by stiffened panels

An analytical study of noise transmission into semi-cylindrical and rectangular acoustic enclosures due to turbulent boundary layer pressure and propeller noise (prop-fan) is presented. The structural noise transmission models include a single panel, discretely stiffened elastic panel and stiffened viscoelastic sandwich panel. Response characteristics of the stiffened panels are evaluated using a transfer matrix procedure. The interior noise field is determined by a Galerkin-like method. The effect on interior noise due to aerodynamic surface flow, cavity back-up pressure, pressurization, mass, stiffness, and damping addition to the structure is investigated. It is shown that stiffened viscoelastic sandwich panels, while providing the same stiffening benefits as an equivalent elastic panel, could significantly reduce vibration levels and subsequently give similar benefits for interior noise control.

Vaicaitis, R.↗
DOE OSTI2021

Complex Oxides under Simulated Electric Field: Determinants of Defect Polarization in AB O 3 Perovskites

Abstract Polarization of ionic and electronic defects in response to high electric fields plays an essential role in determining properties of materials in applications such as memristive devices. However, isolating the polarization response of individual defects has been challenging for both models and measurements. Here the authors quantify the nonlinear dielectric response of neutral oxygen vacancies, comprised of strongly localized electrons at an oxygen vacancy site, in perovskite oxides of the form AB O 3 . Their approach implements a computationally efficient local Hubbard U correction in density functional theory simulations. These calculations indicate that the electric dipole moment of this defect is correlated positively with the lattice volume, which they varied by elastic strain and by A‐site cation species. In addition, the dipole of the neutral oxygen vacancy under electric field increases with increasing reducibility of the B‐site cation. The predicted relationship among point defect polarization, mechanical strain, and transition metal chemistry provides insights for the properties of memristive materials and devices under high electric fields.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Finite domain solution of a hydraulic fracture in a permeable rock

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media.

Domain-based method↗
DOE OSTI2022

Fathoming the anisotropic magnetoelasticity and magnetocaloric effect in GdNi

Intermetallic GdNi adopts a CrB type of crystal structure (space group Cmcm), and it orders ferromagnetically via a second-order phase transition at 70 K, exhibiting unusually strong spontaneous striction along the three independent crystallographic axes in the ferromagnetically ordered state. Here we introduce a microscopic model to describe anisotropic changes of lattice parameters and elastic contribution to magnetocaloric effect of GdNi. In the model, results of density functional theory (DFT) calculations are used as inputs into a Hamiltonian that includes elastic energy of an anisotropic crystal lattice, exchange interactions, and Zeeman effect. The magnetic and elastic Hamiltonians are coupled through an anisotropic Bean-Rodbell model of magnetoelastic interactions. This coupling gives rise to anisotropic changes in the lattice parameters observed experimentally, and the model reveals good to reasonable agreements between the current theoretical results and earlier experimental data, thus validating the model within the limits of assumptions made. We also show that DFT calculations with 4f electrons of Gd treated as core electrons lead to a more adequate estimate of elastic constants of GdNi in comparison with the LDA + U method where 4f electrons are treated as valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗