Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational inelasticity”

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 OSTI2024

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗
DOE OSTI2021

Uranium Carbide and Uranium-Metal Thermal Scattering Law Evaluations and Cross Sections [Slides]

New Uranium Carbide and Uranium Metal evaluations have been developed using the FLASSH platform. Ab initio LD/MD models converged electronic, magnetic, and dynamical properties and demonstrate good agreement with experimental data. Thermal scattering cross sections generate inelastic and elastic thermal neutron scattering libraries in FLASSH . New TSL evaluations occurred in the production of ENDF File 7 data for Uranium Carbide and Uranium Metal submitted to NNDC for ENDF/B-VIII.1 release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1985

Unified constitutive material models for nonlinear finite-element structural analysis

Unified constitutive material models were developed for structural analyses of aircraft gas turbine engine components with particular application to isotropic materials used for high-pressure stage turbine blades and vanes. Forms or combinations of models independently proposed by Bodner and Walker were considered. These theories combine time-dependent and time-independent aspects of inelasticity into a continuous spectrum of behavior. This is in sharp contrast to previous classical approaches that partition inelastic strain into uncoupled plastic and creep components. Predicted stress-strain responses from these models were evaluated against monotonic and cyclic test results for uniaxial specimens of two cast nickel-base alloys, B1900+Hf and Rene' 80. Previously obtained tension-torsion test results for Hastelloy X alloy were used to evaluate multiaxial stress-strain cycle predictions. The unified models, as well as appropriate algorithms for integrating the constitutive equations, were implemented in finite-element computer codes.

Kaufman, A.↗
NASA NTRS1985

Unified constitutive material models for nonlinear finite-element structural analysis

Unified constitutive material models were developed for structural analyses of aircraft gas turbine engine components with particular application to isotropic materials used for high-pressure stage turbine blades and vanes. Forms or combinations of models independently proposed by Bodner and Walker were considered. These theories combine time-dependent and time-independent aspects of inelasticity into a continuous spectrum of behavior. This is in sharp contrast to previous classical approaches that partition inelastic strain into uncoupled plastic and creep components. Predicted stress-strain responses from these models were evaluated against monotonic and cyclic test results for uniaxial specimens of two cast nickel-base alloys, B1900+Hf and Rene 80. Previously obtained tension-torsion test results for Hastelloy X alloy were used to evaluate multiaxial stress-strain cycle predictions. The unified models, as well as appropriate algorithms for integrating the constitutive equations, were implemented in finite-element computer codes.

Kaufman, A.↗
NASA NTRS1988

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗
DOE OSTI2023

Coupled Hydromechanical Modeling and Assessment of Induced Seismicity at FWU: Utilizing Time-Lapse VSP and Microseismic Data

The objective of this work is to utilize integrated geomechanics, field vertical seismic profile (VSP) and microseismic data to characterize the complex subsurface stress conditions at the Farnsworth Unit (FWU). The model is based on a five-spot sector model extracted from a primary geomechanical model. The five-spot well injection pattern is characterized by extensive reservoir characterization data, such well logs, extracted cores and borehole geophone data, to facilitate the detailed examination of stress changes and microseismic event occurrences. The study utilizes field vertical seismic volumes acquired from the injection well 13-10A. The seismic volumes successfully provided snapshots of the behavior of the reservoir at distinct times. The use of VSP and microseismic data provided direct and indirect estimates of the dynamic stress changes occurring in the overburden, reservoir and underburden rock formations. In order to illuminate the stress regions and identify rocks that have undergone inelastic failure, microseismic event occurrences were utilized. Microseismic activity has been detected at the FWU; further study of its locations, timing, and magnitude was needed to deduce the nature of the changing stress state. The results of the study revealed that microseismic events were successfully modeled within the Morrow B formation. Moment magnitudes of seismic events were within the same magnitudes for events in the reservoir, suggesting the suitability of the model. The results of the study showed that the computed moment magnitudes for seismic events were insignificant to warrant safety concerns. The study findings showed the usefulness of coupled hydromechanical models in predicting the subsurface stress changes associated with CO 2 injection. The knowledge gained from this study will serve as a guideline for industries planning to undertake underground CO 2 storage, and characterize the subsurface stress changes.

02 PETROLEUM↗
NASA NTRS2022

A Radiative Transfer Simulator for PACE: Theory and Applications

A radiative transfer simulator was developed to compute the synthetic data of all three instruments onboard NASA’s Plankton Aerosol, Cloud, ocean Ecosystem (PACE) observatory, at the top of the atmosphere (TOA). The instrument suite includes the ocean color instrument (OCI), the HyperAngular Rainbow Polarimeter 2 (HARP2), and the Spectro-Polarimeter for Planetary Exploration 1 (SPEXone). The PACE simulator is wrapped around a monochromatic radiative transfer model based on the successive order of scattering (RTSOS), which accounts for atmosphere and ocean coupling, polarization, and gas absorption. Inelastic scattering, including Raman scattering from pure ocean water, fluorescence due to chlorophyll, and colored dissolved organic matter (CDOM), is also simulated. This PACE simulator can be used to explore the sensitivity of the hyperspectral and polarized reflectance of the Earth system with tunable atmosphere and ocean parameters, which include aerosol and cloud number concentration, refractive indices, and size distribution, ocean particle microphysical parameters, and solar and sensor-viewing geometry. The PACE simulator is used to study two important case studies. One is the impact of the significant uncertainty in pure ocean water absorption coefficient to the radiance field in the ultraviolet (UV) spectral region, which can be as much as 6%. The other is the influence of different amounts of brown carbon aerosols and CDOM on the polarized radiance field at TOA. The percentage variation of the radiance field due to CDOM is mostly for wavelengths smaller than 600 nm, while brown aerosols affect the whole spectrum from 350 to 890 nm, primarily due to covaried soot aerosols. Both case studies are important for aerosol and ocean color remote sensing and have not been previously reported in the literature.

PACE↗
DOE OSTI2026

Semi-inclusive deep-inelastic scattering on a polarized spin-1 target. II. Deuteron and spectator nucleon tagging

We develop the theoretical framework for semi-inclusive deep-inelastic scattering on a polarized spin-1 target and apply it to scattering on the polarized deuteron with spectator nucleon tagging. In Part I (previous article), we present the general form of the semi-inclusive cross section and polarization observables for the spin-1 target. In Part II (this article), we consider deep-inelastic scattering on the polarized deuteron with spectator nucleon tagging as a special case of target fragmentation. Methods of light-front quantization are employed to separate nuclear and hadronic structure in the high-energy process and achieve a composite description. The light-front wave function of the polarized deuteron is obtained from a rotationally covariant three-dimensional wave function in the center-of-mass frame of the proton-neutron system. The tagged structure functions are computed in the impulse approximation. The momentum and spin distribution of the active nucleon are controlled by the deuteron polarization and the detected spectator momentum (𝐷/𝑆 wave ratio). The cross section and spin asymmetries are evaluated for general deuteron polarization (vector and tensor, longitudinal and transverse) as functions of the spectator momentum. Tensor-polarized spin asymmetries of order unity are achieved for spectator momenta of approximately 300 MeV, which select configurations with a large 𝐷 wave. Sum rules for the tagged spin structure functions are derived. The results can be used for simulations of spectator tagging in future polarized fixed-target experiments (Jefferson Lab) or at the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Linear Boltzmann transport for jet propagation in the quark-gluon plasma: Inelastic processes and jet modification

A linear Boltzmann transport (LBT) Monte Carlo model has been developed to describe jet propagation and interaction with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions. A complete set of elastic-scattering processes and medium-induced gluon emissions based on the higher-twist formalism are incorporated for both jet shower and medium recoil partons. It has been employed to describe experimental data on large transverse momentum hadron and jet spectra, correlation and jet substructures in high-energy heavy-ion collisions. Here, we document in detail the structure of the model and validation of the Monte Carlo implementations of the physics processes in LBT, in particular, the inelastic process of medium-induced gluon radiation. We carry out a comprehensive examination of the jet-medium interaction as implemented in LBT through energy loss and momentum broadening of a single hard parton, the energy and transverse momentum transfer from leading partons to medium-induced gluons and jet-induced medium excitation, and medium modification of reconstructed jets in a static and uniform medium. With realistic and event-by-event hydrodynamic medium in heavy-ion collisions, we compute and compare with experimental data on the jet cone-size dependence of the single inclusive jet suppression at both the BNL Relativistic Heavy-Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), the dijet asymmetry at the LHC and γ-jet correlation at RHIC. Effects of medium-induced gluon emissions and jet-induced medium excitation on jet observables are systematically examined. Rescatterings of the radiated gluons and recoil partons with the QGP are found essential to account for the enhancement of soft particle yield toward the edge of the jet cone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Anisotropic inelastic materials milestone 2021

The focus on new fuel and cladding models are drivers for several new capabilities at the MOOSE and BISON levels. Since these materials are anisotropic, they require an extension of our inelastic mechanics capabilities, which currently heavily rely on radial return mapping and the assumption of isotropic stiffness. We plan to add additional code paths to support fully anisotropic materials, without degrading the performance for modeling isotropic materials.

42 ENGINEERING↗
DOE OSTI2025

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Software Verification and Validation Guidelines for Non-Linear Soil-Structure Interaction Analysis

Seismic analysis of structures, systems, and components (SSCs), including the consideration of soil-structure interaction (SSI) effects, is an important and required step in the design and licensing of nuclear power plant SSCs important to safety. Historically, the SSI analysis of nuclear structures has been performed using equivalent linear methods. However, there has been considerable industry investment in alternative seismic design and analysis approaches to reduce the construction cost of new reactors. Toward that goal and in alignment with the Licensing Modernization Project (LMP) framework, the Nuclear Regulatory Commission (US NRC) has proposed a risk-informed, performance-based approach to seismic design that allows inelastic response in those nuclear plant structures not required for confinement. As a complementary effort, reactor designers are exploring nonlinear seismic analysis methods to optimize structural designs and reduce construction costs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Improved control of electron computer-generated holographic grating groove profiles using ion beam gas-assisted etching

In a transmission electron microscope, electrons are described by matter-waves with wavelengths five orders of magnitude smaller than optical electromagnetic waves. Analogous to optical holography, electron wavefronts can be shaped using nanoscale holographic gratings. Here we demonstrate a novel, scalable nanofabrication method for creating off-axis holographic gratings that demonstrate near ideal diffraction efficiencies for binary, sinusoidal, and blazed grating groove profiles. We show that this method can produce up to 50 µm diameter area gratings that diffract up to 68% of the transmitted electron wave into a desired diffraction order with less than 7% into any other order. Additionally, we find that the amount of inelastically scattered electrons from the material gratings remaining in the coherent diffraction orders from the gratings is negligible in the far field.

Johnson, Cameron W.↗
DOE OSTI2022

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations G(r,t), a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor S(k,ω). In this work, we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in S(k,ω) going from the noninteracting limit (U = 0) to strong coupling (U→∞), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find S(k,ω), from which QFI is then calculated. We find that QFI grows with U. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above U=2.5 (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth W(U) and the half-chain von Neumann entanglement entropy. We compute G(r,t) by Fourier transforming S(k,ω). The results indicate a crossover in the short-time short-distance dynamics at low U characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large U featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by U=3. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the G(r,t) dynamics found at low U.

1-dimensional systems↗
NASA NTRS1979

The adiabatic semiclassical perturbation theory for vibrationally inelastic scattering. I - Collinear calculations. II - Three-dimensional treatment

A semiclassical approximation to treat vibrationally inelastic scattering is developed. The vibrational basis set used is adiabatic with respect to a reference potential which is chosen to be as close as possible to the true potential and also gives easily obtainable solutions to the vibrational wave equation. The radial wave functions are obtained using the WKB approximation, and the coupled Schroedinger equations are solved by a first-order perturbation method to yield a phase shift matrix which is exponentiated to give the full scattering matrix. Results were obtained for all the cases computed by Secrest and Johnson and by Clark and Dickinson, and the agreement is better than 10% for half of the cross-sections and rarely off by more than a factor of 2.

Cross, R. J., Jr.↗