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 379 records · Page 21

Large-momentum effective theory

Since the parton model was introduced by Feynman more than 50 years ago, much has been learned about the partonic structure of the proton through a large body of high-energy experimental data and dedicated global fits. However, limited progress has been made in calculating partonic observables such as the parton distribution function (PDFs) from the fundamental theory of strong interactions, quantum chromodynamics (QCD). Recently some advocated for a formalism, large-momentum effective theory (LaMET), through which one can extract parton physics from the properties of the proton traveling at a moderate boost factor such as γ ~ 2 – 5 . The key observation behind this approach is that Lorentz symmetry allows the standard formalism of partons in terms of light-front operators to be replaced by an equivalent one with large-momentum states and time-independent operators of a universality class. With LaMET, the PDFs, generalized PDFs or generalized parton distributions, transverse-momentum-dependent PDFs, and light-front wave functions can all be extracted in principle from lattice simulations of QCD (or other nonperturbative methods) through standard effective field theory matching and running. Future lattice QCD calculations with exascale computational facilities could help one to understand the experimental data related to the hadronic structure, including those from the upcoming electron-ion colliders dedicated to exploring the partonic landscape of the proton. In this work, the progress made in the past few years in the development of the LaMET formalism and its applications is reviewed, with an emphasis on a demonstration of its effectiveness from initial lattice QCD simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

Progress on Associate-Particle Imaging Algorithms, 2020

The present work describes progress on the development of imaging algorithms that use fast neutron signatures acquired using the associated-particle imaging (API) method. The present work complements ongoing work to develop neutron source and detector hardware to enable field inspection by investigating algorithms that are capable of discriminating between critical materials or extracting three-dimensional geometrical information from single-sided or transmission measurements. The present work is divided into three approaches:(1)Iterative reconstruction of inelastic gamma-ray emissions to perform three-dimensional time-of-flight imaging in a single view in either transmission or backscatter configurations. Iterative reconstruction enables image resolution better than the inherent TOF resolution.(2)Decomposition of registered neutron and x ray radiographs into an assumed material list for each pixel in the image.(3)Material identification using full spectral analysis that includes the emergent neutron and gamma ray energies, times, and angles.Progress for each approach is summarized for fiscal year 2020.

97 MATHEMATICS AND COMPUTING↗

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Abstract]

Historically, the free gas approximation has been used to treat the thermal scattering of neutrons with energies below a few electron-volts (eV) in unevaluated materials. However, this method inadequately reproduces neutron scattering at these energies. Until recently, only a limited number of materials had available thermal scattering law (TSL) files/libraries in the ENDF nuclear data libraries in this energy range. With advancements in atomistic modeling techniques, such as molecular dynamics, ab-initio molecular dynamics, and density functional theory, TSL libraries have become available for many more materials. This is particularly relevant due to the rising interest in several advanced reactor systems that require novel moderator and reflector materials. While quasi-integral and integral benchmarks have been designed to validate historically important moderator materials (such as light water and polyethylene), there is currently a lack of standard validation methods for TSLs, especially when multiple conflicting TSLs exist. To address this issue, the Oak Ridge National Lab Nuclear Data group has been working on utilizing inelastic neutron scattering (INS) measurements combined with transmission (i.e., total cross section) measurements to evaluate and validate TSLs for different materials. We plan to demonstrate how this method has worked on materials such as polyethylene, lucite, and polystyrene. In addition, we will compare the newly created libraries to ENDF libraries for these materials and explain why integral benchmarks should not be used for validation when multiple TSLs exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. In this work, the results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a phonon-based sampling method for thermal neutron scattering data

Simulations of reactor systems require access to accurate nuclear data. For many systems, thermal neutron scattering data can have large effects on the eigenvalue and neutron flux distributions. Inelastic thermal neutron scattering can excite or de-excite vibrational, rotational, and translational modes in a material, so thermal scattering evaluations are often obtained by summing over the number of phonons created/destroyed by a scattering event. In recent years, the thermal scattering cross sections and angular distributions have greatly improved in accuracy, but the format in which this data is delivered to simulation codes has remained virtually unchanged. Thermal scattering data is typically either compiled into large tables and sorted by incoming neutron energy, outgoing neutron energy, scattering angle, and material temperature, or represented as cumulative distribution functions of momentum exchange or energy exchange. Either method can be quite memory intensive when fine bins are used. In an effort to decrease the amount of space that processed thermal scattering data requires, an alternate format is proposed. The phonon-based sampling method introduced here can sample the number of phonons excited for each collision, the change in neutron energy, and the scattering angle while avoiding pre-computed angular bins and limiting the amount of data that is dependent on incoming energy. Through this method, the generation and storage of large interpolation tables is avoided, which could have benefits in both memory storage and accuracy. While the initial implementation of this method is slower than current alternatives, it is significantly more resistant to grid coarseness errors and has good potential for improvement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A super-resolution technique to analyze single-crystal inelastic neutron scattering measurements using direct-geometry chopper spectrometers

Direct-geometry time-of-flight chopper neutron spectroscopy is instrumental in studying dynamics in liquid, powder, and single crystal systems. We report here that real-space techniques in optical imagery can be adapted to obtain reciprocal-space super resolution dispersion for phonon or magnetic excitations from single-crystal neutron spectroscopy measurements. The procedure to reconstruct super-resolution energy dispersion of excitations relies on an accurate determination of the momentum and energy-dependent point spread function and a dispersion correction technique inspired by an image disparity calculation technique commonly used in stereo imaging. Here by, applying these methods to spinwave dispersion data from a virtual neutron experiment demonstrates ~5-fold improvement over nominal energy resolution.

47 OTHER INSTRUMENTATION↗

Optimization of an Energy Tuning Assembly for High Explosives Detection

The Portable Isotopic Neutron Spectroscopy (PINS) system, employs neutron-induced gamma-ray spectroscopy and provides a nondestructive method for high explosives detection. In standard operation it uses Californium-252 as a neutron source. Operating PINS with a deuterium-tritium (DT) neutron generator has some advantages over Cf-252, including lifetime and ability to produce high-energy inelastic scattering gamma rays. However, current systems using DT neutron generators suffer from a high environmental background and reduced ability to induce neutron capture, reducing spectral quality and limiting nitrogen sensitivity. Here, this study presents the development of an energy-tuning assembly (ETA) designed to optimize the DT neutron energy spectrum to increase nitrogen reaction rates in a target, thereby improving high explosive detection capabilities. A metaheuristic optimization framework, MultiGNOWEE, coupled with MCNP, was employed to generate two ETA configurations: a single-objective ETA, which maximizes nitrogen capture reactions, and a multi-objective ETA, which balances neutron capture and inelastic scattering. Simulations demonstrated the optimized configurations achieved up to a 10-fold improvement in nitrogen capture rates compared to the bare configuration. Experimental validation was conducted using a DT neutron generator and a high-purity germanium (HPGe) detector. Two prototype ETAs were constructed and assessed on a melamine simulant. Measurements demonstrated improved nitrogen detection for both prototype ETA configurations when compared to the standard system.

97 MATHEMATICS AND COMPUTING↗

Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components

Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.

42 ENGINEERING↗

Ammonia synthesis on BaTiO 2.5 H 0.5 : computational insights into the role of hydrides

Perovskite oxyhydrides such as BaTiO 2.5 H 0.5 have been found to be able to catalyze NH3 synthesis, but the mechanism and the role of the catalyst's lattice hydrides in the catalytic reaction remain unknown. Here we employ first principles density functional theory to investigate the mechanism of ammonia synthesis and the role of lattice hydrides on a prototypical perovskite oxyhydride, BaTiO 2.5 H 0.5 (BTOH). Two mechanistic hypotheses, the distal and alternating pathways, have been tested on the Ti2O 2 termination of the BTOH (210) surface, previously determined to be the most stable surface termination under the reaction conditions considered. In the distal pathway, H atoms hydrogenate N 2 to form the *N–NH x key intermediates, followed by N–N bond breaking. In the alternating pathway, H atoms hydrogenate N 2 in an alternating fashion to form the *NH x –NH y intermediates before N–N bond breaking and formation of co-adsorbed *NH x /*NH y on the surface. We find that the subsurface hydride vacancy formed after reaction of *N 2 with the lattice hydride is key to the distal pathway, leading to surface nitride formation after breaking the *N–NH 3 bond, while the neighboring surface Ti sites are key to bridging and stabilizing the *NNH intermediate in the alternating pathway. In both pathways, desorption of NH 3 is the most uphill in energy. Our results provide important insights into the role of hydrides and surface vacancies in hydrogenation reactions over BTOH, which will be useful to guide future spectroscopic experiments such as operando IR and inelastic neutron scattering to verify the key intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proceedings of RIKEN BNL Research Center Workshop: Small-x Physics in the EIC Era [Slides]

Understanding the high energy limit of hadronic and nuclear collisions is at the forefront of nuclear and particle physics. When boosted to ultrahigh energies, all hadrons and nuclei eventually transform into a universal form of matter called the Color Glass Condensate (CGC). The CGC is characterized by the high density (saturation) of small-x gluons which leads to distinct experimental signatures. Tantalizing hints of the CGC have been observed at HERA, RHIC and the LHC, but the prospects for the discovery of the CGC are more promising at future experiments such as the Electron-Ion Collider (EIC) at BNL Indeed, according to the National Science Academy report published in 2019, one of the three major goals of the EIC is to address the nature of the gluon saturation. (The other two are the mass and spin structure of the nucleons.) With this in mind, the small-x community is gearing up to meet the challenges of the EIC era. Over the past several years, there has been impressive progress in the next-to-leading order (NLO) calculations in the CGC framework of various observables such as single hadron production in proton-nucleus collisions, inclusive and exclusive dijet and trijet production in Deep Inelastic Scattering (DIS), jet-plus-photon production in DIS, etc. We expect that NLO calculations will be the standard tool to confront future experimental data at the EIC. Another emerging trend of the community is the interplay between small-x physics and spin physics. The RHIC result for the gluon helicity ΔG has underscored the necessity to understand the longitudinal spin structure of the proton at small-x. There have been theoretical indications that a significant fraction of spin and orbital angular momentum is stored in the small-x region. As for the transversely polarized proton, a surprising new connection between the gluon Sivers function at small-x and the QCD Odderon has been pointed out and its implications at the EIC has been discussed. In view of these developments, we think it is timely to organize a dedicated workshop on small-x physics to summarize the present status of the field and to discuss future directions. A major focus of this workshop will be to identify outstanding problems that could significantly benefit from collaborative efforts amongst scientists working on formal, phenomenological, and computational aspects of small-x physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The FLUKA code for space applications: recent developments

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Non-NASA Center↗

Warming nuclear pasta with dark matter: kinetic and annihilation heating of neutron star crusts

Neutron stars serve as excellent next-generation thermal detectors of dark matter, heated by the scattering and annihilation of dark matter accelerated to relativistic speeds in their deep gravitational wells. However, the dynamics of neutron star cores are uncertain, making it difficult at present to unequivocally compute dark matter scattering in this region. On the other hand, the physics of an outer layer of the neutron star, the crust, is more robustly understood. We show that dark matter scattering solely with the low-density crust still kinetically heats neutron stars to infrared temperatures detectable by forthcoming telescopes. We find that for both spin-independent and spin-dependent scattering on nucleons, the crust-only cross section sensitivity is 10 -43 - 10 -41 cm 2 for dark matter masses of 100 MeV - 1 PeV, with the best sensitivity arising from dark matter scattering with a crust constituent called nuclear pasta (including gnocchi, spaghetti, and lasagna phases). For dark matter masses from 10 eV to 1 MeV, the sensitivity is 10 -39 - 10 -34 cm 2 , arising from exciting collective phonon modes in a neutron superfluid in the inner crust. Furthermore, for any s-wave or p-wave annihilating dark matter, we show that dark matter will efficiently annihilate by thermalizing just with the neutron star crust, regardless of whether the dark matter ever scatters with the neutron star core. This implies efficient annihilation in neutron stars for any electroweakly interacting dark matter with inelastic mass splittings of up to 200 MeV, including Higgsinos. We conclude that neutron star crusts play a key role in dark matter scattering and annihilation in neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Testing the lepton content of the proton at HERA and EIC

Although protons are baryons with an overall vanishing lepton number, they possess a non-trivial leptonic content arising from quantum fluctuations which can be described by lepton parton distribution functions (PDFs) of the proton. These PDFs have been recently computed and can be used to define lepton-induced processes at high-energy colliders. In this article, we propose a novel way to test the computation of lepton PDFs of the proton by analyzing both non-resonant di-lepton and resonant Z gauge boson production processes induced by leptons within the proton at proton-electron colliders like HERA and EIC. Despite the fact that lepton PDFs of the proton are known to be small, this work demonstrates that both processes imply a measurable yield of events at HERA and EIC, which could be used to test these PDFs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Application of Gaussian Bayes classifier to differentiate chlorine-based chemical agents

The Portable Isotopic Neutron Spectroscopy (PINS) is a commercialized system developed by Idaho National Laboratory (INL) to examine chemical warfare agents (CWA) non-destructively, utilizing Prompt Gamma Neutron Activation Analysis (PGNAA) techniques. The PINS system takes advantage of a high-resolution gamma-ray spectrum from a mechanically-cooled high-purity germanium (HPGe) detector. One of the difficult technical challenges is to discriminate the chlorine-based chemical agents. Especially, CN, CNB, CNS and CG have similar chemical compositions to make it hard to discriminate them with a higher confidence. Current identification algorithms for PINS systems with 252-Cf sources have been improved and updated continuously as more field data became available, and new algorithms was studied to complement the current algorithms by adopting the Gaussian Bayes classifier. These new algorithms were intended to be applied to a subset of chlorine-based chemical agents, and their main goal is discriminate CN, CNB, CNS and CG with their ratios of the chlorine neutron inelastic 1763keV peak to the chlorine thermal neutron capture 1959keV peak, which is referred to as the “Cl i/c” or “clic” ratio in this study. The Cl i/c ratios were assumed to follow Gaussian distributions with the means and the standard deviations unique to their corresponding chemical agents. The prior probabilities of these four chemical agents were optimized with a collection of field data to achieve the best performance in terms of precision or positive predictive value (PPV). Finally, their posterior probabilities as functions of the Cl i/c ratio were implemented in the current version of PINS analysis software in order to be tested with more field data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Slides]

This presentation finds that INS measurements are not as variable based on sample composition, because phonon spectra is a properly of the bulk material. Additionally, INS and transmission measurements are and should be the most significant part of TSL validation, as they are the fundamental quantities used by Monte Carlo neutronics codes. Relaying solely on theoretical atomistic calculations for TSL evaluation is not enough, they are highly dependent on inter-atomic potentials (which are not the most accurate depending on the quantity trying to be reproduced) and evaluators knowledge and understanding of the material being studied. When dealing with Crit. benchmarks - While extremely useful for validation of nuclear data at all energies, critical benchmarks are not the best tool to provide a definitive answer on conflicting TSLs. Finally, INS and transmission measurements need to be the basis of validation of TSLs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗