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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Equations of motion of slung-load systems, including multilift systems

General simulation equations are derived for the rigid body motion of slung-load systems. This work is motivated by an interest in trajectory control for slung loads carried by two or more helicopters. An approximation of these systems consists of several rigid bodies connected by straight-line cables or links. The suspension can be assumed elastic or inelastic. Equations for the general system are obtained from the Newton-Euler rigid-body equations with the introduction of generalized velocity coordinates. Three forms are obtained: two generalize previous case-specific results for single-helicopter systems with elastic and inelastic suspensions, respectively; and the third is a new formulation for inelastic suspensions. The latter is derived from the elastic suspension equations by choosing the generalized coordinates so that motion induced by cable stretching is separated from motion with invariant cable lengths, and by then nulling the stretching coordinates to get a relation for the suspension forces. The result is computationally more efficient than the conventional formulation, is readily integrated with the elastic suspension formulation, and is easily applied to the complex dual-lift and multilift systems. Results are given for two-helicopter systems; three configurations are included and these can be integrated in a single simulation. Equations are also given for some single-helicopter systems, for comparison with the previous literature, and for a multilift system. Equations for degenerate-body approximations (point masses, rigid rods) are also formulated and results are given for dual-lift and multilift systems. Finally, linearlized equations of motion are given for general slung-load systems are presented along with results for the two-helicopter system with a spreader bar.

Cicolani, Luigi S.↗

NASA NTRS2002

Incorporation of Mean Stress Effects into the Micromechanical Analysis of the High Strain Rate Response of Polymer Matrix Composites

The results presented here are part of an ongoing research program, to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. A micromechanics approach is employed in this work, in which state variable constitutive equations originally developed for metals have been modified to model the deformation of the polymer matrix, and a strength of materials based micromechanics method is used to predict the effective response of the composite. In the analysis of the inelastic deformation of the polymer matrix, the definitions of the effective stress and effective inelastic strain have been modified in order to account for the effect of hydrostatic stresses, which are significant in polymers. Two representative polymers, a toughened epoxy and a brittle epoxy, are characterized through the use of data from tensile and shear tests across a variety of strain rates. Results computed by using the developed constitutive equations correlate well with data generated via experiments. The procedure used to incorporate the constitutive equations within a micromechanics method is presented, and sample calculations of the deformation response of a composite for various fiber orientations and strain rates are discussed.

Goldberg, Robert K.↗

DOE OSTI2025

Libration of hydroxyl groups in layered aluminum (oxy)hydroxides and other material analogs: insights from inelastic neutron scattering and theory

We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.

Wang, Hsiu-Wen [Oak Ridge National Laboratory (ORN↗

DOE OSTI2023

Improving dynamic collision frequencies: Impacts on dynamic structure factors and stopping powers in warm dense matter

Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, in this work, we systematically investigate variations in dynamic electron–ion collision frequencies ν(ω) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν(ω). These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1979

Transport properties in the atmosphere of Jupiter

Activities reported include: (1) testing of the computer program used to obtain transport properties for the Hulburt-Hirschfelder potential; (2) calculation of transport properties for the C2-C interaction; (3) preliminary calculations for the C2-C2 interaction; (4) calculation of transport properties for the C2H-He interaction; (5) consideration of the effect of inelastic collisions on the transport properties; and (6) the use of the Hulburt-Hirschfelder potential to model ion-atom interactions.

Biolsi, L., Jr.↗

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 a 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 NTRS1996

Optimization of residual stresses in MMC's through the variation of interfacial layer architectures and processing parameters

The objective of this work was the development of efficient, user-friendly computer codes for optimizing fabrication-induced residual stresses in metal matrix composites through the use of homogeneous and heterogeneous interfacial layer architectures and processing parameter variation. To satisfy this objective, three major computer codes have been developed and delivered to the NASA-Lewis Research Center, namely MCCM, OPTCOMP, and OPTCOMP2. MCCM is a general research-oriented code for investigating the effects of microstructural details, such as layered morphology of SCS-6 SiC fibers and multiple homogeneous interfacial layers, on the inelastic response of unidirectional metal matrix composites under axisymmetric thermomechanical loading. OPTCOMP and OPTCOMP2 combine the major analysis module resident in MCCM with a commercially-available optimization algorithm and are driven by user-friendly interfaces which facilitate input data construction and program execution. OPTCOMP enables the user to identify those dimensions, geometric arrangements and thermoelastoplastic properties of homogeneous interfacial layers that minimize thermal residual stresses for the specified set of constraints. OPTCOMP2 provides additional flexibility in the residual stress optimization through variation of the processing parameters (time, temperature, external pressure and axial load) as well as the microstructure of the interfacial region which is treated as a heterogeneous two-phase composite. Overviews of the capabilities of these codes are provided together with a summary of results that addresses the effects of various microstructural details of the fiber, interfacial layers and matrix region on the optimization of fabrication-induced residual stresses in metal matrix composites.

Pindera, Marek-Jerzy↗

NASA NTRS2000

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗

DOE OSTI2021

Neutron thermalization in nuclear graphite: A modern story of a classic moderator

Graphite has been of great interest since the dawn of the nuclear era. In this work, we will focus on the thermal neutron scattering of two types of nuclear graphite that have different microstructures, namely, the historic PGA which has a very anisotropic medium-grained needle coke, and the new G347A which is a near-isotropic fine-grained graphite. We report the results of room temperature inelastic neutron scattering measurements for the scattering function, double differential scattering cross-section, and phonon density of states. The obtained results are discussed in conjunction with results extracted from first-principles lattice dynamics simulations of crystalline graphite, using the incoherent approximation and by also including the coherent-one phonon contribution. Excellent agreements have been found between the measured and calculated scattering functions and phonon densities of states. The inclusion of the coherent-one phonon scattering shows a significant improvement of the calculated scattering functions, double differentials, and inelastic scattering cross-sections. We also highlight the physical concerns about the current thermal neutron scattering ENDF/B-VIII of nuclear graphite library and its corresponding phonon density of states. This work aims to provide a new presentation of neutron thermalization by graphite and offers a new set of high quality experimental and computational data for consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1980

Vibration-rotation line shifts for 1 sigma g + H2/V,J/-1S/0/ He computed via close coupling - Temperature dependence

The density shifting of vibration-rotation transitions of H2 perturbed by He was computed (as a function of temperature) with no adjustable parameters. The calculation was carried out using the framework of the impact theory of Baranger with S-matrix elements obtained via close coupling calculations which incorporated the ab initio H2-H2 system potential of Tsapline et al.(1977). Vibrational and rotational inelasticity were neglected in the calculations; nevertheless good agreement with experimental data was obtained, up to moderate temperatures, for the density shift. A much poorer comparison was obtained for the density broadening.

Hahne, G. E.↗

NASA NTRS1988

A crystallographic model for nickel base single crystal alloys

The purpose of this research is to develop a tool for the mechanical analysis of nickel-base single-crystal superalloys, specifically Rene N4, used in gas turbine engine components. This objective is achieved by developing a rate-dependent anisotropic constitutive model and implementing it in a nonlinear three-dimensional finite-element code. The constitutive model is developed from metallurgical concepts utilizing a crystallographic approach. An extension of Schmid's law is combined with the Bodner-Partom equations to model the inelastic tension/compression asymmetry and orientation-dependence in octahedral slip. Schmid's law is used to approximate the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response and strain-rate sensitivity of the single-crystal superalloys. Methods for deriving the material constants from standard tests are also discussed. The model is implemented in a finite-element code, and the computed and experimental results are compared for several orientations and loading conditions.

Dame, L. T.↗

DOE OSTI2022

Validating first-principles phonon lifetimes via inelastic neutron scattering

Phonon lifetimes are a key component of quasiparticle theories of transport; yet first-principles lifetimes are rarely directly compared with inelastic neutron scattering (INS) results. Existing comparisons show discrepancies even at temperatures where perturbation theory is expected to be reliable. In this paper, we demonstrate that the reciprocal space voxel (q voxel), which is the finite region in reciprocal space required in INS data analysis, must be explicitly accounted for within theory in order to draw a meaningful comparison. Here, we demonstrate accurate predictions of peak widths of the scattering function when accounting for the q voxel in CaF 2 and ThO 2 . Passing this test implies high fidelity of the phonon interactions and the approximations used to compute the Green's function, serving as a critical benchmark of theory and indicating that other material properties should be accurately predicted, which we demonstrate for thermal conductivity.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Validity of approximate methods in molecular scattering - Thermal HCl-He collisions

Accurate close coupling scattering calculations are presented for thermal energy HCl-He collisions. The interaction potential is obtained from the Gordon-Kim electron gas model, adjusted to have the correct long-range multipole form. A variety of phenomenological cross sections are computed from the close coupling S matrix, and these are compared with results from several commonly employed approximate methods. In particular, it is found that the total integral, total differential, and gas kinetic cross sections are accurately predicted by the central field approximation which retains just the spherical average of the interaction. Integral inelastic cross sections are represented quite accurately by the coupled states approximation of McGuire and Kouri, but only qualitatively by the effective potential method of Rabitz.

Green, S.↗

NASA NTRS1996

Absolute Inelastic Differential Electron Scattering Cross Sections for He, Xe, N2, and CO at Near Threshold Impact Energies for 90(deg) Scattering Angle

There is a great deal of need for accurate differential cross sections (DCSs) associated with electron impact excitation of various atomic and molecular species. These cross sections are needed for modeling of various plasma systems (ranging from lasers and material processing plasmas to planetary and astrophysical plasmas) and for guiding development of theoretical and computational schemes.

DCS beam scattering beam-beam scattering↗

NASA NTRS1995

Shielding materials for highly penetrating space radiations

Interplanetary travel involves the transfer from an Earth orbit to a solar orbit. Once outside the Earth's magnetosphere, the major sources of particulate radiation are solar cosmic rays (SCR's) and galactic cosmic rays (GCR's). Intense fluxes of SCR's come from solar flares and consist primarily of protons with energies up to 1 GeV. The GCR consists of a low flux of nuclei with energies up to 10(exp 10) GeV. About 70 percent of the GCR are protons, but a small amount (0.6 percent) are nuclei with atomic numbers greater than 10. High energy charged particles (HZE) interact with matter by transferring energy to atomic electrons in a Coulomb process and by reacting with an atomic nucleus. Energy transferred in the first process increases with the square of the atomic number, so particles with high atomic numbers would be expected to lose large amounts of energy by this process. Nuclear reactions produced by (HZE) particles produce high-energy secondary particles which in turn lose energy to the material. The HZE nuclei are a major concern for radiation protection of humans during interplanetary missions because of the very high specific ionization of both primary and secondary particles. Computer codes have been developed to calculate the deposition of energy by very energetic charged particles in various materials. Calculations show that there is a significant buildup of secondary particles from nuclear fragmentation and Coulomb dissociation processes. A large portion of these particles are neutrons. Since neutrons carry no charge, they only lose energy by collision or reaction with a nucleus. Neutrons with high energies transfer large amounts of energy by inelastic collisions with nuclei. However, as the neutron energy decreases, elastic collisions become much more effective for energy loss. The lighter the nucleus, the greater the fraction of the neutron's kinetic energy that can be lost in an elastic collision. Thus, hydrogen-containing materials such as polymers are most effective in reducing the energy of neutrons. Once neutrons are reduced to very low energies, the probability for undergoing a reaction with a nucleus (the cross section) becomes very high. The product of such a reaction is often radioactive and can involve the release of a significant amount of energy. Thus, it is important to provide protection from low energy neutrons during a long duration space flight. Among the light elements, lithium and boron each have an isotope with a large thermal neutron capture cross section, Li-6 and B-10. However, B-10 is more abundant in the naturally-occurring element than Li-6, has a thermal neutron capture cross section four times that of Li-6, and produces the stable products, He-4 and Li-7 in the interaction while Li-6 produces radioactive tritium (H-3). Thus, boron is the best light-weight material for thermal neutron absorption in spacecraft. The work on this project was focused in two areas: computer design where existing computer codes were used, and in some cases modified, to calculate the propagation and interactions of high energy charged particles through various media, and materials development where boron was incorporated into high performance materials.

Kiefer, Richard L.↗

DOE OSTI2022

Analysis of Microwave Resonant Cavity Transducer Degradation Mechanisms

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Stress Intensity Thresholds for Development of Reliable Brittle Materials

Brittle material failure in high consequence systems can appear random and unpredictable at subcritical stresses. Gaps in our understanding of how structural flaws and environmental factors (humidity, temperature) impact fracture propagation need to be addressed to circumvent this issue. A combined experimental and computational approach composed of molecular dynamics (MD) simulations, numerical modeling, and atomic force microscopy (AFM) has been undertaken to identify mechanisms of slow crack growth in silicate glasses. AFM characterization of crack growth as slow as 10 -13 m/s was observed, with some stepwise crack growth. MD simulations have identified the critical role of inelastic relaxation in crack propagation, including evolution of the structure during relaxation. A numerical model for the existence of a stress intensity threshold, a stress intensity below which a fracture will not propagate, was developed. This transferrable model for predicting slow crack growth is being incorporated into mission-based programs.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Second Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. Because the transducer is intended for immersion in a high temperature corrosive fluid, understanding of material degradation is crucial for estimation of transducer performance lifetime, and development of measurement interpretation algorithms. We conducted a preliminary computational investigation of relevant damage mechanisms of a stainless steel 316 cylindrical resonator in FLiBe salt. The two main damage mechanisms, creep and corrosion, were modeled using multiphysics COMSOL software. Degradation was modeled for a temperature range 500°C to 700°C. Coupling of the damage mechanisms was not considered. These models predict significant inelastic deformation at most temperatures due to creep, and qualitatively predict chromium depletion both along the liquid/solid interface and along the grain boundaries. An algorithmic approach for compensation of these degradation effects during fluid flow measurements will be developed in the future work. To validate sensor physics, we have performed proof-of-principle test of flow sensing in water. For this test, we have developed a cylindrical resonator for K-band, which was machined from brass. The cavity was excited through WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder. To increase the spectral signal visibility, we developed a signal processing method for baseline subtraction. A flow loop for proof-of-principle test of transducer performance in water was assembled. A commercial flow meter was installed in the loop for reference measurements. Cylindrical cavity was excited in the TE 011 mode with resonant frequency f ≈ 17.8GHz. Frequency shift of cavity spectral response was obtained by gradually increasing water flow rate from 0 to 60gpm. Corresponding monotonic increase of resonant frequency shift by several MHz was observed. Approximate figure of merit of sensitivity to flow rate is 100KHz/GPM. In addition, we have identified an existing liquid sodium experimental setup for demonstration of flow sensing in environment similar to that of an advanced reactor. The setup consists of a cylindrical vessel and center feed line, where transducer inserted through the lid will measure velocity of the impinging liquid jet. As a calibration experiment, we have assembled a water vessel with center feed with the same dimensions as those of the liquid sodium setup. We have also developed and insertion probe consisting of a 50cm brass waveguide enclosed in protective SS316 tube. Using the water loop, we have demonstrated feasibility of sensing the impinging liquid jet in the vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗