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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An interplay between a hydrogen atmosphere and dislocation characteristics in BCC Fe from time-averaged molecular dynamics

The interplay between hydrogen and dislocations (e.g., core and elastic energies, and dislocation–dislocation interactions) has implications on hydrogen embrittlement but is poorly understood. Continuum models of hydrogen enhanced local plasticity have not considered the effect of hydrogen on dislocation core energies. Energy minimization atomistic simulations can only resolve dislocation core energies in hydrogen-free systems because hydrogen motion is omitted so hydrogen atmosphere formation can’t occur. Additionally, previous studies focused more on face-centered-cubic than body-centered-cubic metals. Discrete dislocation dynamics studies of hydrogen–dislocation interactions assume isotropic elasticity, but the validity of this assumption isn’t understood. Here, we perform time-averaged molecular dynamics simulations to study the effect of hydrogen on dislocation energies in body-centered-cubic iron for several dislocation character angles. We see atmosphere formation and highly converged dislocation energies. We find that hydrogen reduces dislocation core energies but can increase or decrease elastic energies of isolated dislocations and dislocation–dislocation interaction energies depending on character angle. We also find that isotropic elasticity can be well fitted to dislocation energies obtained from simulations if the isotropic elastic constants are not constrained to their anisotropic counterparts. These results are relevant to ongoing efforts in understanding hydrogen embrittlement and provide a foundation for future work in this field.

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

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Giant microwave absorption in the vortex lattice in s -wave superconductors

In this article we study microwave absorption in superconductors in the presence of a vortex lattice. Here, we show that in addition to the conventional absorption mechanism associated with the vortex core motion, there is another mechanism of microwave absorption, which is caused by the time dependence of the quasiparticle density of states outside the vortex cores. This mechanism exists even in the absence of vortex motion. It provides the dominant contribution to microwave absorption at sufficiently small magnetic fields. At low frequencies, the dissipative part of the microwave conductivity σ⁡(ω) is proportional to the inelastic relaxation time τ in , which is typically much larger than the elastic relaxation time τ el . At high frequencies σ⁡(ω) is proportional to the quasiparticle diffusion time across the intervortex distance τ D , which is still larger than τ el .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Field emission energy distributions of hydrogen and deuterium on the /100/ and /110/ planes of tungsten.

Total energy distributions of field emitted electrons from the tungsten (110) and (100) planes as a function of coverage by hydrogen and deuterium have been recorded utilizing a spherical deflection energy analyzer. The elastic tunneling resonance spectrum gives a plot of the 'local density of states' in the adsorbate. The inelastic tunneling spectrum reveals those discrete excitation energies available in the adsorbate-substrate complex. These spectroscopic data have been used to infer the chemical nature of the binding states which have been observed in the flash desorption spectrum of hydrogen from tungsten.

Plummer, E. W.↗
NASA NTRS1983

Edge effects in composites by moire interferometry

The very high sensitivity of moire interferometry has permitted the present edge effect experiments to be conducted at a low average stress and strain level, assuring linear and elastic behavior in the composite material samples tested. Sensitivity corresponding to 2450 line/mm moire was achieved with a 0.408 micron/fringe. Simultaneous observations of the specimen face and edge displacement fields showed good fringe definition despite the 1-mm thickness of the specimens and the high gradients, and it is noted that the use of a carrier pattern and optical filtering was effective in even these conditions. Edge effects and dramatic displacement gradients were confirmed in angle-ply composite laminates.

Czarnek, R.↗
NASA NTRS1987

A quadrilateral shell element using a mixed formulation

A simple quadrilateral shell element consisting of five nodes, four corner nodes and a central node, is developed for linear elastic analysis of thin as well as moderately thick shells. Based on a modified Hellinger-Reissner principle, finite element equations are derived from the assumed displacement and strain fields. By carefully choosing appropriate strain terms, all kinematic deformation modes are suppressed. Although the present element is similar to a displacement-based degenerated shell, no locking is experienced when it is applied to thin shell problems. Five examples are given to illustrate the analysis capability of the shell element. Numerical results indicate that the element shows fast mesh convergence and gives excellent stress predictions.

Saleeb, A. F.↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2020

Wellbore injectivity response to step-rate CO 2 injection: Coupled thermo-poro-elastic analysis in a vertically heterogeneous formation

Safe and permanent carbon geological storages require an elaborate analysis of geomechanical stability. Subsurface injection of CO 2 changes the local temperature and pore pressure and, further, alters the stress due to thermo-poro-elastic responses. A permanent CO 2 storage testing was conducted in the water leg of Cranfield reservoir in Mississippi, USA, where hosted CO 2 enhanced-oil recovery activities. During CO 2 injection, the injection rate was ramped up twice but the bottom-hole pressure did not increase with the imposed injection rates as expected. This unexpected field observation suggests the possibility of an open-mode fracture development at the injector. However, the injector response and potential fracture development have not been rigorously interpreted with detailed near-wellbore temperature and pore pressure change upon the CO 2 injection. In this study, we performed history matching of CO 2 injection using coupled thermo-poro-elastic reservoir simulation to examine the possibility of fracturing. We built a reservoir model including vertical heterogeneity in both petrophysical and geomechanical properties estimated from well-logging analysis and laboratory experiments. The simulation results show that CO 2 injection changes stresses and support the hypothesis of development of an open-mode fracture at the injector during the Cranfield test. The near-injector region exhibits a large temperature reduction up to 55 °C with ensuing effective horizontal stress reduction up to 9.1 MPa. However, a caprock integrity issue is unlikely because of the horizontal stress contrast within layers and low hydraulic communication with the injection zone. This work indicates that the injectant temperature should be considered in the design of high-rate CO 2 injector.

58 GEOSCIENCES↗
DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗
DOE OSTI2022

Two–Photon Printing of Shape–Memory Microstructures and Metasurfaces via Radical–Mediated Thiol–Vinyl Hydrothiolation

Shape-memory resists capable of high-resolution curing into arbitrarily designed structures are increasingly demanded for soft robotics, optical sensors, microscale manufacturing, and biomedicine. Amorphous, shape-memory thiol-vinyl networks were printed using two-photon polymerization (2PP) curing of a simple resin formulated with commercially available reagents. The ability to print high-resolution feature sizes down to 200 nm is attributed to the use of radicalmediated, thiol-vinyl step-growth polymerization that quickly crosslinks the resin, limiting diffusive transport. The thermomechanical behavior of the 2PP-cured material analyzed in compression, tension and three-point bending is similar to the behavior of the UV-polymerized samples. To demonstrate the ability to design, field, and test 4D responsive microstructures, an array of nine springs with coil diameters of 330 μm was printed. Following compressive shapefixing, printed arrays can release 11 μJ of stored elastic strain energy when reheated. Further, a new concept of dichroic-memory of a metamaterial device is demonstrated by printing a twisted woodpile structure with circular dichroism as characterized by Mueller Matrix ellipsometry. Furthermore, our results demonstrate how combining high-resolution 2PP curing with stimuli-responsive molecular architectures can further the engineering of responsive microstructures and metamaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Morphology and Dynamics of Bottlebrush Polymers

This thesis makes significant advances to the understanding of bottlebrush polymers. While bottlebrushes have received much attention due to the recent discovery of their unprecedented properties, including supersoftness, ultra-low viscosity, and hyperelasticity, this thesis is the first fundamental investigation at the molecular level that comprises structure and dynamics. Neutron scattering experiments, detailed within, reveal spherical or cylindrical shapes, instead of a random coil conformation. Another highlight is the analysis of the fast dynamics at the sub nm-length scale. The combination of three neutron spectrometers and the development of a new analysis technique enabled the calculation of the mean-square displacement over seven orders of magnitude in time scale. This unprecedented result can be applied to a broad class of samples, including polymers and other materials. The thesis is accessible to scientists from other fields, provides the reader with easily understandable guidelines for applying this analysis to other materials, and has the potential to make a significant impact on the analysis of neutron scattering data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2023

Determining Stress in Metallic Conducting Layers of Microelectronics Devices Using High Resolution Electron Backscatter Diffraction and Finite Element Analysis

Delayed failure due to stress voiding is a concern with some aging microelectronics, as these voids can grow large enough to cause an open circuit. Local measurements of stress in the metallic layers are crucial to understanding and predicting this failure, but such measurements are complicated by the fact that exposing the aluminum conducting lines will relieve most of their stress. In this study, we instead mechanically thin the device substrate and measure distortions on the thinned surface using high resolution electron backscatter diffraction (HREBSD). These measurements are then related to the stresses in the metallic layers through elastic simulations. This study found that in legacy components that had no obvious voids, the stresses were comparable to the theoretical stresses at the time of manufacture (≈300 MPa). Distortion fields in the substrate were also determined around known voids, which may be directly compared to stress voiding models. In conclusion, the technique presented here for stress determination, HREBSD coupled with finite element analysis to infer subsurface stresses, is a valuable tool for assessing failure in layered microelectronics devices.

HREBSD↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
NASA NTRS1982

On ultrasonic factors and fracture toughness

Recent experimental and theoretical studies on ultrasonics have shown that the scattering of elastic waves by material defects yields data which characterize crack properties, such as size and orientation, and also the mechanical properties of the given material. In the present study, elastodynamic fields due to the presence of a pair of inhomogeneities in a material of plate geometry are investigated by the method of equivalent inclusions. The stress amplitude change of the plates during the passage of plane time-harmonic waves is found, and the relation between fracture toughness and ultrasonic factors is determined. The approach used does not assume the existence of a sharp crack in the material.

Fu, L. S.↗
NASA NTRS1986

Design of anisotropic plates for improved damage tolerance

An analytical study is presented showing the effects of the notch tip geometry on the location and direction of crack growth from an existing notch in a unidirectional fibrous composite modeled as a homogeneous, anisotropic, elastic material. Anisotropic elasticity and the normal stress ratio theory are used to study crack growth from elliptical notches in unidirectional composites. Sharp cracks, circular holes, and ellipses are studied under far-field tension and shear loading. The capabilities of a previously developed design code was upgraded to handle more generalized plate geometries and laminates under a more generalized loading and boundary conditions. Discussion of the developments of the design code is presented.

Guerdal, Zafer↗