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

Neutron Scattering (NS) Spectroscopy

A wide range of advanced experimental methods has been used in catalysis science to understand chemical transformations at molecular level. Among these, neutron scattering not only gives catalytic information that is highly complementary to other microscopic scattering techniques such as electrons (microscopy and diffraction) and photons from visible light to synchrotron X-rays but also often provides unique insights into catalysis. This chapter aims to provide a general overview of neutron scattering and its applications for heterogeneous catalysis, emphasizing chemistry at the gas/solid interface. It starts with an introduction to the theory of neutron scattering, with pros and cons for catalyst research, followed by a description of three main neutron scattering techniques: neutron diffraction (ND), inelastic neutron spectroscopy (INS), and quasi-elastic neutron scattering (QENS). It includes recent applications in each technique and then a short introduction of other neutron techniques that are less frequently used in the catalysis field. It ends with a summary and a future outlook.

Braatz, Jisue↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗
DOE OSTI2026

Design and commissioning of the Correlated Gamma-Neutron Array for sCattering (CoGNAC)

The Correlated Gamma-Neutron Array for sCattering (CoGNAC) is a dual $n$-$γ$ detector array fielded at the Weapons Neutron Research facility at the Los Alamos Neutron Science Center. Here, the detector array utilizes elpasolite detectors and liquid scintillators to perform measurements of neutron elastic and inelastic scattering cross sections and emitted particle angular distributions. This manuscript details the design and experimental setup of CoGNAC as well as some of the preliminary studies performed to select the elpasolite detectors deployed in the array. Examples of the analysis methods available to CoGNAC are given in the context of past CoGNAC measurements of neutron scattering reactions on 28 Si and 56 Fe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
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↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗