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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Investigating the Accuracy of Water Models through the Van Hove Correlation Function

In this study, we present molecular-simulation-based calculations of the Van Hove correlation function (VHF) of water using multiple modeling approaches: classical molecular dynamics with simple three-site nonpolarizable models, with a polarizable model, and with a reactive force field; density functional tight-binding molecular dynamics; and ab initio molecular dynamics. Due to the many orders of magnitude difference in the computational cost of these approaches, we investigate how small and short the simulations can be while still yielding sufficiently accurate and interpretable results for the VHF. We investigate the accuracy of the different models by comparing them to recently published inelastic X-ray scattering measurements of the VHF. We find that all of the models exhibit qualitative agreement with the experiments, and in some models and for some properties, the agreement is quantitative. This work lays the foundation for future simulation approaches to calculating the VHF for aqueous solutions in bulk and under nanoconfinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1972

Evaluation of solution procedures for material and/or geometrically nonlinear structural analysis by the direct stiffness method.

This paper presents an assessment of the solution procedures available for the analysis of inelastic and/or large deflection structural behavior. A literature survey is given which summarized the contribution of other researchers in the analysis of structural problems exhibiting material nonlinearities and combined geometric-material nonlinearities. Attention is focused at evaluating the available computation and solution techniques. Each of the solution techniques is developed from a common equation of equilibrium in terms of pseudo forces. The solution procedures are applied to circular plates and shells of revolution in an attempt to compare and evaluate each with respect to computational accuracy, economy, and efficiency. Based on the numerical studies, observations and comments are made with regard to the accuracy and economy of each solution technique.

Stricklin, J. A.↗

DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗

NASA NTRS1974

Boltzmann equation for the electron gas of a nondegenerate plasma

The collision terms in the Boltzmann equation associated with various processes are derived. For processes having a Fokker-Planck (F-P) limit, the associated F-P operator is derived by means of physical arguments to determine the form of the operator; its multiplying constant is fixed by calculating the total energy exchange rate and comparing with the rate computed by other means. In this manner, the F-P operator is derived for electron-ion scattering, electron-electron scattering in the high-velocity limit, electron-atom elastic scattering, Compton scattering, and the high-velocity limit for inelastic scattering. Other processes considered are bremsstrahlung, radiative-recombination, photoionization, collisional ionization of atoms, and suprathermal-particle ionization of atoms.

Gould, R. J.↗

DOE OSTI2023

Neutron Leakage Spectra of the EUCLID Experiment [Abstract]

Integral experiments with sub-critical and critical configurations of special nuclear material are performed in support of nuclear data validation and adjustment. Different nuclear data evaluations may have different values for individual cross sections due to uncertainties in (or lack of) differential experiments, but compensating errors in these data sets can lead to the same k eff results for one application while vastly different for another application. One example of this is 239 Pu, where both ENDF/B-VIII.0 and JEFF-3.3 correctly compute k eff of the Jezebel critical assembly, but the individual contributions from each reaction are vastly different. To help resolve this specific case, the Experiments Underpinned by Computational Learning for Improvements in Nuclear Data (EUCLID) project used machine learning to design a set of experiments to help resolve the compensating errors in 239 Pu. A total of six responses were measured during the experimental campaign, which constrain the data in ways that k eff alone cannot and will be used for adjustment of the nuclear data. One of these responses is the neutron leakage spectrum, which recent work has shown to be useful for constraining the prompt fission neutron spectrum and inelastic scattering. The neutron leakage spectra were measured utilizing a 3 in. right cylinder EJ-301D detector. The measured signal in the detector was deconvoluted using spectrum unfolding techniques, which are presented and compared to simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Laser diagnostics for small rockets

Two nonintrusive flowfield diagnostics based on spectrally-resolved elastic (Rayleigh) and inelastic (Raman) laser light scattering were developed for obtaining local flowfield measurements in low-thrust gaseous H2/O2 rocket engines. The objective is to provide an improved understanding of phenomena occurring in small chemical rockets in order to facilitate the development and validation of advanced computational fluid dynamics (CFD) models for analyzing engine performance. The laser Raman scattering diagnostic was developed to measure major polyatomic species number densities and rotational temperatures in the high-density flowfield region extending from the injector through the chamber throat. Initial application of the Raman scattering diagnostic provided O2 number density and rotational temperature measurements in the exit plane of a low area-ratio nozzle and in the combustion chamber of a two-dimensional, optically-accessible rocket engine. In the low-density nozzle exit plane region where the Raman signal is too weak, a Doppler-resolved laser Rayleigh scattering diagnostic was developed to obtain axial and radial mean gas velocities, and in certain cases, H2O translational temperature and number density. The results from these measurements were compared with theoretical predictions from the RPLUS CFD code for analyzing rocket engine performance. Initial conclusions indicate that a detailed and rigorous modeling of the injector is required in order to make direct comparisons between laser diagnostic measurements and CFD predictions at the local level.

Zupanc, Frank J.↗

DOE OSTI2023

SIDIS-RC EvGen: A Monte-Carlo event generator of semi-inclusive deep inelastic scattering with the lowest-order QED radiative corrections

SIDIS-RC EvGen is a C++ standalone Monte-Carlo event generator for studies of semi-inclusive deep inelastic scattering (SIDIS) processes at medium to high lepton beam energies. In particular, the generator contains binary and library components for generating SIDIS events and calculating cross sections for unpolarized or longitudinally polarized beam and unpolarized, longitudinally or transversely polarized target. The structure of the generator incorporates transverse momentum-dependent parton distribution and fragmentation functions, whereby we obtain multi-dimensional binned simulation results, which will facilitate the extraction of important information about the three-dimensional nucleon structure from SIDIS measurements. In order to build this software, we have used recent elaborate QED calculations of the lowest-order radiative effects, applied to the leading order Born cross section in SIDIS. Here, in this paper, we provide details on the theoretical formalism as well as the construction and operation of SIDIS-RC EvGen, e.g., how we handle the event generation process and perform multi-dimensional integration. We also provide example programs, flowcharts, and numerical results on azimuthal transverse single-spin asymmetries.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Threshold resummation for computing large-x parton distribution through large-momentum effective theory

Parton distribution functions (PDFs) at large x are poorly constrained by high-energy experimental data, but extremely important for probing physics beyond standard model at colliders. We study the calculation of PDFs at large-x through large-momentum P z expansion of the lattice quasi PDFs. Similar to deep-inelastic scattering, there are two distinct perturbative scales in the threshold limit where the matching coefficient can be factorized into a space-like jet function at scale P z |1 – y| and a pair of heavy-light Sudakov form factors at scale P z . The matching formula allows us to derive a full renormalization group resummation of large threshold logarithms, and the result is consistent with the known calculation to the next-to-next to leading order (NNLO). This paves the way for direct large-x PDFs calculations in lattice QCD. As by-products, we find that the space-like jet function is related to a time-like version calculated previously through analytic continuation, and the heavy-light Sudakov form factor, calculated here to NNLO, is a universal object appearing as well in the large momentum expansion of quasi transverse-momentum-dependent PDFs and quasi wave-function amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A Multiphysics Multiscale Simulation Platform for Damage, Environmental Degradation, and Life Prediction of CMCs in Extreme Environments

This project successfully developed a multiphysics, multiscale computational framework to enhance the design and development of CMCs, with a focus on modeling highly nonlinear, time-dependent damage mechanisms and material degradation under extreme conditions, such as those experienced in turbine service environments. The project made significant advances in improving our understanding of progressive damage, oxidative degradation, and time-dependent inelastic deformation in CMCs, with particular attention to the role of uncertainties in predictions. Key outcomes include the integration of advanced material characterization, uncertainty quantification, and multiphysics constitutive models to predict the behavior of CMCs over their service life. A novel multiscale methodology was employed, which integrated microscale constituent behaviors with structural-scale responses, enabling the manufacturing defects in the microstructure that are prone to damage nucleation. Through the development of DL algorithms, the project advanced the prediction of damage initiation and crack propagation, taking into account the defect morphology and statistical variations across multiple scales. The framework was rigorously validated using thermomechanical experiments, which tested CMCs under various mechanical loadings at elevated temperatures, further enhancing the model's predictive capability. Overall, the research outcomes have provided a more accurate, reliable method for predicting CMC component life, significantly advancing material design, and improving component reliability in extreme environments. This work has strong implications for the optimization of turbine components and other high-performance applications where CMCs are used.

03 NATURAL GAS↗

DOE OSTI2024

Elucidating the Role of Dimensionality on the Electronic Structure of the Van der Waals Antiferromagnet NiPS 3

Abstract The sustained interest in investigating magnetism in the 2D limit of insulating antiferromagnets is driven by the possibilities of discovering, or engineering, novel magnetic phases through layer stacking. However, due to the difficulty of directly measuring magnetic interactions in 2D antiferromagnets, it is not yet understood how intra layer magnetic interactions in insulating , strongly correlated, materials can be modified through layer proximity. Herein, the impact of reduced dimensionality in the model van der Waals antiferromagnet NiPS 3 is explored by measuring electronic excitations in exfoliated samples using Resonant Inelastic X‐ray Scattering (RIXS). The resulting spectra shows systematic broadening of NiS 6 multiplet excitations with decreasing layer count from bulk down to three atomic layers (3L). It is shown that these trends originate from a decrease in transition metal‐ligand and ligand–ligand hopping integrals, and by charge‐transfer energy evolving from Δ = 0.83 eV in the bulk to 0.37 eV in 3L NiPS 3 . Relevant intralayer magnetic exchange integrals computed from the electronic parameters exhibit a decrease in the average interaction strength with thickness. This study underscores the influence of inter layer electronic interactions on intra layer ones in insulating magnets, indicating that magnetic Hamiltonians in few‐layer insulating magnets can greatly deviate from their bulk counterparts.

36 MATERIALS SCIENCE↗

NASA NTRS1980

Neutrino, gamma-ray, electron, and positron production in an ultrarelativistic plasma

Neutrino, gamma-ray, electron and positron production resulting from inelastic proton-proton collisions in a highly relativistic plasma such as may exist in extragalactic radio or gamma-ray burst sources is examined. The source functions of primary (pions, kaons, and neutrons) and secondary (photons, electrons, positrons and neutrinos) products of relativistic nuclear collisions are computed for the cases of power law and Maxwellian relativistic proton distributions. It is shown that in plasma which is optically thin to interactions between the plasma and secondary gamma-rays, electrons and positrons, only a small fraction of the initial energy is emitted in the observable form of neutrinos and gamma rays. In an optically thick plasma on the other hand, most of the energy of the relativistic protons is found to be equally divided between gamma rays and neutrinos, although only the neutrinos may escape freely to be observed.

Marscher, A. P.↗

DOE OSTI2023

Competing incommensurate spin fluctuations and magnetic excitations in infinite-layer nickelate superconductors

The recently discovered infinite-layer nickelates show great promise in helping to disentangle the various cooperative mechanisms responsible for high-temperature superconductivity. However, lack of antiferromagnetic order in the pristine nickelates presents a challenge for connecting the physics of the cuprates and nickelates. Here, by using a quantum many-body Green’s function-based approach to treat the electronic and magnetic structures, we unveil the presence of many two- and three-dimensional magnetic stripe instabilities that are shown to persist across the phase diagram of LaNiO 2 . Our analysis indicates that the magnetic properties of the infinite-layer nickelates are closer to those of the doped cuprates, which host a stripe ground state, rather than the undoped cuprates. The computed longitudinal-spin, transverse-spin, and charge spectra of LaNiO 2 are found to contain an admixture of contributions from localized and itinerant carriers. Theoretically obtained dispersion of magnetic excitations (spin-flip) is found to be in good accord with the results of recent resonant inelastic X-ray scattering experiments. Our study gives insight into the origin of strong magnetic competition in the infinite-layer nickelates and their relationship with the cuprates.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Nonlinear structural analysis of cylindrical thrust chambers using viscoplastic models

The results are presented of a viscoplastic stress-strain analysis of a cylindrical thrust chamber used in experiments that were designed to simulate the SSME operating conditions. The inelastic strain was calculated by using a viscoplastic model, and a quasi-three-dimensional structural analysis was performed by using a finite element program MARC. The temperatures and pressures were calculated by using the loading cycles of experiments, and these temperatures and pressures were used in the computations. The deformed shape of the component was predicted after the end of each loading cycle. The predicted shape qualitatively replicated the deformed shape of the component as observed in experiments. The results indicate that the use of viscoplastic models for structural analysis may lead to more realistic life assessments of experimental thrust chambers.

Arya, Vinod K.↗

NASA NTRS2014

A Technique for Mapping Characteristic Lengths to Preserve Energy Dissipated via Strain Softening in a Multiscale Analysis

It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).

softening↗

NASA NTRS1989

Thermo-viscoplastic analysis of hypersonic structures subjected to severe aerodynamic heating

A thermoviscoplastic computational method for hypersonic structures is presented. The method employs unified viscoplastic constitutive model implemented in a finite element approach for quasi-static thermal-structural analysis. Applications of the approach to convectively cooled hypersonic structures illustrate the effectiveness of the approach and provide insight into the transient inelastic structural behavior at elevated temperatures.

Thornton, Earl A.↗

NASA NTRS1990

Time-dependent wave-packet method for the complete determination of S-matrix elements for reactive molecular collisions in three dimensions

An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.

Judson, Richard S.↗