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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A collinear perspective on the Regge limit

The high energy (Regge) limit provides a playground for understanding all loop structures of scattering amplitudes, and plays an important role in the description of many phenomenologically relevant cross-sections. While well understood in the planar limit, the structure of non-planar corrections introduces many fascinating complexities, for which a general organizing principle is still lacking. We study the structure of multi-reggeon exchanges in the context of the effective field theory for forward scattering, and derive their factorization into collinear operators (impact factors) and soft operators. We derive the structure of the renormalization group consistency equations in the effective theory, showing how the anomalous dimensions of the soft operators are related to those of the collinear operators, allowing us to derive renormalization group equations in the Regge limit purely from a collinear perspective. The rigidity of the consistency equations provides considerable insight into the all orders organization of Regge amplitudes in the effective theory, as well as its relation to other approaches. Along the way we derive a number of technical results that improve the understanding of the effective theory. We illustrate this collinear perspective by re-deriving all the standard BFKL equations for two-Glauber exchange from purely collinear calculations, and we show that this perspective provides a number of conceptual and computational advantages as compared to the standard view from soft or Glauber physics. We anticipate that this formulation in terms of collinear operators will enable a better understanding of the relation between BFKL and DGLAP in gauge theories, and facilitate the analysis of renormalization group evolution equations describing Reggeization beyond next-to-leading order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Phase-Field Modeling of Damage Evolution in Ceramic Matrix Composite (CMC) and Environmental Barrier Coating (EBC)

Ceramic matrix composites (CMCs) protected by environmental barrier coatings (EBCs) present a promising materials solution for next generation gas turbines. Developments of more robust and efficient EBCs and mechanically tougher CMCs are thus of significant technological importance. Here we develop a phase-field modeling framework that incorporates the thermally grown oxide (TGO), recognized as a critical factor for degradation and failure of EBCs. We simulate crack growth in the TGO and the potential extension into the bond coat / CMC substrate. The model efficiently takes account of the large inelastic deformation induced by the severe volume expansion of TGO, thanks to our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm. A phase-field model is built for damage evolution in CMCs including crack growth and interfacial sliding. The effects of fiber layout and interfacial sliding on the macroscopic toughness of CMCs are revealed by large-scale simulations and compared to experiments.

advanced energy systems and materials↗
DOE OSTI2025

Direct NeTS sampling of nuclear graphite $S(α, β, T)$ in Serpent

For advanced reactor applications, Neural Thermal Scattering (NeTS) modules were developed to predict the thermal scattering law (TSL or $S(α, β, T)$) of a nuclear graphite neutron moderator. NeTS are multi-layer, feedforward artificial neural networks, which act as universal function approximators designed for TSL datasets. In this case, a 4-layer neural network with 164 neurons per layer is trained using FLASSH evaluated data in PyTorch and serialized as a torchscript dictionary to predict $S(α, β, T)$ on-the-fly. Relative, absolute and maximum percent deviations of NeTS from File 7 data generated using the FLASSH code are on the order of 0.01%, 0.1% and 1%, respectively, with low inference latencies of 0.000172 s per $S(α, β, T)$ at a given temperature. Capturing the full dimensionality of possible inelastic neutron-lattice interactions, NeTS functionality is embedded in the Serpent Monte Carlo code, where $S(α, β, T)_{NeTS}$ sampling is conducted on-the-fly and compared to ACE look-up-tables for predicting TREAT criticality. k-eff differences between sampling algorithms of 6 pcm are observed and are within the order of Monte Carlo uncertainty. Compared to discrete and continuous-energy ACE files (30 MB and 131 MB per temperature), the NeTS format is on the order of 200–300 kB for a continuous-temperature, interpolation-free representation of $S(α, β, T)$ and cross sections. NeTS-in-Serpent runtimes comparable with ACE look-up tables are achieved by scaling NeTS for high performance computing architectures with hybrid OpenMP + MPI parallelization. This work validates a novel, self-contained reactor physics framework for predictive cross sections, and demonstrates a general methodology for embedding modern machine learning libraries within existing neutronic analysis frameworks.

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2024

Universality of Collective Density Fluctuations in Liquids at and away from Equilibrium – An Integrated Neutron Scattering, Theoretical, and Computational Study (Final Technical Report)

Liquids, ubiquitous on earth, are prototypical disordered condensed matter. As the mediate phase between non-structured gases and ordered solids, liquids p ossess intriguing complexities that are absent on either extreme counter phase, especially when they are driven away from equilibrium. Despite decades of studies, the physics of liquids still keeps surprising us and challenges our understanding of condensed matter, emergent phenomena, and complex systems. Among the many interesting properties, the nature of collective dynamics in liquids remains elusive because of both the presence of strong interactions and the absence of translational invariance. This project studied the universal features of collective density fluctuations (or collective modes), which include both phonon-like excitations and collective relaxations, in the Q-dependent generalized hydrodynamic regime in liquids and liquid-like matter both at and away from equilibrium. We performed systematic studies of the fast collective dynamics (including both phonon-like excitations and collective relaxations) of four classes of liquids characterized by different fragilities and interactions using synergistically integrated coherent neutron scattering experiments (Inelastic Neutron Scattering (INS), Quasi-Elastic Neutron Scattering (QENS), and Neutron Spin Echo (NSE)), a ViscoElastic Hydrodynamics theory we have been developing, andMolecular Dynamics (MD) simulations and our analysis package LiquidLib. The transformative knowledge of the collective dynamics in liquids is important to understand and ultimately control the transport of mass, energy, and charge in liquids, both at and away from thermal equilibrium, which may will likely impact a wide range of disciplines.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic Scale Crystal Field Mapping of Polar Vortices in Oxide Superlattices

Polar vortices in oxide superlattices can be utilized as potential candidates for data storage applications due to their unique polarization topologies. The structure and dipole arrangement in polar vortices has been studied via X-ray scattering techniques, (scanning) transmission electron microscopy ((S)/TEM) and computational calculations. However, the fundamental correlation between the atomic structure and the electronic structure (which is manifested in the chemical bonding) has heretofore not been explored. Further, The hybridization between nominally empty $\textit{d}$ orbitals on the B-site with the occupied O 2$\textit{p}$ orbitals favors the condensation of a polar (ferroelectric) state in $\textit{AB}$O 3 perovskite oxides. The complex, continuously rotating local polarization texture of the vortices, in turn, can result in especially intricate $\textit{d}$- orbital interactions. Soft X-ray spectroscopy can probe these interactions at the transition metal $\textit{L}$-edge, but these techniques do not have the spatial resolution to resolve variations within one vortex (~5 nm region). Electron energy loss spectroscopy (EELS) in the STEM mode uses inelastically scattered electrons to probe the core-shell excitations (empty density of states) of transition metals at atomic resolution.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

TNSL support in GNDS 2.0 and beyond [Slides]

This presentation begins by discussing the TNSL format options that went through a major overhaul in GNDS-2.0 and it examines the changes in 2.0. Additionally, it discusses three issues with further changes that should be considered. The first issue is that the project needs some guidance on what to expect when evaluations are performed with coherent inelastic. The second issue is that GNDS-2.0 does not provide a way to clearly specify in the evaluation how to switch to ‘standard’ incident neutron evaluations for energies or temperatures outside the TNSL domain. The third issue is that when GNDS-2.0 was designed, it was assumed that S(α,β) would always be given on a uniform interpolation grid. The presentation concludes by discussing how New JENDL-5 TNSL evaluations have some complications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Laser Rayleigh and Raman Diagnostics for Small Hydrogen/oxygen Rockets

Localized velocity, temperature, and species concentration measurements in rocket flow fields are needed to evaluate predictive computational fluid dynamics (CFD) codes and identify causes of poor rocket performance. Velocity, temperature, and total number density information have been successfully extracted from spectrally resolved Rayleigh scattering in the plume of small hydrogen/oxygen rockets. Light from a narrow band laser is scattered from the moving molecules with a Doppler shifted frequency. Two components of the velocity can be extracted by observing the scattered light from two directions. Thermal broadening of the scattered light provides a measure of the temperature, while the integrated scattering intensity is proportional to the number density. Spontaneous Raman scattering has been used to measure temperature and species concentration in similar plumes. Light from a dye laser is scattered by molecules in the rocket plume. Raman spectra scattered from major species are resolved by observing the inelastically scattered light with linear array mounted to a spectrometer. Temperature and oxygen concentrations have been extracted by fitting a model function to the measured Raman spectrum. Results of measurements on small rockets mounted inside a high altitude chamber using both diagnostic techniques are reported.

Degroot, Wilhelmus A.↗
DOE OSTI2024

Target mass corrections in lepton–nucleus DIS: Theory and applications to nuclear PDFs

Motivated by the wide range of kinematics covered by current and planned deep-inelastic scattering (DIS) facilities, we revisit the formalism, practical implementation, and numerical impact of target mass corrections (TMCs) for DIS on unpolarized nuclear targets. An important aspect is that we only use nuclear and later partonic degrees of freedom, carefully avoiding a picture of the nucleus in terms of nucleons. After establishing that formulae used for individual nucleon targets (p,n), derived in the Operator Product Expansion (OPE) formalism, are indeed applicable to nuclear targets, here we rewrite expressions for nuclear TMCs in terms of re-scaled (or averaged) kinematic variables. As a consequence, we find a representation for nuclear TMCs that is approximately independent of the nuclear target. We go on to construct a single-parameter fit for all nuclear targets that is in good numerical agreement with full computations of TMCs. We discuss in detail qualitative and quantitative differences between nuclear TMCs built in the OPE and the parton model formalisms, as well as give numerical predictions for current and future facilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1982

The 3-D inelastic analysis methods for hot section components: Brief description

Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate yet more cost-effective analysis of components subjected to severe thermal gradients and loads in the presence of mechanical loads, with steep stress and strain gradients are being developed. Anisotropy, time and temperature dependent plasticity and creep effects are also addressed. The approach is to develop four different theories, one linear and three higher order theories (polynomial function, special function, general function). The theories are progressively more complex from linear to general function in order to provide streamlined analysis capability with increasing accuracy for each hot section component and for different parts of the same component according to the severity of the local stress, strain and temperature gradients associated with hot spots, cooling holes and surface coating cracks. To further enhance the computational effectiveness, the higher order theories will have embedded singularities (cooling passages, for example) in the generic modeling region. Each of the four theories consists of three formulation models derivable from independent theoretical formulations. These formulation models are based on: (1) mechanics of materials; (2) special finite elements; and (3) an advanced formulation to be recommended by the contractor.

Chamis, C. C.↗
DOE OSTI2022

Proposed Generalized Header File for TSLs [Slides]

This presentation discusses several recent issues have been raised over the past few years regarding inconsistent definitions for the number of scattering atoms between incoherent elastic and inelastic (MT2 and MT4 files, respectively) in the same MF7 file and isotopic distribution of primary scattering atom. It also discusses how processing codes handle these potential issues in different ways. Additionally, it discusses how the preferred solution would be automation in order to reduce possibility of user error and a proposed solution in the form of a MF7/MT1 file is examined with preliminary examples. The presentation concludes with a question-and-answer discussion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Analysis of Organic Electron Materials using Neutron Scattering (Final Report)

Neutron scattering is a powerful suite of measurement techniques that enable unique analysis of structure and dynamics of a wide variety of materials. Specifically, neutrons interact strongly with protons, which enables excellent contrast in light element organic materials that cannot be obtained with any other measurements. This series of three proposals focused on using inelastic neutron scattering (INS) to study the vibrations (dynamics) of the organic molecules involved in organic electronic devices. Our main goals were to demonstrate that INS is capable of elucidating new, useful, and unique knowledge about organic electronic samples that could be used to improve the function of devices or to solve outstanding problems that could not be solved using other measurement techniques. Since all of the goal involved using INS or other neutron techniques to do something “new” a large part of the effort was focused on complex data interpretation and development of modeling workflows that could be used by other members of the organic electronics community to study their materials. To facilitate the goal of making INS measurements more accessible to a broader community and for a wider range of materials, we optimized and engineered the computational workflows to be well documented, publicly available, and user friendly. The individual titles of the three sequential DOE-BES awards are: (1) Engineering doping profiles in organic semiconducting materials, (2) Measurement and modeling of molecular dopant structure and dynamics, and (3) Quantification of Dynamic Disorder in Electronic and Excitonic Organic Materials. For readability, this report is organized to highlight the progress of the research thematically. Many of the ideas started in one grant period and then spilled over into later awards and collaborative research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A physics-informed multi-agents model to predict thermo-oxidative/hydrolytic aging of elastomers

This paper introduces a novel physics-informed multi-agents constitutive model to propose prediction in quasi-static constitutive behavior of cross-linked elastomer and the loss of mechanical performance during environmental aging. The presented model is used to simulate the effect of single-mechanism chemical aging (i.e. thermal-inducedor hydrolytic aging) on the behavior of the material in this hybrid framework. Those environmental single-mechanism damages change the polymer matrix over time due to massive chain scission, chain formations, and changing the arrangement of molecules in the polymer matrix. Here we propose a data-driven super-constrained machine-learned engine to represent damage in the polymer matrix and capture the changes in material behavior, including its inelastic features such as Mullins effect and permanent set in the course of aging. We have simplified the 3D stress–strain tensor mapping problem into a small number of super-constrained 1D mapping problems by means of a sequential order reduction. An assembly of multiple replicated conditional neural-network learning-agents (L-agents) is trained to systematically simplify the high-dimensional mapping problem into multiple 1D problems, each represented by a different type of agent. Our hybrid framework is designed to capture the effect of deformation history, aging time, and aging temperature. The model is validated with respect to a comprehensive set of experiments specifically designed to benchmark model capabilities and also against available data in the literature. Thermodynamic consistency and frame independency have been verified. Besides acceptable predictive abilities, a significant reduction of computational cost to predict behavior at multiple states of deformation is the most significant feature of this model.

42 ENGINEERING↗
DOE OSTI2022

Thermal Neutron Scattering Improvements and Fixes for MCNP6.3 [Slides]

Major error (10,000+ pcm) discovered when using uranium-dioxide or uranium-nitride included in ENDF80SaB2 release. Improvements are being made towards preparing the code for ENDF/B-VIII.1 (ACE) format updates. A minor error (0-10’s pcm) was discovered while working on new capability development for Lab Directed Research & Development (LDRD) project. Previously, only a single coherent or incoherent elastic channel (along with an inelastic channel) was allowed in a thermal neutron scattering evaluation and subsequent processed ACE file. Based on what is planned in ENDF/B-VIII.1 release, there will be (a) new thermal scattering evaluation(s) which uses this mixed-mode coherent and incoherent elastic scattering blocks. MCNP6.3 has been modified to be able to handle such data in the future. Some reorganization has been made in the cross section calculation (acetot, sabcol, colidn) in the presence of TSL data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Grain dynamics in zero gravity

The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. One interesting limit for which it was recently possible to construct a theory is that where the grain-grain interactions are dominated by binary collisions. The kinetic model of granular systems if similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicitly. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high-density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities, where an elastic description of the collision is still valid, is investigated. This will allow direct interpretation by dynamical computer simulations as well as by kinetic theory.

Werner, B. T.↗
NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗
NASA NTRS1985

Local strain redistribution corrections for a simplified inelastic analysis procedure based on an elastic finite-element analysis

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 proposes. 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 predicts stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load-cycled problems. Neuber-type corrections were derived and incorporated in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was used on a mechanically load-cycled benchmark notched-plate problem. The predicted material response agrees well with the nonlinear finite-element solutions for the problem. The simplified analysis computer program was 0.3% of the central processor unit time required for a nonlinear finite-element analysis.

Kaufman, A.↗
DOE OSTI2023

Capturing dynamical correlations using implicit neural representations

Abstract Understanding the nature and origin of collective excitations in materials is of fundamental importance for unraveling the underlying physics of a many-body system. Excitation spectra are usually obtained by measuring the dynamical structure factor, S ( Q , ω ), using inelastic neutron or x-ray scattering techniques and are analyzed by comparing the experimental results against calculated predictions. We introduce a data-driven analysis tool which leverages ‘neural implicit representations’ that are specifically tailored for handling spectrographic measurements and are able to efficiently obtain unknown parameters from experimental data via automatic differentiation. In this work, we employ linear spin wave theory simulations to train a machine learning platform, enabling precise exchange parameter extraction from inelastic neutron scattering data on the square-lattice spin-1 antiferromagnet La 2 NiO 4 , showcasing a viable pathway towards automatic refinement of advanced models for ordered magnetic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetization dynamics fingerprints of an excitonic condensate t 2 g 4 magnet

The competition between spin-orbit coupling (SOC) λ and electron-electron interaction U leads to a plethora of novel states of matter, extensively studied in the context of t 2 g 4 and t 2 g 5 materials, such as ruthenates and iridates. Excitonic magnets—the antiferromagnetic state of bounded electron-hole pairs-are prominent examples of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground state of SOC t 2 g 4 Hubbard models. Here we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the λ = 0 state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wave-vector q → 0 . In the excitonic condensation regime at large U , we also have found a novel high-energy π mode composed solely of orbital excitations. Overall, these features do not appear all together in any of the neighboring states in the phase diagram and thus constitute unique fingerprints of the t 2 g 4 excitonic magnet, of importance in the analysis of neutron and resonant inelastic x-ray scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗