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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The 3D inelastic analysis methods for hot section components

A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗
NASA NTRS1975

Vibration-translation energy transfer in anharmonic diatomic molecules. 2: The vibrational quantum number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom was used to predict the variation of thermally averaged vibration-translation rate coefficients with temperature and initial-state quantum number. Multiple oscillator states were included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model was also used as a basis for evaluating several less complete but analytic models. Two computationally simple analytic approximations were found that successfully reproduced the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations were also identified. The relative rates of multiple-quantum transitions from excited states were evaluated for several molecular types.

Mckenzie, R. L.↗
NASA NTRS1976

Vibration-translation energy transfer in anharmonic diatomic molecules. II - The vibrational quantum-number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom is used to predict the variation of thermally averaged vibrational-translational rate coefficients with temperature and initial-state quantum number. Multiple oscillator states are included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model is also used as a basis for evaluating several less complete, but analytic, models. Two computationally simple analytic approximations are found that successfully reproduce the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations are identified, and the relative rates of multiple-quantum transitions from excited states are evaluated for several molecular types.

Mckenzie, R. L.↗
NASA NTRS1986

On 3D inelastic analysis methods for hot section components

Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗
NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗
NASA NTRS1985

Cyclic structural analyses of anisotropic turbine blades for reusable space propulsion systems

Turbine blades for reusable space propulsion systems are subject to severe thermomechanical loading cycles that result in large inelastic strains and very short lives. These components require the use of anisotropic high-temperature alloys to meet the safety and durability requirements of such systems. To assess the effects on blade life of material anisotropy, cyclic structural analyses are being performed for the first stage high-pressure fuel turbopump blade of the space shuttle main engine. The blade alloy is directionally solidified MAR-M 246 alloy. The analyses are based on a typical test stand engine cycle. Stress-strain histories at the airfoil critical location are computed using the MARC nonlinear finite-element computer code. The MARC solutions are compared to cyclic response predictions from a simplified structural analysis procedure developed at the NASA Lewis Research Center.

Manderscheid, J. M.↗
NASA NTRS1985

Cyclic structural analyses of anisotropic turbine blades for reusable space propulsion systems

Turbine blades for reusable space propulsion systems are subjects to severe thermomechanical loading cycles that result in large inelastic strains and very short lives. These components require the use of anisotropic high temperature alloys to meet the safety and durability requirements of such systems. To assess the effects on blade life of material anisotropy, cyclic structural analyses are being performed for the first stage high pressure fuel turbopump blade (HPFTB) of the space shuttle main engine (SSME). The blade alloy is directionally solidified MAR-M 246 alloy. The analyses are based on a typical test stand engine cycle. Stress-strain histories at the airfoil critical location are computed using the MARC nonlinear finite element computer code.

Manderscheid, J. M.↗
NASA NTRS2001

The Applicability of the Generalized Method of Cells for Analyzing Discontinuously Reinforced Composites

The paper begins with a short overview of the recent work done in the field of discontinuous reinforced composites, focusing on the different parameters which influence the material behavior of discontinuous reinforced composites, as well as the various analysis approaches undertaken. Based on this overview it became evident that in order to investigate the enumerated effects in an efficient and comprehensive manner, an alternative approach to the computationally intensive finite-element based micromechanics approach is required. Therefore, an investigation is conducted to demonstrate the utility of utilizing the generalized method of cells (GMC), a semi-analytical micromechanics-based approach, to simulate the elastic and elastoplastic material behavior of aligned short fiber composites. The results are compared with simulations using other micromechanical based mean field models and finite element (FE) unit cell models found in the literature given elastic material behavior, as well as finite element unit cell and a new semianalytical elastoplastic shear lag model in the inelastic range. GMC is shown to definitely have a window of applicability when simulating discontinuously reinforced composite material behavior.

Pahr, D. H.↗
NASA NTRS1990

Theoretical line shapes for rotational spectra of HCl in Ar

An interaction potential for HCl-Ar recently derived from spectra of van der Waals complexes (Hutson,1988) has been used to obtain converged close-coupling scattering S matrices from which pressure broadening (linewidth, line shift, and line coupling) cross sections have been computed within the Fano-Ben-Reuven formalism. Linewidths agree well with experimental data with the exception of the 0-1 line at low temperatures (an error of about 15 percent at 125 K). Line shifts, especially the largest, for the 0-1 line, are also in reasonable accord with experimental data. Line coupling is predicted to be quite small. The coupled-states molecular scattering approximation is found to be accurate except at very low collision energies. The infinite-order sudden approximations, on the other hand, is not reliable for this system. Inelastic scattering out of the spectroscopic levels accounts for only a fraction of the linewidth at thermal energies, especially in the lower rotational levels.

Green, Sheldon↗
NASA NTRS2001

The Applicability of the Generalized Method of Cells for Analyzing Discontinuously Reinforced Composites

The paper begins with a short overview of the recent work done in the field of discontinuous reinforced composites, focusing on the different parameters which influence the material behavior of discontinuous reinforced composites, as well as the various analysis approaches undertaken. Based on this overview it became evident, that in order to investigate the enumerated effects in an efficient and comprehensive manner, an alternative approach to the computationally intensive finite-element based micromechanics approach is required. Therefore, an investigation is conducted to demonstrate the utility of utilizing the generalized method of cells (GMC), a semi-analytical micromechanics-based approach, to simulate the elastic and elastoplastic material behavior of aligned short fiber composites. The results are compared with (1) simulations using other micromechanical based mean field models and finite element (FE) unit cell models found in the literature given elastic material behavior, as well as (2) finite element unit cell and a new semianalytical elastoplastic shear lag model in the inelastic range. GMC is shown to definitely have a window of applicability when simulating discontinuously reinforced composite material behavior.

Pahr, D. H.↗
DOE OSTI2020

Temperature-dependent phonon lifetimes and thermal conductivity of silicon by inelastic neutron scattering and ab initio calculations

In this paper, inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\overset{→}{Q}$ and frequencies ω in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

SHIVER - Spectroscopy HIstogram Visualizer for Event Reduction

Visualizing data from neutron scattering experiments is the first step in understanding the physics. The program is intended to generate and plot cuts and slices, through the four dimensional single crystal inelastic datasets, measured on direct geometry neutron spectrometers at the Spallation Neutron Source (ARCS, CNCS, HYSPEC, SEQUOIA).

Savici, AndreiT [Oak Ridge National Laboratory (OR↗
DOE OSTI2023

Validating ENDF/B-VIII.1beta1 with LLNL Pulsed-Sphere Neutron-Leakage Spectra [Slides]

In this presentation 6 Li Distinct changes are observed. 235 U Discussion shows Improvements (seen in VIII.1beta0) are coming from new PFNS informed by Chi-Nu experimental data. Additionally, 239 Pu large changes coming from inelastic scattering. 9 Be ENDF/B-VIII.0 shows to be better right after elastic peak while 16 O had slight changes. 19 F had small improvements in ENDF/B-VIII.1 and 56 Fe had noticeable changes in thicker spheres. Examples of negligible changes for the sake of completeness are given.

16-Oxygen↗
NASA NTRS2011

PING Gamma Ray and Neutron Measurements of a Meter-Sized Carbonaceous Asteroid Analog

Determining the elemental composition of carbonaceous (spectral type C) asteroids is still one of the basic problems when studying these objects. The only main source of elemental composition information for asteroids is from their optical, NIR and IR properties, which include their spectral reflectance characteristics, albedo, polarization, and the comparison of optical spectroscopy with meteorite groups corresponding to asteroids of every spectral type. Unfortunately, these sources reflect observations from widely contrasting spatial scales that presently yield a void in the continuum of microscopic and macroscopic evidence, a lack of in situ measurement confirmation, and require deeper sensing techniques to discern the nature of these asteroids. The Probing In situ with Neutrons and Gamma rays (PING) instrument is ideally suited to address this problem because it can be used to determine the bulk elemental composition, H and C content, the average atomic weight and density of the surface and subsurface layers of C-type asteroids, and can provide measurements used to determine the difference between and distinguish between different types of asteroids. We are currently developing the PING instrument that combines gamma ray and neutron detectors with a 14 Me V pulsed neutron generator to determine the in-situ bulk elemental abundances and geochemistry of C-type asteroids with a spatial resolution of 1 m down to depths of tens of cm to 1 m. One aspect of the current work includes experimentally testing and optimizing PING on a known meter-sized Columbia River basalt C-type asteroid analog sample that has a similar composition and the same neutron response as that of a C-type asteroid. An important part of this effort focuses on utilizing timing measurements to isolate gamma rays produced by neutron inelastic scattering, neutron capture and delayed activation processes. Separating the gamma ray spectra by nuclear processes results in higher precision and sensitivity elemental composition measurements. Using gated data acquisition techniques allows for the unambiguous identification of gamma ray lines from different isotopes and nuclear processes, especially in situations when limited detector resolution results in overlapping gamma ray lines that cannot be individually resolved. In this paper, we will present the PING basalt layering experimental data, taken at the test facility at NASA Goddard Space Flight Center using the time tagged event-by-event data analysis technique, compared to our MCNPX computer simulation results for the C-type asteroid and basalt layering simulant models. Comparison of these data will show the advantages, validity, and measurement sensitivity of PING's nuclear interrogation methods to obtain more precise and sensitive in situ bulk elemental composition and density measurements of the subsurface of asteroids.

Bodnarik, J.↗
DOE OSTI2021

Application of Machine Learning Algorithms to Identify Problematic Nuclear Data

In this work we aim to show that Machine learning algorithms are promising tools for the identification of nuclear data that contribute to increased errors in transport simulations. We demonstrate this through an application of a machine learning algorithm (Random Forest) to the Whisper/MCNP6 criticality validation library to identify nuclear data that are associated with an increase of the bias (simulated - experimental $k_{eff}$) in the calculations. Specifically, the $k_{eff}$ sensitivity profiles (w.r.t. nuclear data) of 233 U solution benchmarks are used to predict the bias and Shapley Additive Explanations (SHAP) are used to explain how the sensitivities are related to the predicted bias. The SHAP values can be interpreted as sensitivity coefficients of the machine learning model to the $k_{eff}$ sensitivities which are used to make predictions of bias. Using the SHAP values we can identify specific subsets of nuclear data which have the highest probability of influencing bias. We demonstrate the utility of this method by showing how SHAP values were used to identify an inconsistency in the 19 F inelastic scattering nuclear data. The methodology presented here is not limited to transport problems and can be applied to other simulations if there are experimental measurements to compare against, simulations of those experimental measurements, and the ability to calculate sensitivities of the model output with respect to the data inputs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations

Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.

42 ENGINEERING↗
DOE OSTI2020

Understanding dynamic properties of materials using neutron spectroscopy and atomistic simulation

Recent developments in first-principles lattice dynamics and classical force field based molecular dynamics are revolutionising the field of neutron spectroscopy. Herein we present a short review of these methods, their critical role in the supporting of cutting-edge experiments, and how they are improved by matching experimental data. We begin with a brief overview of how lattice dynamics calculations can be compared to inelastic neutron scattering (INS) and molecular dynamics simulations to both INS and quasi-elastic neutron scattering (QENS). We then provide a series of exemplar applications where lattice dynamics and molecular dynamics have been used in conjunction with neutron spectroscopy to bring significant understanding to topical areas of materials science namely: (i) lattice dynamics and INS for the study of hybrid organic-inorganic perovskites (ii) lattice dynamics and INS for the study of flexible porous solids and (iii) molecular dynamics and QENS for probing molecular behaviour in zeolite catalysis. In all three cases, the understanding gained through the synergy of experiment and computation would have been significantly reduced using either in isolation. Finally, we consider the current state of the art, describing outstanding challenges and suggesting future directions in this exciting and fertile area of physical science.

36 MATERIALS SCIENCE↗