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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

One Dimensional Penstock Flow Models for Hydropower Digital Twin

This report details the inelastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The flow model is coupled to a mechanistic turbine model to compute the rotation rate of the turbine, which is dependent on the flow rate, gravitational head, and various parameters of the machine, such as inlet vane angle. Results are given for an example hydropower facility in which the turbine is subjected to a specified time-varying load, and the results on the flow rate are shown. Future applications for the flow model are also briefly reviewed.

13 HYDRO ENERGY↗
DOE OSTI2023

Elastic nucleon-pion scattering at m π = 200 MeV from lattice QCD

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π = 200 MeV. The s-wave scattering lengths with both total isospins I = 1/2 and I = 3/2 are inferred from the finite-volume spectrum below the inelastic threshold together with the I = 3/2 p-wave containing the Δ(1232) resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels, enabling a determination of the I = 3/2 scattering length with statistical errors below 5%, while the I = 1/2 scattering length is somewhat less precisely evaluated. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a step toward future computations down to physical light quark masses with multiple lattice spacings and volumes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Anharmonic Origin of the Giant Thermal Expansion of NaBr

ll phonons in a single crystal of NaBr are measured by inelastic neutron scattering at temperatures of 10, 300, and 700 K. Even at 300 K, the phonons, especially the longitudinal-optical phonons, show large shifts in frequencies and show large broadenings in energy owing to anharmonicity. Ab initio computations are first performed with the quasiharmonic approximation (QHA) in which the phonon frequencies depend only on V and on T only insofar as it alters V by thermal expansion. This QHA is an unqualified failure for predicting the temperature dependence of phonon frequencies, even 300 K, and the thermal expansion is in error by a factor of 4. Ab initio computations that include both anharmonicity and quasiharmonicity successfully predict both the temperature dependence of phonons and the large thermal expansion of NaBr. The frequencies of longitudinal-optical phonon modes decrease significantly with temperature owing to the real part of the phonon self-energy from explicit anharmonicity originating from the cubic anharmonicity of nearest-neighbor NaBr bonds. Anharmonicity is not a correction to the QHA predictions of thermal expansion and thermal phonon shifts but dominates the behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Boron-epoxy-reinforced titanium aircraft landing-gear drag strut

The structural performance of a boron-epoxy-reinforced titanium drag strut, containing a bonded scarf joint and designed to the criteria of a large commercial transport, has been evaluated experimentally and analytically. The strut was exposed to two lifetimes of fatigue loading and was statically loaded to the tensile and compressive design ultimate loads. Throughout the test program no evidence of any damage in the drag strut was detected by strain-gage measurements, ultrasonic inspection, or visual observation. The bonded joint was analyzed using the NASTRAN computer program. A comparison of the strains predicted by the NASTRAN computer program with the experimentally determined values shows excellent agreement. An analytical study indicated that the nonlinear behavior of a structural spacer at each end of the strut could be explained by the inelastic behavior and possible creep of the adhesive.

Howell, W. E.↗
NASA NTRS1990

Equations of motion of slung load systems with results for dual lift

General simulation equations are derived for the rigid body motion of slung load systems. These systems are viewed as consisting of several rigid bodies connected by straight-line cables or links. The suspension can be assumed to be elastic or inelastic, both cases being of interest in simulation and control studies. Equations for the general system are obtained via D'Alembert's principle and the introduction of generalized velocity coordinates. Three forms are obtained. Two of these generalize previous case-specific results for single helicopter systems with elastic or inelastic suspensions. The third is a new formulation for inelastic suspensions. It is derived from the elastic suspension equations by choosing the generalized coordinates so as to separate motion due to cable stretching from motion with invariant cable lengths. The result is computationally more efficient than the conventional formulation, and is readily integrated with the elastic suspension formulation and readily applied to the complex dual lift and multilift systems. Equations are derived for dual lift systems. Three proposed suspension arrangements can be integrated in a single equation set. The equations are given in terms of the natural vectors and matrices of three-dimensional rigid body mechanics and are tractable for both analysis and programming.

Cicolani, Luigi S.↗
NASA NTRS2001

General Multimechanism Reversible-Irreversible Time-Dependent Constitutive Deformation Model Being Developed

Since most advanced material systems (for example metallic-, polymer-, and ceramic-based systems) being currently researched and evaluated are for high-temperature airframe and propulsion system applications, the required constitutive models must account for both reversible and irreversible time-dependent deformations. Furthermore, since an integral part of continuum-based computational methodologies (be they microscale- or macroscale-based) is an accurate and computationally efficient constitutive model to describe the deformation behavior of the materials of interest, extensive research efforts have been made over the years on the phenomenological representations of constitutive material behavior in the inelastic analysis of structures. From a more recent and comprehensive perspective, the NASA Glenn Research Center in conjunction with the University of Akron has emphasized concurrently addressing three important and related areas: that is, 1) Mathematical formulation; 2) Algorithmic developments for updating (integrating) the external (e.g., stress) and internal state variables; 3) Parameter estimation for characterizing the model. This concurrent perspective to constitutive modeling has enabled the overcoming of the two major obstacles to fully utilizing these sophisticated time-dependent (hereditary) constitutive models in practical engineering analysis. These obstacles are: 1) Lack of efficient and robust integration algorithms; 2) Difficulties associated with characterizing the large number of required material parameters, particularly when many of these parameters lack obvious or direct physical interpretations.

Saleeb, A. F.↗
DOE OSTI2021

Magnetic field induced phase transitions of the triangular spin model

In this work, motivated by a recent inelastic neutron scattering experiment on YbMgGaO 4 , I study the homogeneous spin model on the triangular lattice using the cluster mean-field theory and exact diagonalization. Phase diagrams with and without magnetic field are computed. The results show that the magnetic field can induce the phase transition from the spin liquid state (or the stripe state) to the 120 ° antiferromagnetic state. These phase transitions are suppressed by the next-nearest-neighbor exchange interaction J 2 / J 1 and vanish as J 2 / J 1 > 0.11 . I analyze a parameter space at J 2 / J 1 = 0.1 , where a field induced phase transition can occur, and find that the deviation of the theoretical spin excitation energy from the experimental data is only about 0.054. These results imply that an effective homogeneous spin model still applies to YbMgGaO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Criticality Experiments to Reduce Compensating Errors in Plutonium Nuclear Data

Compensating errors between nuclear data observables in a library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors in nuclear data. A new criticality experiment, described in this work, was designed with the specific target nuclear data of 239 Pu fission, inelastic scattering, elastic scattering, capture, nu-bar, and prompt fission neutron spectrum (PFNS). This work will focus on the design and execution of the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). The criticality experiment includes two different configurations with very different geometries: one is cube-like to minimize neutron leakage while the other is slab-like to maximize leakage. Having these two widely varying configurations allows the scattering sensitivities of 239 Pu to the neutron multiplication factor to be greatly changed while minimally impacting the other cross section sensitivities. Both configurations utilize the Pu ZPPR (Zero Power Physics Reactor) plates as fuel. The experiments were designed using a D-Optimality criteria, which is an optimization method minimizing the log-determinant of the adjusted nuclear data covariance for the target reactions. These experiments include not only inference of k eff , as done in all critical benchmark experiments, but several other responses as well, such as neutron multiplication measurements and reaction rate ratios. After analysis of the measured data is complete, adjustment of nuclear data will be performed to assess whether the new experimental data successfully reduced compensating errors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Neural network ensembles and uncertainty estimation for predictions of inelastic mechanical deformation using a finite element method-neural network approach

The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.

42 ENGINEERING↗
DOE OSTI2020

Wave optics simulations for new soft x-ray beamlines at NSLS-II

In this work, we performed fully- and partially-coherent synchrotron emission and propagation simulations with the "Synchrotron Radiation Workshop" computer code to analyze the performance of two soft X-ray beamlines under development at the National Synchrotron Light Source II: Soft X-ray Nanoprobe (SXN), and Angle-Resolved Photoemission Spectroscopy (ARPES) and Resonant Inelastic X-ray Scattering (RIXS) Imaging (ARI). The SXN beamline intends to provide high flux and high spatial resolution coherent soft X-ray imaging capabilities using both zone plate and lensless coherent imaging techniques. The ARI beamline aims to perform high flux ARPES and RIXS experiments with a focal spot size at the sample approaching 100 nm using highly-demagnifying mirrors in Kirkpatrick-Baez geometry. To accurately calculate the resolution and the degree of X-ray coherence provided by the two state-of-the-art beamlines, partial coherence effects are required to be taken into account in wave optics simulations for these two beamlines. In this talk, beamline performance parameters such as spot size, degree of coherence, flux, and energy resolution at the sample are presented. The effects of mirror surface slope errors on beamline performance were studied and some suggestions for further optimization are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Charge carrier dynamics in two-dimensional hybrid perovskites: Dion–Jacobson vs. Ruddlesden–Popper phases

Two-dimensional (2D) halide perovskites are promising materials for environmentally stable next-generation optoelectronic device applications. Besides the widely investigated Ruddlesden–Popper (RP) phases, Dion–Jacobson (DJ) phases are attracting considerable attention due to their rapid emergence as efficient solar cell materials. However, there is very little atomistic understanding of the charge carrier dynamics under ambient conditions for these DJ-phases, limiting the possibilities to tune their optoelectronic performances through compositional engineering routes. In this study, by combining nonadiabatic molecular dynamics with time-domain density functional theory methods at room temperature, we compare the dominant non-radiative carrier recombination and dephasing processes in RP and DJ monolayered lead halide perovskites. Our systematic study demonstrates that performance-limiting nonradiative carrier recombination processes greatly depends on the electron–phonon interactions induced by structural fluctuations and instantaneous charge localization in these materials. The stiffer interlayer packing due to the presence of single spacer dications, which separates the lead iodide slabs, reduces the thermal fluctuations in the DJ-phase to a greater extent than that in the RP-phase 2D-perovskites. Specific electronic coupling between the closely spaced lead iodide layers enhances the delocalization of band-edge charge densities in DJ-phase systems. Compared to the RP-phase, reduced inelastic electron–phonon scattering in DJ-phase perovskites significantly limits intrinsic non-radiative recombination processes. The consequent enhancement in the photogenerated charge carrier lifetime makes DJ-phase perovskites potentially suitable for various optoelectronic devices. The computational insights gained from this study allow us to outline a set of robust design principles for DJ-phase perovskites to strategically tune their optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Nonlinear structural and life analyses of a combustor liner

Three dimensional, nonlinear finite element structural analyses were performed for a simulated combustor liner specimen to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic creep response of the one half scale component. Results indicate continued cyclic hardening and ratcheting while experimental data suggested a stable stress strain response after only a few loading cycles. The computed stress strain history at the critical location was put into two life prediction methods, strainrange partitioning and a Pratt and Whitney combustor life prediction method to evaluate their ability to predict cyclic crack initiation. It is found that the life prediction analyses over predicted the observed cyclic crack initiation life.

Moreno, V.↗
NASA NTRS1986

Constitutive modeling of cyclic plasticity and creep, using an internal time concept

Using the concept of an internal time as related to plastic strains, a differential stress-strain relation for elastoplasticity is rederived, such that (1) the concept of a yield-surface is retained; (2) the definitions of elastic and plastic processes are analogous to those in classical plasticity theory; and (3) its computational implementation, via a 'tangent-stiffness' finite element method and a 'generalized-midpoint-radial-return' stress-integration algorithm, is simple and efficient. Also, using the concept of an internal time, as related to both the inelastic strains as well as the Newtonian time, a constitutive model for creep-plasticity interaction, is discussed. The problem of modeling experimental data for plasticity and creep, by the present analytical relations, as accurately as desired, is discussed. Numerical examples which illustrate the validity of the present relations are presented for the cases of cyclic plasticity and creep.

Watanabe, O.↗
DOE OSTI2021

Investigation into the transport properties of planetary interiors through inelastic X-ray scattering experiments and quantum molecular dynamics (Final Technical Report)

This was a two-year research project for the period 08/15/2018 - 08/14/2020 conducted at the University of Nevada, Reno using Pronghorn, a newly built high-performance cluster located at the University. The short-term project goal was to provide computational support for the experimental campaigns at LCLS through data analysis of previous and upcoming scheduled experiments, along with performing state-of-the-art atomistic simulations of dense plasmas. The longer-term goal was to create a new computational high energy density physics group located at the University of Nevada, Reno. Guided by results from past and future high-resolution scattering experiments, we aimed to research and develop non-equilibrium and non-adiabatic atomistic simulations of dense plasmas that went beyond the Born-Oppenheimer approximation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structure and spectroscopy of graphite monofluoride

The structure of graphite monofluoride, (CF) n , has been debated since its discovery in 1934. In this work, we investigate a commercial graphite monofluoride by vibrational spectroscopy (infrared, Raman and the first inelastic neutron scattering spectra of this material). The spectroscopy shows that the material contains unreacted graphite and the partially fluorinated product dicarbon fluoride, (C 2 F) n , We evaluate the previously proposed $P\bar{6}m2$ and $P\bar{3}m1$ structures using computational methods and find F···F contacts render the $P\bar{6}m2$ structure dynamically unstable. We propose two alternative structures, $Cmc2_1$ $P6_3mc$, generated by displacement of one layer relative to another and find that $Cmc2_1$ is also dynamically unstable. The calculations are validated by comparison of calculated and observed INS spectra

Density functional theory↗
DOE OSTI2022

State-to-state rate coefficients for HCS+ in rotationally inelastic collisions with H2 at low temperatures

ABSTRACT HCS+ ions have been detected in several regions of the interstellar medium (ISM), but an accurate determination of the chemical-physical conditions in the molecular clouds where this molecule is observed requires detailed knowledge of the collisional rate coefficients with the most common colliders in those environments. In this work, we study the dynamics of rotationally inelastic collisions of HCS+ + H2 at low temperature, and report, for the first time, a set of rate coefficients for this system. We used a recently developed potential energy surface for the HCS+–H2 van der Waals complex and computed state-to-state rotational rate coefficients for the lower rotational states of HCS+ in collision with both para- and ortho-H2, analysing the influence of the computed rate coefficients on the determination of critical densities. Additionally, the computed rate coefficients are compared with those obtained by scaling the ones from HCS+ in collision with He (an approximation that is sometimes used when data is lacking), and large differences are found. Furthermore, the approximation of using the rates for the HCO+ + H2 collision as a rough approximation for those of the HCS+ + H2 system is also evaluated. Finally, the complete set of de-excitation rate coefficients for the lowest 30 rotational states of HCS+ by collision with H2 is reported from 5 to 100 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Computational investigation of hydrogen-induced phonon changes in carbon fiber

Optical vibrational spectroscopy has shown promise as a noninvasive means of monitoring the mechanical properties of carbon fiber (CF), which is increasingly used for industrial and consumer purposes. However, interpretation of optical vibrational spectra for solid materials is inferential, particularly when defects are present. Because inelastic neutron scattering (INS) spectroscopy is not subject to selection rules, the full vibrational spectra can be measured. And, identifying correlations between INS features and tensile properties can assist in the interpretation of spectra from more commonly used optical vibrational spectroscopic techniques, such as Raman and infrared (IR) spectroscopy. Recent INS experiments on high-performance commercial carbon fibers showed features near 900 and 1100 cm –1 in addition to a broad feature near 3000 cm –1 that increased in intensity with decreasing tensile strength. These features were assigned to hydrogen defects. Here, we use density functional theory to simulate the INS spectra of several hydrogen defect geometries in graphite as a model for carbon fiber structure units, confirming the experimental assignment of these peaks to hydrogen modes and providing insights into the structure and lattice dynamics of the defects.

36 MATERIALS SCIENCE↗